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evaporation  anO  the  square  root  of  tlTe.  In  the  second  stud;,  'Sragq' 
soybeans  were  grown  in  330.  660  and  990  ppm  COj  concentration  under 
well  watered  conditions  to  obtain  a data  base  to  evaluate  the  ability 
of  the  model  to  predict  water  use  during  canopy  development.  Average 
dally  error  in  the  predicted  water  use  was  less  than  6t  and  the 
differences  between  the  measured  and  the  simulated  values  were  not 
significant  at  the  95*  level  of  confidence. 

The  model  was  used  to  study  effect  of  increasing  atmospheric  COg 
concentration  from  330  to  990  ppm  ( "Green  House'  effect  ) on  water 
use  of  field  grown  soybeans.  If  the  predicted  increase  In  air 
temperature  was  Ignored,  higher  CO^  may  reduce  seasonal  water  use  by 
10*.  However,  when  the  effect  of  increased  00^  on  stomatal 
resistance,  leaf  area  development  and  a 5 *C  temperature  increase  were 
considered,  the  model  predicted  a L9*  increase  in  water  use  and 
possibly  a 3 °C  rise  in  leaf  temperaure. 

Effects  of  variations  In  weather,  irrigation  frequencies,  and 
planting  dates  on  soybean  crop  coefficients  were  evaluated.  The  study 
showed  that  the  use  of  published  crop  coefficients  values  for  varying 
climate  and  growing  practices  may  rsult  in  over  or  under  prediction  of 
crop  water  requirements  by  as  mucn  as  15  to  20*.  When  better  accuracy 
Is  required,  the  model  can  be  used  to  adjust  crop  coefficients  to  more 
closely  fit  local  weather  and  management  conditions. 


CHAPTER  I 
INTRODUCTION 


trinsp’ration  nay  lead  to  confusion  wten  attempts  are  made  to  estimate 
crop  water  requirements.  Transpiration  1s  a necessary  eipenditure  of 
water  1n  crop  production.  Evaporation,  on  the  other  hand,  is  often  an 
avoidable  water  loss.  An  irrigation  strategy  which  minimites 
evaporation  from  the  soil,  nay  produce  savings  of  water. 

The  separation  of  evapotranspiration  Into  evaporation  and  trans- 
piration and  their  nutual  interrel atlonship  toms  an  Important  problem 
in  water  conservation  research  and  management.  During  the  tine  course 
of  an  annual  crop,  the  crop  system  changes  from  one  in  which 
evapotrahspiratioh  Is  entirely  soil  evaporation  to  one  in  whicn  it  is 
Phostly  transpiration,  and  finally  to  one  in  which  both  transpi ration 
and  evaporation  are  affected  by  crop  senescence  (Jones  et  al., 
1994).  Soil  evaporation  constitutes  a major  loss  of  soil  water  when 
the  soil  goes  through  frequent  cycles  of  wetting  and  drying  and  when 
plant  canopies  are  Incomplete.  Models  are  being  developed  using 
various  approaches  and  levels  of  detail  to  Improve  the  prediction  of 
evapotranspiration  under  such  situations. 

Irrigation  management  requires  short  tern  estimates  of 
evapotranspiration  to  allow  farmers  to  make  decisions  concerning 
timing  and  amount  of  irrigation.  In  one  approach  to  predicting 
evapotranspiration,  the  whole  realm  of  crop  responses  to  the  soi1- 
plant-climate  and  management  factors  Is  accounted  for  by  crop 
coefficients  at  different  stages  of  growth.  The  actual  crop  water 
needs  are  detemined  by  multiplying  reference  crop  evapotranspiration 
by  the  stage  and  crop  specific  coefficient  (Jenson  et  al.,  1970; 
Ooorenbos  anO  Pruitt,  1977).  Recently.  Ritchie  (1972)  developed  an 
enplrica!  model  to  estimate  soli  evaporation  and  transpi ration  from 


well  watered  crops  wliere  the  soil  goes  through  cycles  of  wetting  and 
drying.  In  this  model  transpiration  was  related  to  leaf  area 
deveTopnent  and  potential  evapotranspi ration  using  an  empirical 
function  of  leaf  area  to  account  for  the  degree  of  crop  cover.  The 
soil  evaporation  was  computed  in  two  stages.  Following  a rain  or 
Irrigation  event  and  until  a certain  cumulative  evaporation  occurred, 
evaporation  was  calculated  as  a function  of  net  radiation  reaching  the 
surface.  From  this  time  on,  evaporation  was  determined  using  equations 
developed  by  Black  et  al.  (1969)  relating  cumulative  evaporation  to 


the  square  root  of  time. 

The  evapotranspiration  process  Is  primarily  dependent  upon  incom- 
ing solar  radiation  when  water  is  non-limiting.  In  Che  absence  of 
advective  energy,  evapotranspiration  proceeds  at  a rate  that  is  about 
equal  to  the  evaporation  equivalent  of  Che  net  radiation  level. 
Pronounced  advection  from  wind  In  sem1-arld  and  arid  climate  can  occur 
as  warm  dry  air  mass  passes  over  irrigated  fields  (Doorenbos  and 
Pruitt,  1977),  Therefore,  empirical  methods  oust  be  corrected  to 
improve  their  accuracy  under  local  conditions. 

A promising  alternative  approach  is  where  the  soil-plant- 
atmosphere  is  treated  as  an  energy  and  mass  exchange  system.  This 
approach  takes  into  account  the  energy  sources  and  physical  cha- 
racteristics of  the  soil  and  plant.  The  principle  driving  force  for 
water  vapor  is  the  difference  between  the  vapor  pressure  at  the  evapo- 
rating surface  and  the  bulk  air.  The  resulting  rate  of  water  loss  is 
proportional  to  the  energy  gradient  that  develops  in  the  soil-plaht- 
system  and  is  inversely  proportional  to  flow  resistances  in  the  path- 


flow. 


HontsUh  (1965)  combined  energy  balance,  aerodynamic  and  surface 
resistances  and  developed  a model  to  compute  water  loss  from  a closed 
canopy.  He  computed  the  aerodynamic  resistances  from  wind  speed  and 
canopy  roughness  cnaracteristtcs,  whereas  crop  resistance  was  computed 
by  considering  stomatal  resistance  of  various  leaf  layers  acting  In 
parallel.  Hontelth's  model  could  not  be  extended  to  the  problems  of 
soil  evaporation  because  soil  resistance  estimates  or  measurements 
were  lacking  except  for  situations  of  wet  and  dry  soils.  Waggoner  and 
Relfsnyder  (1965).  Stewart.  1970,  Shawcroft  et  a1.  (1974),  and 
Goudriaan  (1977)  have  built  complex  micrometeorological  models  to 
define  crop  and  soli  responses  to  environment.  The 

micrometeorological  models  are  most  suitable  for  crop  responses  to  its 
environment.  However,  In  these  models  the  soil  was  usually  simplified 
as  fully  wet  (resistance  • 0)  or  dry,  I.e.  very  high  resistance. 
Tnese  assmiptfons  are  not  suitable  In  real  crop  situations  where  the 
soil  surface  goes  through  repeated  cycles  of  wetting  and  drying. 

Therefore  the  objectives  of  this  study  were; 

1.  Develop  and  test  a model  of  son  evaporation  process  that 
takes  into  account  tne  dynamics  of  soil  wetting  and  drying 
including  intermittent  conditions  of  soil  wetness. 

Z.  Develop  an  integrated  evapotranspiration  model,  that  treats 
soil  evaporation  and  transpl ration  as  a separate  processes, 
for  predicting  the  water  use  and  microclimate  of  a well 
watered  developing  crop,  where  the  soil  surface  may  go 
through  cycles  of  wetting  and  drying. 


OeveTap  salutians  for  tne  inecrianlstlc  inodol  that  will  allow 
Its  use  with  normally  measured  climate  and  crop  variables  so 
that  It  can  be  used  for  practical  applications  and  evalua- 

Conduct  experiments  In  outdoor,  controlled  environment 
chambers  to  develop  a data  base  to  evaluate  the  ability  of 
the  model  to  predict  evaporation  and  tnanspiratlon  from  a 
developing  well  watered  soybean  crop  grown  at  330.  660  and 
990  ppm  00^  concentration  environments. 

Apply  the  model  to  study  the  ramification  of  CO^  induced 
changes  In  stomatal  resistance,  leaf  area  development  and 
temperatures  on  water  use  of  field  grown  soybeans. 

Determine  the  effect  of  variations  in  tne  climate,  cultural 
practices  and  irrigation  management  strategies  on  crop  coef- 
ficients for  a soybean  crop  grown  on  sandy  soils. 


CHAPTER  II 

A HECHANlSriC  HOOEL  OF  EVAPORATION  FRW  DRYING  SOIL 
IntroNuctiON 

Evaporation  (E)  anO  transpi ration  (T)  are  physIcallY  alitllar 
processes.  Evaporation  Is  the  conversion  of  water  Into  vapor  and  Its 
transfer  from  a soil  to  tNe  atmosphere,  wfiereas  transpiration  is  evap- 
oration from  a plant  surface.  It  1s  extremely  difficult  to  measure 
evaporation  separately  from  transpiration  in  a soil  In  which  plants 
are  growing.  Therefore,  the  two  processes  are  considered  together 
under  one  name,  evapotranspiratlon,  or  ET, 

Evaporation  varies  In  response  to  soil  water  status,  the  extent 
of  crop  cover  and  the  niicrocl Irnate  of  the  canopy.  For  some  time  after 
plant  emergence,  the  water  use  of  a crop  may  be  dominated  by  the  evap- 
oration from  a bare  soil  surface.  As  the  plant  continues  to  grow, 
transpiration  becomes  more  significant  relative  to  soH  evaporation 
and  for  a period  of  time  the  evapotranspiratlon  Is  related  to  the 
degree  of  crop  cover.  The  separation  of  ET  Into  £ and  T and  their 
mutual  Interrelationship  form  Important  but  unresolved  problems,  espe- 
cially in  water  conservation  research  and  practice.  Hodern  precision 
techniques  of  Innigatlon  application  and  methods  of  evaporation  and 
transpiration  control  now  being  developed,  require  this  Information 
to  allow  their  effect  on  evapotranspiratlon  to  be  evaluated. 

Recently  there  has  been  considerable  progress  In  the  development 
of  meteorological  methods  for  estimating  the  Individual  components  of 
evapotranspiratlon.  Models 


are  being  developed  using 


acHes  and  levels  of  details.  One  of  the  wire  promising  approaches  is 
that  ahere  the  soil -plant-atmosphere  is  Created  as  an  energy  and  mass 
eschange  system.  This  approach  takes  into  account  the  enengy  sources 
and  the  physical  characteristics  of  the  soil,  plant  and  the 
atmosphere.  hotaPle  among  these  models  are  those  developed  by 
Waggoner  and  Relfsoyder  (1968),  Sha»croft  et  al.  (1974).  and  Soudriaan 
(1977)  where  the  soil  as  well  as  the  plant  canopy  was  treated  as  a 
multi-layer  system.  In  these  models,  the  soil  has  been  simplified 
usually  to  be  wet  or  dry  {E^O  ).  These  assumptions  are  not  suitable 
for  real  crop  systetas  as  plants  emerge  from  bare  soH  and  grow  through 
the  season  subjected  to  natural  rainfall  as  well  as  irrigation 
cycles.  A physical  model  of  the  soil  evaporation  process  which  acco- 
unts for  fluctuating  soil  water  conditions  was  needed  so  that  it  can 
be  combined  with  soil-canopy  models  for  predicting  ET. 

Review  of  Literature 
Physics  of  soil  evaporation 

The  evaporation  at  the  surface  is  characterlaed  by  the  mass  and 
the  momentum  balance  of  an  air  layer  above  the  surface.  The  supply  of 
energy,  a favoraOle  vapor  pressure  gradient  and  the  supply  of  water  to 
the  site  of  evaporation  are  the  three  necessary  conditions  for  evapor- 
ation  to  occur.  Energy  is  required  to  supply  the  latent  heat  for 
evaporation,  and  can  be  derived  from  the  evaporating  body  itself,  thus 
causing  it  to  cool,  or  from  outside  as  radiated  or  advective  energy. 
The  vapor  pressure  difference  between  Che  atmosphere  and  the  evaporat- 
ing surface  provides  the  driving  force  for  water  loss.  The  vapor 
pressure  in  the  atmosphere  surrounding  the  evaporating  body  must  rema- 
in lower  than  the  vapor  pressure  at  the  surface  of  that  body,  and  the 
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vapor  tfusc  Pe  transported  away,  by  diffuston  or  convection,  or  botn. 
These  two  conditions  are  influenced  by  meteorological  factors  such  as 
a1r  temperature,  bumldlty,  wind  velocity,  and  radiation  wbicn  together 
constitute  atmospheric  evaporatlvlty.  Lastly,  there  needs  to  be  a 
continual  supply  of  water  through  the  interior  of  the  body  to  the 
evaporating  surface.  This  transport  depends  upon  both  the  content  and 
the  potential  of  water  in  the  body  as  well  as  upon  its  conductive 
properties.  The  water  transport  takes  place  both  in  the  liquid  and 
the  vapor  phase,  involving  all  or  some  combination  of  pressure,  gravi- 
tational, temperature,  and  osmotic  potential  gradients.  Accordingly, 
the  actual  evaporation  rate  Is  determined  either  by  external  evapora- 
tlvity  or  by  the  soil's  ability  to  deliver  water,  which  ever  Is 
smaller. 

Three  Stages  Df  Drying 

In  laboratory  studies  where  evaporatlvlty  was  held  constant,  the 
soil  drying  process  nas  been  observed  to  occur  In  three  recognizable 
stages  (Pearse  et  al.,  1949).  These  are  the  constant  rate  stage,  the 
falling  rate  stage  and  the  slow  rate  stage. 

The  constant  rate  stage  occurs  early  in  the  process,  when  the 
soil  Is  wet  and  conductive  enough  to  supply  water  to  the  site  of 
evaporation  to  meet  the  evaporative  demand.  During  this  stage,  evapor- 
ation Is  considered  to  occur  at  potential  rate  and  is  limited  by,  and 
hence  controlled  by,  the  external  conditions  such  as  radiation,  wind 
and  humidity.  Duration  of  the  constant  rate  stage  depends  upon  the 
Intensity  of  mteorological  factors,  as  well  as  upon  the  soil 
conductive  properties.  Under  similar  external  conditions,  the  first 
stage  will  be  sustained  longer  In  a clayey  than  in  a sandy  soil,  since 


clayey  soil  retains  liiglier  wetness  and  conductivity  values  as  suction 
develops  1n  the  upper  tone  of  the  profile. 

In  the  falling  race  stage  the  evaporation  rate  falls  below  the 
potential  rate.  In  this  stage  the  evaporation  rate  is  limited  or 
dictated  by  the  rate  at  which  the  gradually  drying  soil  profile  can 
deliver  inolsture  toward  the  evaporating  tone.  This  stage  may  persist 
for  a much  longer  period  than  the  first  stage. 

The  second  stage  blends  Into  the  third  stage  so  gradually  that 
these  two  cannot  be  separated  easily  (Srover,  1972).  The  third  slow 
rate  stage  may  exist  at  a nearly  steady  race  for  many  days.  This 
apparently  occurs  when  the  surface  zone  has  become  so  desiccated  that 
transfer  of  liguld  water  through  it  effectively  ceases.  Water  trans- 
nition  through  the  desiccated  soil  thereafter  occurs  primarily  by  the 
slow  process  of  vapor  diffusion,  and  it  is  affected  by  the  vapor  diff- 
usivity  of  Che  dried  surface  zone  and  by  Che  adsorptive  forces  acting 
over  molecular  distances  at  the  particle  surfaces.  This  stage  can  be 
important  where  the  surface  layer  becomes  quickly  desiccated. 

A qualitative  explanation  of  the  occurrence  of  these  stages  was 
given  by  hillel  (1980)  using  Oarcy*s  law,  Ouring  the  first  stage,  the 
soil  surface  gradually  dries  out  and  soil  water  is  drawn  upward  In 
response  to  steepening  evaporation  induced  soil  water  potential  grad- 
ients. The  rate  of  evaporation  can  r^aln  nearly  constant  as  long  as 
the  soil  water  gradients  toward  Che  surface  ccmipensaCe  for  the 
decreasing  soil  hydraulic  conductivity.  Howeven,  when  the  soil 
sunface  appnoaches  equilibrium  with  Che  atmosphere,  the  soil  water 
content  gradient  coward  the  surface  cannot  increase  any  more,  and  in 
fact,  tends  to  decrease  as  the  soil  in  depth  loses  more  and  mre 


«3ter.  As  tne  evaporation  process  conclniies,  boCP  ctie  gradients  and 
the  conductivities  at  each  depth  near  the  surface  decrease  simultane- 
ously and  as  a result,  the  flux  toward  the  surface  and  the  evaporation 
rate  decrease  as  well.  Continuation  of  the  evaporation  process  for  a 
prolonged  period  1s  sometimes  accompanied  Oy  the  downward  movement  of 
a "drying  front*  into  the  profile  and  the  development  of  a distinct 
desiccated  zone,  through  which  water  can  move  from  the  still  moist 
underlying  layers  only  by  vapor  diffusion. 

Stage  specific  methods  have  evolved  for  computing  evaporation. 
They  range  from  simple  empirical  methods  to  more  complex  numerical 
models  requiring  large  consiuters.  In  the  next  section  these  methods 
are  discussed. 

Constant  rate  stage.  Evaporation  during  the  constant  rate  stage 
is  somewhat  analogous  to  evaporation  from  a free  water  surface  but  may 
usually  be  less  (van  Bavel  and  Hillel.  1976).  The  rate  of  evaporation 
during  this  stage  can  be  calculated  using:  (1)  an  empirical  function 
of  time  (Penman,  igai);  (2)  a constant  fraction  of  pan  evaporation 
(Ritchie,  1972;  Kanemasu  et  al.,  1976),  or  (3)  using  physical  models 
such  as  the  combination  method  based  on  measurements  in  the  atmosphere 
(Penman.  19*8;  Konteith,  1965;  Conaway  and  van  Bavel,  1967).  In  pra- 
ctical applications  of  these  procedures,  the  duration  of  this  stage 
duration  must  be  known  ahead  of  time.  The  duration  may  vary  from  less 
than  a day  to  a few  days  after  rainfall  depending  on  the  micrometeoro- 
logical  conditions,  Ritchie  (1972)  found  that  cumulative  eveporation 
following  a rainfall  event,  rather  than  time,  better  described  durat- 
ion of  the  first  stage.  According  to  references  cited  by  Ritchie,  the 
first  stage  ended  when  cumulative  evaporation  was  12  nm  In  case  of 


Melanto  clay  loam,  9 mn  for  Yolo  loam,  and  6 nm  1n  casa  of  Houston 
blacK  clay  and  ?la1nfleld  sand  respectively.  Ritchie  {1972),  Tanner 
and  Jury  (1975),  Kanemasu  et  a)  {1976)  and  Al-Khafaf  et  al.  (1978) 
have  used  this  procedure. 

Models  of  the  falling  rate  ttaee  based  on  analytical  solution. 
To  seme  extent.  It  is  believed  that  evaporation  during  the  second 
stage  is  not  affected  by  environmental  changes,  and  that  evaporation 
depends  on  the  soil  »ater  content  distribution  (Ritchie,  1972;  Hlllel, 
1980).  Ritchie  and  Burnett  (1971)  found  evaporation  to  be  closely 
related  to  the  eater  content  of  the  top  0 to  30  me  depth  in  Houston 
Black  clay.  Evaporation  decreased  rapidly  until  water  content  reached 
20t,  and  thereafter  the  dally  evaporation  rates  ranged  from  0,50  to 
0.25  nin/day  and  were  practically  Independent  of  above  ground  condit- 
ions. Therefore  a favored  approach  has  been  to  describe  the  cumulative 
evaporation  (CE)  during  this  stage  as  a function  of  tine  only  (Penman, 
1941;  Wiegand  and  Taylor,  1961;  Philips  ,1957;  Gardner, 1959),  Philips 
(1957)  and  Gardner  (1959)  analytically  solved  the  eguatlon  of  unsatu- 
rated flow  with  diffusivlty  as  an  exponential  function  of  water 
content  and  showed  that  CE  was  proportional  to  the  sguare  root  of 
tine.  It  Is  believed  that  this  approach  nay  apply  under  few  relati- 
vely Invariant  conditions  such  as  under  fully  developed  canopy  or  In 
laboratory  situations. 

Due  to  the  non-periodic  boundary  condition  used  In  deriving  this 
model,  its  use  was  limited  to  dally  evaporation  predictions,  Covey 
and  Bloodworth  (1966)  have  questioned  the  assumptions  Chat  external 
evaporativlty  was  unimportant.  The  assumption  of  fixed  boundary  cond- 
itions at  the  onset  of  an  evaporation  cycle  may  also  be  questioned. 
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Theoretically  the  solution  of  periodic  boundary  condition  under  iso* 
thermal  flow  condition  can  be  separated  into  a steady  state  tern  and  a 
periodic  tern.  However,  when  considering  evaporation  over  tines 
longer  than  2<  hours,  Che  steady  state  solution  is  valid  if  the  appro- 
priate boundary  condition  was  used.  This  was  the  rationale  of  Gardner 
(19S9).  Blact  et  a1.  (1969)  predicted  cumulative  evaporation  within 
of  true  values  for  the  Plainfield  sand;  however,  they  noted  that 
such  a good  agreement  nay  have  been  circumstantial  and  cautioned  that 
evaporation  from  finer  textured  soils  cannot  in  general  be  described 
Oy  such  a simple  nodel.  In  their  study,  the  evaporation  rate  was  so 
low  that  for  depths  of  wetting  and  time  involved,  the  soil  profile 
behaved  as  though  it  were  wet  infinitely  deep.  They  observed  that  C£ 
departed  from  the  linear  relationship  10  days  after  Irrigation, 
Gardner  and  Gardner  (1969)  while  working  with  finer  textured  soil  also 
made  similar  observations.  Kijne  (1975)  found  that  CE  during  the 
first  four  days  was  best  described  by  using  a coefficient  of  0.67, 
whereas  subsequent  period  CE  was  best  described  by  using  0,4.  It  was 
not  surprising  to  see  these  differences,  since  the  time  dependent 
model  cannot  account  for  variations  in  evaporativity,  irrigation  appl- 
ied, or  status  of  soil  water.  In  spite  of  the  problems  of  the  time 
dependent  mode!,  it  has  the  appeal  of  simpllcicy  and  has  been  used  by 
mtchie  (1972),  Ksnemasu  et  al.  (1976),  and  Al-Khafaf  et  al.  (1978). 

The  solution  of  soil  evaporation  problem  following  the  first 
stage  results  in  cumulative  evaporation  as  a linear  function  of  square 
root  of  time.  However,  this  solution  and  its  experimental  verifica- 
tion refer  to  idealized  situations.  In  this  approach  the  Initial  soil 
water  content  was  assumed  to  be  uniform.  But  the  water  content 
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aistribuUon  at  the  beginning  of  the  second  stage  is  rarely  uniform, 
and  it  Is  dear  that  suhseguent  evaporation  rates  must  depend  upon 
It.  At  the  end  of  the  first  stage  tne  soil  surface  is  more  or  less  at 
its  final  constant  water  content.  However,  the  soil  water  content 
distribution  in  the  profile  would  depend  upon  previous  evaporation 
rates  and  on  the  duration  of  the  first  stage,  fiardner  and  Hillei 
(1962)  observed  in  a laboratory  study  that  the  effect  of  initial  soil 
water  content  distribution  tends  to  become  small  after  some  time 
beyond  which  the  evaporation  rate  depends  only  on  the  total  remaining 
water.  Therefore,  the  use  of  a strictly  time  dependent  model  may  not 
give  true  evaporation  predictions  in  naturally  varying  environmental 
conditions. 

Numerical  models  of  tne  falling  rate  stage.  Perhaps  the  most 
satisfactory  answer  to  computing  evaporation  lies  within  combined  heat 
and  water  transfer  models.  Although  the  energy  and  mass  balances  are 
closely  interrelated,  most  previous  work  did  not  fully  take  this 
coupling  into  account,  [luring  the  last  two  decades  attention  has  Peer 
given  to  the  phenomenon  of  combined  heat  and  mass  transfer  in  soils 
(Philip  and  Oe  Vries,  1957;  Cary,  1965).  Solutions  of  the  flow  epua- 
tion  under  isothermal  flow  conditions  predict  many  of  the  essential 
features  of  the  first  and  second  stage  of  evaporation.  According  to 
Philip  (1967),  the  isothermal  model  can  be  expected  to  simulate  evapo- 
ration process  ( within  lOt)  so  long  as  Che  surface  is  not  extremely 
dry.  FritCdn  et  al.  (1970)  concluded  that  the  isothermal  flow  model 
cannot  describe  the  soil  water  distribution  in  the  case  of  nonisother- 
mal  fIcMS,  such  as  near  the  surface,  and  simultaneous  neat  and  mass 
transfer  analysis  would  give  better  approximations  to  the  experimental 


ddta.  Thus,  the  conDined  processes  of  heat  and  mass  transfer  have 
speclaT  relevance  to  evapotransplratlen  of  the  arid  regions,  and  temp- 
erature regulation  of  seed  beds  (Bloom  et  a1.,  1972). 

Hillel  (1975)  developed  a dynamic  simulation  model  based  on  a 
nimericai  solution  of  the  water  flow  equation  under  isothermal  condi- 
tions to  assess  the  effects  of  constant  and  varying  evaporativity  on 
the  cumulative  evaporation  and  soil  water  content  distribution  in 
space  and  time.  Hillel  used  evaporativity  as  an  external  forcing 
function  and  evaporation  was  computed  as  the  minimum  of  externally 
induced  evaporation  and  transmission  of  soil  water  permitted  by  the 
soil.  The  model  ignored  surface  radiation  balance,  conduction  of  heat 
into  the  soil,  exchange  of  heat  between  the  soil  and  the  atmosphere, 
and  tnermdlly  induced  transport  of  vapor  in  the  surface  zone.  In  a 
laten  development  Van  Bavel  and  KiHel  (1975)  considered  the  surface 
energy  balance  of  the  top  10  m of  soil  and  combined  it  with  the  tra- 
nsport of  heat  and  water  in  the  atmospheric  boundary  layer  and  in  the 
soil.  The  results  of  the  model  suggested  that  even  during  the  early 
stages  of  evaporation,  soil  acts  as  an  active  participant  In  determin- 
ing the  surface  energy  balance  and  hence  daily  evaporation  rates. 

Numerical  models  of  simultaneous  heat  and  mass  balance  based  on 
Pick's  Law  or  some  modification  of  it  have  been  developed  by  a number 
of  researchers  (Philip  and  De  Vries.  1957;  Cary,  1965;  Grover,  1972; 
Van  Keulen,  1975;  Van  Loon  and  Mosten,  1979).  Cary  and  Taylor  (1962) 
and  Cary  (1965)  Introduced  the  concept  of  irreversible  thermodynamics 
for  studying  simultaneous  transfer  of  heat  and  water.  However,  most 
of  these  studies  were  simulation  studies,  with  limited  or  no  testing 
under  field  situations.  These  procedures  have  greatly  helped  advance 
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tne  multlpnaee  soil  water  transgort  phenomena, 
however  they  are  complei,  and  their  utility  to  solve  field  problems 
remains  mostly  unchecked.  Additionally,  simultaneous  analysis  of  soil 
water  and  heat  flow  In  soils  requires  detailed  characterization  of 
water  content  and  thermal  properties  of  soil  materials.  Obtaining 
this  Information  Is  laborious  and  the  accuracy  of  the  values  Is  often 
questionable.  Also,  this  type  of  model  operates  with  time  Intervals 
of  the  order  of  minutes  or  maybe  seconds,  making  their  application 
unsuitable  in  seasonal  models  where  simulation  Is  performed  for  a 
large  number  of  days. 

Atmospheric  models  of  the  falling  rate  stage.  In  the  foregoing 
discussion,  emphasis  was  on  the  methods  evolved  from  solutions  of  the 
flow  equation  for  isothermal  and  non-lsothermal  conditions.  In  the 
atmospheric  models,  evaporation  from  a surface  1s  computed  using  the 
vapor  pressure  deficit  between  the  surface  and  Its  sunroundlngs  and 
the  surface  resistance.  The  driving  force  for  evaporation  Is  the 
vapor  pressure  gradient  between  the  a1r  Into  which  vapor  Is  moving  and 
Che  soil  from  which  evaporation  Is  taking  place.  If  e(T^]  Is  tne 
vapor  pressure  of  the  air  at  dew  point  temperature  T^,  and  e(Tj)  Is 
Che  vapor  pressure  at  the  soil  surface  with  temperature  T , then  eva- 
poration rate  E can  be  defined  as. 


where  x Is  the  latent  heat  of  vaporization  of  water,  Cp  is  the  volu- 
metric heat  capacity  of  air  at  constant  pressure,  t Is  the  psychro- 
metric  constant,  RS  Is  the  resistance  of  the 


laminar  boundary  layer 
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adjacent  to  soil,  and  US  Is  fie  soil  resistance  to  vapor  transfer.  RB 
depends  open  wind  velocity  and  roughness  characteristics  of  the  soil 
surface,  whereas  R6  1s  a function  of  soil  hydraulic  properties  and  Is 
characterlted  by  the  soil  water  content  status  and  depends  on  the 
history  of  evaporation.  Therefore  RS  Is  approximately  zero  as  long  as 
the  soil  surface  allows  free  evaporation  (E  controlled  by  environment 
only)  and  becomes  several  orders  of  nagnitude  higher  than  RB  when  the 
soil  Is  dry.  Atmospheric  models  have  heeh  used  by  Conaway  and  van 
Ravel  (1967),  Tanner  and  Fuchs  (196B),  and  Novak  and  Black  (1985)  to 
predict  evaporation  from  bare  well-irrigated  soil  surfaces,  and  by 
Soudriaan  (1977),  and  Shawcroft  et  a1.  (1974)  to  compute  evaporation 
from  soil  underneath  a canopy. 

Statement  Of  Bie  Problem 

Soil  evaporation  and  transpiration  models  are  combined  to 
estimate  amounts  of  water  lost  from  a crop  as  soil  evaporation  and 
transpiration.  Uhen  the  plants  are  In  the  early  growth  stage  with 
little  vegetative  cover,  the  evapotransplratlon,  ET,  from  the  entire 
crop  Is  dominated  by  soil  evaporation.  As  the  canopy  cover  Increases, 
transpiration  becomes  the  dominant  component  of  ET,  In  earlier 
atteepts  (Black  et  al.,  1969;  Ritchie,  1972)  evaporation  and  transp- 
iration were  computed  Independently  and  added  together  algebraically 
to  give  ET.  In  essence.  It  was  assianed  that  evaporation  was 
independent  of  the  microclimate  surrounding  the  plant.  It  has  been 
shown  that  the  presence  of  even  a thin  layer  of  mulch  to  suppress  E 
leads  to  significant  changes  In  the  crop's  radiation  balance.  The 
elimination  of  evaporation  will  reduce  the  humidity  of  the  air  within 
the  crop  canopy,  which  may  lead  to  a compensating  Increase  In  Che 


transplrstfon  rate.  These  ohservations  suggest  that  evapotranspiratloe 
(ET)  Thorn  a crop  is  controlled  by  both  the  soil  and  plant  factors  as 
•rell  as  the  nlcrocl inate  surrounding  ttiern.  Therefore  one  of  the  more 
promising  approaches  in  computing  ET  is  where  the  so1 1-p1ant-3t<Bos- 
phere  was  treated  simultaneously  as  an  energy  and  mass  exchange 

Shawcroft  et  a1.  (1971),  and  Soudriaan  (1977)  have  achieved  the 
coupling  between  E and  T in  their  micrometeorological  models  by 
solving  Egn.  (2.1)  for  the  canopy  and  the  soil  simultaneously.  In 
these  studies  it  was  assumed  that  when  soil  was  saturated  RS  was  zero, 
and  it  approaches  a large  value  when  soil  becomes  completely  dry.  Ilo 
attempt  has  been  made  to  determine  RS  for  Intermediate  soil  water 
conditions.  Due  to  this  difficulty  the  scope  of  application  of  these 
models  may  be  limited  to  eithen  fully  saturated  or  dry  soil 
conditions.  In  a field  situation,  soils  go  through  repeated  cycles  of 
wetting  and  drying  following  rain  andfor  irrigation.  Therefore,  a 
physical  model  of  the  soil  evaporation  process  which  accounts  for 
fluctuating  soil  water  content  conditions  was  needed  so  that  it  can  be 
llnxed  with  soil-canopy  models  for  predicting  ET. 

Eqn.  (2.1)  is  routinely  used  to  compute  evaporation  from  a single 
leaf  or  from  a canopy  with  RS  replaced  by  canopy  resistance  (Honteitn, 
19SS).  For  a single  leaf.  RS  is  estimated  using  a steady  state  para- 
meter, where  It  is  internally  computed  from  a measured  evaporation 
rate,  computed  e(T,),  e(T^),  and  R8.  for  a leaf  this  Is  a valid  appr- 
oach since  the  leaf  surface  temperature  is  easily  measured,  and  as 
long  as  a leaflet  is  thin.  It  can  be  assumed  to  be  at  the  surface 
measurement  techniques  for  RS  have  not  been 


temperature.  Such 


developed  because  tbe  assumption  of  unlfom  body  temperature  cannot  be 
extended  to  soil  due  to  tne  finite  thickness  of  drying  and  large  heat 
storage  capacity.  Therefore,  an  alternate  computational  procedure  was 
proposed.  Epn.  (2.1)  can  be  rearranged  to  solve  for  RV; 


where  RV  ■ RS  * RS  1s  the  total  soil  resistance  to  vapor  transport. 
All  variables  of  Eqn.  2.2  can  be  measured  except  Tj.  when  evaporation 
occurs  at  depths  below  the  surface.  Therefore,  to  facilitate  a pract- 
ical solution  of  Eqn.  2.2,  it  was  assumed  that  the  vapor  pressure  at 
the  true  evaporating  boundary  may  be  approximated  by  the  saturated 
vapor  pressure  at  the  soil  surface  temperature.  With  this  approximat- 
ion, Eqn.  (2.2)  can  be  solved  during  many  drying  cycles  to  give  a tine 
series  in  R».  This  time  series  would  then  be  developed  Into  a 
suitable  model  of  RV. 

following  the  first  stage  drying,  the  soil  surface  is  more  or 
less  at  tne  Its  final  constant  water  content,  and  the  vapor  pressure 
Is  some  constant  fraction  of  the  saturated  surface  vapor  pressure. 
This  fraction  depends  upon  soil  water  natric  potential  and  may  be 
coiiv)uted  1f  tne  water  content  is  known.  However,  as  long  as  surface 
water  content  does  not  vary  appreciably,  errors  introduced  by  tbe  use 
of  surface  temperature  would  be  consistent  throughout  and  can  be  acc- 
ounted by  an  empirical  model.  Hence,  the  assumption  is  not  likely  to 
cause  serious  errors  In  the  soil  evaporation  predictions. 
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Objectives 

The  objectives  of  this  stiioy  were  to: 

1.  Design  and  conduct  experiments  to  measure  E,  T,,  and  Tj,  and 
compute  RV  for  a drying  sandy  soil  subjected  Co  different 
Irrigation  management  practices  and  exposed  to  several 
levels  of  evaporatlvHy. 

2.  Obtain  an  empirical  model  of  the  total  soil  resistance 
(RV  • RB  ♦ RS),  applicable  over  a Hide  range  of  soil  nacer 
content  and  environmental  conditions. 

3.  Verify  tne  soil  evaporation  model  (Eqn.  2.1)  Hith  RV  given 
by  the  model  developed  In  (2),  by  conparing  the  predicted 
and  measured  evaporation  rates  as  well  as  cumulative  evapo- 
ration over  the  drying  cycles. 

«.  Perform  error  and  sensitivity  analysis  of  the  model  to  nigh- 
light  strengths  and  neatnesses  and  Identify  further  research 

Research  Method 

Evaporativity  determines  the  maximum  flux  at  which  the  atmosphere 
can  vaporiie  water  from  a free  water  surface.  The  magnitude  of  evapo- 
ratlvity  is  determined  by  the  air  temperature,  relative  humidity,  wind 
velocity,  and  radiation.  On  the  other  hand,  evaporation  from  a soli 
1s  determined  by  evaporativity  and/or  by  soil’s  ability  tb  supply 
water  to  the  evaporation  site.  Jackson  et  a1.  (1973)  found  chat  water 
content  of  0 - 0.09  m of  soil  layer  fluctuated  1n  response  to  the 
evaporativity  up  to  37  days  after  irrigation.  Hlllei  (1980)  stated 
that  the  length  of  time  and  the  duration  of  the  constant  rate  stage 
depends  upon  evaporative  intensity  as  well  as  soli  conductive 
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properties.  Results  of  a laboratory  study  by  Gardner  (1963)  suggest 
ttiat  higher  initial  soil  »ater  content  would  result  In  higher 
cumulative  soil  evaporation  over  the  same  time  period.  Therefore,  It 
was  concluded  that  variations  in  day-to-day  evaporativlty.  rainfall 
amounts  and  frequencies,  soil  water  concent,  and  soil  hydraulic  prope- 
rtles  were  the  key  variables  for  the  soil  evaporation  study  planned. 
Equipment  and  instivnentatlon 

To  simulate  a diversity  In  evaporative  conditions  In  a relatively 
short  time  period,  bare  soil  evaporation  experiments  were  conducted  In 
outdoor  controlled  environment  chambers  so  that  factors  establishing 
levels  of  evaporativlty  (I.e.  temperature,  humidity,  and  wind  speed) 
can  be  controlled  and  monitored  accurately  at  desired  levels.  This 
facility  was  located  at  the  Irrigation  Research  and  Education  Park  on 
the  University  of  Florida  campus  In  Gainesville,  Florida  (29.4  °N, 
82,2  °U,  altitude  66.5  HSL).  Each  chamber  consisted  of  an  aerial 
acrylic  top  measuring  1x2x1. 6 m supported  by  frames  and  an  underground 
lysimeter  1 m deep.  The  lysimeters  were  filled  with  a deep,  well- 
drained  Arredondo  fine  sand,  a member  of  the  Hyperthermic,  coated 
family  of  Typic  Ijuartilpsaninents,  brought  from  an  adjacent  field. 
Five,  L ffl  long  alianinum  tubes  (0.06  n diameter]  were  horizontally 
Installed  In  the  lysimeter  at  0.15  m Intervals  from  the  top  of  the 
lysimeter  for  water  contents  measurement  with  a neutron  probe.  The 
lysimeters  were  also  equipped  with  a vaccum  type  drainage  facility  at 
the  bottom.  Each  chamber  was  Instrumented  with  a net  radiometer 

(mounted  0.6  m above  soil  surface  during  this  experiment), 
thermocouples  to  monitor  air  entry  and  exit  tanperacures, 

thermocouples  to  measure  soil 


temperatures 


0.01,  0.02. 


0.10 


tipping 


raingpuge 


collect  and  meaiure  water  condensed  from  the  outgoing  air.  All 
sensors  nere  monitored  with  a dedicated  mini  computer  operating  in 
real  tine.  By  suitaoly  selecting  air  dry  bulo  and  dew  point  tempera- 
tures. desired  relative  humidity  levels  were  maintained.  Air  in  the 
chamber  was  continuously  circulated  by  a fined  capacity  fans  that 
maintained  uniform  wind  velocity  over  the  soil  surface.  The  outputs 
from  various  sensors  was  stored  at  five  minute  Intervals  on  a magnetic 
drive,  printed  on  paper  and  the  summary  was  displayed  on  Che  screen 
for  visual  monitoring  of  the  system.  The  details  of  the  control 
system,  its  design  and  operations,  and  Information  on  various 
auxiliary  instrumentation  was  published  by  Jones  et  al.  (19B4b). 

Design  of  experiments 

Experiments  were  designed  to  provide  information  on  the  effects 
of  soil  water,  frequency  and  amounts  of  irrigation,  and  evaporativlty 
on  soli  evaporation  and  development  of  the  soil  resistance.  The  exp- 
erimental design  (Tables  2.1  - 2.3)  consisted  of  thnee  humidity  level 
experiments.  EXPl,  EXP2,  and  EXP3  and  tnree  Irrigation  treatments  In 
each  experiment.  Each  humidity  Irrigation  combination  was  replicated 
three  times  In  succession  (REPl,  REP2,  and  REPS).  An  earlier  prelimi- 
nary study  had  Indicated  that  when  the  soil  profile  Initially  at  Meld 
capacity  was  allowed  to  dry,  evaporation  approached  almost  zero  by  the 
seventh  day,  therefore,  duration  of  each  replication  was  one  week. 
During  the  entire  study  each  relative  humidity  level  experiment  was 
conducted  simultaneously  In  three  chambers  so  that  the  effects  of  all 
three  Irrlgatlor  treatments  could  be  evaluated  together.  The  three 
Irrigation  treatments  were:  a)  Irrigation  every  morning  to  bring  the 
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•.  2.l.NicrcMi«teoro1a9lcal  conditions  maintained  between  300  - 
2000  nouns  Ouninj  various  soil  evaporation  experiments  and 
replications  In  the  environment  controlled  chamber  under 
daily  irrigation  treatment.  Tabulated  values  include  date, 
-f  irrigation  (Tirr).  amount  of  Irrigation  (Irr,  mnl, 
"t  radiation  (RN,  m).  dally  cumulative  evaporation 
: air  temperature  (T.)  and  dew  point  temperature 

°C,  and  relative  humidity  (Rh,  t)  with  associated 
ndard  deviations,  [nmm  milH  meters). 


dally  n 

llil! 


Table  Z. 
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HlcrcMieteorologlcal  conditions  Tsintained  between  SQO  • 
2000  hours  during  various  soil  evaporation  experiments 
and  replications  in  the  environment  controlled  chamber 
under  irrigation  treatment  3.  Tabulated  values  include 
dace,  time  of  irrigation  (TIrr),  amount  of  Irrigation 
(Irr,  imi] , dally  net  radiation  (Rll,  imi),  daily 
cumulative  evaporation  (£^,  im);  air  temperature  (T,) 
and  dew  point  temperature  (Tg)  In  *0,  and  relative 


EXP  DATE  TIrr  Irr  RN 

REP  hr  on  nun  4 


Ig 


Table  ?.3  - continued. 


EXP  DATE  Tier  Irr  SN  E^  T,  sd  T.  sd  Rn  sd 
REP  hr  mi  rr.  ^ ’i 


entire  soil  profile  to  tne  fleTd  capacity  (11.5*  by  voliiiee);  b)  First 


day  irrigation  to  bring  entire  soil  profile  to  field  capacity  followed 
by  one  week  of  drying;  and  c)  Similar  to  (b)  Out  irrigation  amount  was 
doubled.  Following  the  first  three  ekperlments  another  one  week  etp- 
eriment  (EIP4)  was  conducted  to  study  evaporation  under  low  and  relat- 
ively constant  radiation  load  simulating  situations  of  a fully  develo- 
ped canopy.  To  simulate  this  situation  all  chambers  were  covered  with 

The  dally  mean  and  standard  deviation  of  the  day  time  (7  a.m,  to 
8 p.m.)  mlcrometeorologlcal  conditions  such  as  dry  bulb  temperature, 
the  dew  point  temperature,  and  the  relative  humidity,  as  well  as  the 
dally  measured  cumulative  evaporation  and  net  radiation  amounts  are 
listed  in  Tables  ?.l  - 2,3  for  each  of  the  three  Irrigation  treatments 
respectively.  As  a part  of  on-going  previous  projects,  the  software 
controlling  the  chambers  was  designed  to  flush  each  chamber  at  nights 
with  ambient  air  and  did  not  control  dew  point  or  dry  bulb 
temperatures.  Therefore  the  night  time  mlcrometeorologlcal  Inputs 
were  Ignored,  but  the  evaporation  rates  were  recorded  at  all  times  to 
determine  cumulative  molstune  loss.  This  loss  of  data  during  night 
hours  did  not  cause  any  problems  since  vapor  pressure  deficits  and 
night  time  moisture  loss  was  very  low. 

Prior  to  each  irrigating,  the  water  content  of  the  profile,  by 
depth,  was  determined  using  e neutron  probe  , and  the  calibration 
curve  reported  by  Jones  et  al.  (1985b).  Due  to  the  small  floor  area 
of  the  chambers  (2  m2)  jj  feasible  to  take  gravimetric  soil 
samples.  The  difference  between  measured  water  content  and  field 
capacity  water  content  determined  the  amount  of  irrigation  needed  to 


field  capacity.  Once 
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repUnlst)  tne  whsle  profile  at  least  to  the 
amount  was  determined,  sprinkier  cans  were  used  to  apply  the  irrigat- 
ion. Uditionai  water  content  measurements  were  made  every  other 
day.  The  soil  temperatures  below  the  surface  were  recorded  every  hour 
using  a Campbell  Scientific  data  logger  C9-S.  In  addition,  daily 
evaporation  was  measured  by  weighing  the  condensed  vapor  that  was 
collected  from  each  chamber. 

At  the  end  of  the  evperimental  period  the  data  were  transferred 
to  the  Agricultural  Engineering's  PRIhE  550-11  computer.  The  data 
were  then  averaged  (evaporation  and  net  radiation  sunned)  for  every  30 
minute  period  for  further  analysis. 

Overview  of  evoeriments 

The  experimental  period  was  characterized  by  a wide  array  of 
evaporativlty  conditions.  The  daily  irrigation  treatment,  TRl,  rec- 
eived largest  share  of  net  radiation,  while  the  magnitudes  of  RN  from 
other  two  treatments  were  lower  due  to  lower  surface  water  content  and 
higher  albedo.  As  the  initially  saturated  soil  surface  lost  water  by 
evaporation,  surface  albedo  Increased,  reflecting  increasingly  more 
and  more  shortwave  radiation.  Decrease  In  the  surface  soil  moisture 
also  caused  higher  soil  temperature  and  emitted  more  long  wave  radiat- 
ion thereby  reducing  net  radiation.  It  was  observed  that  the  ampli- 
tudes of  surface  temperatures  from  TR2  and  TR3  were  twice  that  of  TR  1 
on  the  last  day  of  evaporation;  for  Intermediate  periods  the  ampli- 
tudes were  steadily  increasing.  The  dally  evaporation  rates  for  TRl 
varied  in  accordance  with  the  deily  net  radiation,  while  for  TR2  and 
TR3  the  evaporation  rates  decreased  exponentially  initially,  with  more 
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gradual  rate  of  fall  later.  However,  periodic  Oehavior  was  visible 
even  on  the  seventh  day,  although  it  was  very  small  in  amplitude.  The 
data  also  showed  the  effect  of  large  heat  storage  capacity  of  the  soil 
on  the  phase  difference  between  the  evaporation  flux  and  surface 
ten^eratures.  Hhen  the  surface  was  saturated  the  phase  lag  was 
almost  zero,  aod  became  progressively  larger  as  soil  drying  continued. 
Hally  and  cijiulailve  evaporation 

The  daily  evaporation  rates  and  cumulative  evaporation  losses  for 
all  three  irrigation  treatments  were  plotted  as  a function  of  time. 
Some  of  the  observations  made  here  showed  disagreement  with  the 
conihpnly  held  notion  regarding  stages  of  soil  drying  and  the  effect  of 
atmospheric  conditions  on  the  course  of  soli  evaporation.  It  was 
found  that  initiation  of  a drying  stage,  transition  from  one  stage 
Into  neit,  and  stage  duration  was  a function  of  evaporativlty, 
available  soil  water,  and  the  history  of  evaporation,  and  not  simply  a 
time  function  as  is  often  assumed. 

The  nesults  of  EXPl  are  plotted  in  Fig.  a.l.  During  EXPUREPl, 
data  for  only  first  three  days  was  usable  due  to  problems  with  the 
control  and  measurement  system.  The  first  day  evaporation  from  each 
treatment  varied  in  accordance  with  the  irrigation  amount  applied  ; 
TR3  had  received  highest  irrigatloh  amount  and  evaporated  IDX  more 
than  TRE  and  361  more  than  TRl  during  EXRI-REP3.  This  probably 
occurred  due  to  longer  duration  of  the  first  stage  for  TR3  than  1n  the 
other  two  treatments.  In  relation  with  TSE  and  TR3.  evaporation  rates 
for  TRI  following  day  one  represented  day-to*day  fluctuations  In  the 
evaporativlty,  where  as  the  evaporation  rates  from  TRE  and  TR3  varied 
In  relation  to  the  irrigation  amount.  E for  TR3  were  coosisteotly 
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F1g,  2,1  Daily  and  cumulative  evapcraCIcin  under  the  three 
Irrigation  treatments  during  EliPl 
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nijnef  tnan  Tl)2  and  they  decreased  linearly  for  the  first  three  days, 
while  Z for  TR2  dropped  exponentially  following  the  first  day. 
Evaporation  rates  remained  practically  unchanged  i days  after 
Irrigation,  but  their  magnitudes  were  different  during  REP2  and  REPS. 

The  linearity  of  cumulative  evaporation,  CE,  for  TRl  depended  on 
day  to  day  variability  in  E.  uniformly  high  evaporativlty  during  REP3 
resulted  in  highest  CE  among  all  replication  of  EXPl.  The  CE  for  TR2 
lagged  CE  from  TRl  following  the  first  day,  whereas  CE  for  TR3  were 
about  same  up  to  two  days.  This  criteria  was  used  by  Ritchie  (1972) 
to  determine  transition  time  from  stage  1 of  drying  into  stage  2. 
although  the  TRl  was  not  a truly  well  watered  treatment,  the  results 
showed  that  the  transition  point  was  delayed  with  more  irrigation. 
For  Houston  Black  clay,  Ritchie  (1972)  found  that  stage  1 ended  when 
CE  was  4 rm  following  an  irrigation  or  rainfall  event  if  soil  was  at 
field  capacity  initially.  With  this  criteria,  the  sandy  soil  in  this 
experiment  required  two  limits,  one  for  each  irrigation  treatment. 

During  EXP2.  Fig.  2.2,  relative  htinidity  was  lowered  to  Increase 
evaporative  potential  above  that  of  naturally  occurring  evap- 
orativity.  The  natural  evaporativlty  during  REPl  and  REP3  was  fairly 
uniform  and  the  CE  as  well  as  E behaved  similarly.  On  the  first  day  E 
from  TS3  was  slightly  higher  than  TR2;  on  the  second  day  the  diff- 
erence was  highest,  and  it  decreased  rapidly  over  the  next  two  days, 
and  finally  became  almost  constant  after  fourth  day.  During  EXPI,  the 
drying  was  more  gradual  and  changes  in  day  to  day  evaporation  were 
less  severe.  There  were  large  variations  in  the  day  to  day  evaporativ- 
lty during  REP2  that  resulted  In  some  interesting  observations.  The  E 
for  TR2  decreased  linearly  up  to  four  days  unlike  all  previous  cases. 
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Kliere  as  E for  TH3,  the  decrease  was  linear  for  six  days.  A leak  due 
to  rain  Increased  E for  both  the  treatments  on  the  last  day.  The  two 
cloudy  days  during  EXPZ  - REP2  following  the  first  day  of  Intense 
radiation,  probably  allowed  soil  rehydration  and  maintained  soil  near 
the  surface  relatively  moist.  The  fourth  day  increase  in  evaporation 
load  caused  a sudden  depression  in  E for  TRZ.  Illustrating  the  influ- 
ence of  above  ground  conditions.  Additionally  It  can  De  also  seen  that 
CE  for  both  the  treatments  had  not  levelled  off  even  on  the  last  day 
wnich  again  Illustrates  the  point  that  under  a right  combination  of 
conditions,  the  above  ground  conditions  could  affect  the  evaporation 
process  well  after  the  stage  I was  over. 

Results  of  EXP3  [Fig.  2.3)  and  EXP4  (Fig.  2.4a)  snowed  many  of 
the  behaviors  seen  earlier.  Very  stable  environmental  conditions 
during  EXP3-HEP1  produced  results  as  If  experiments  were  conducted  in 
laboratory.  Low  evaporatlvUy  load  during  E4R4-REP1,  resulted  in  slow 
drying  and  the  drop  In  E was  very  gradual.  Stage  1 duration  was  about 
two  days  longer  than  all  previous  experiments. 

The  cumulative  evaporation  data  for  TR2  and  TR3  from  all  experi- 
ments were  combined  together  (Fig.  2.4b)  to  see  If  It  could  be 
described  by  a time  function  as  suggested  by  Gardner  (1959).  There 
was  wide  scatter  in  the  dally  evaporations  mounts.  Tlie  seventh  day 
CE  for  TR2  ranged  from  10.5  fim  to  22.5  nm  whereas  for  TR3  the  range 
was  13.5  tin  to  26.5  tin  which  amounts  to  more  Chan  lOOt  variation  In 
both  the  cases.  If  one  looks  at  the  Individual  CE  curve  it  may  seem 
possible  to  represent  them  by  a time  function,  but  it  would  be  very 
difficult  to  specify  one  general  relationship  that  would  be  satis- 
factory across  all  experiments  and  treatments.  The  scatter  In  the 
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fig.  2.3  Daily  and  cumulative  evaporation  under  ttie  tnree 
Irrigation  treatments  during  EXP3 


Fig.  2.4a  Daily  ana  cunnilatlve  evaparailon  iinfler  the  three 
rrigatlon  treatmenti  during  EXP4  - SEPl 


P1g.  2.4Q  Cumulative  evaporatien  data  for  the  entire  study 
period  for  irrigation  treatments  2 and  3 


data  was  the  result  of  different  irrigation  amounts,  the  initial  soil 
water  content  at  the  tine  of  Inrigation,  and  variations  In  the  evapor- 
ativlty.  Describing  these  data  using  a time  function  would  clearly 
be  a gross  approximation.  In  summary,  following  conclusions  were 
drawn  from  the  experiments; 

1.  The  duration  of  stage  I was  related  to  both  the  cumulative 
evaporation  following  irrigation  and  the  initially  available 
water  for  evaporation. 

Z,  Evaporation  rates  as  well  as  cumulative  evaporation 
amounts  were  a strong  function  of  amount  of  water  applied 
and  cannot  be  simply  described  as  a tine  function. 

3.  The  rate  of  decrease  in  evaporation  rates  depended  upon 
irrigation  amount  and  the  history  of  evaporation. 

4,  Evaporation  rates  were  influenced  by  radiation  when  soil  was 
subjected  to  cycles  of  high  and  low  radiation  loads. 

9.  The  physical  model  (Eqn.  2.1)  may  eliminate  some  or  all  of 
these  problems  if  soil  resistance  SV  1s  appropriately 
calibrated  to  account  for  conditions  in  the  soil  as  well  as 
In  the  atmosphere. 

Analysis  of  soil  resistance 

The  time  series  of  RV  obtained  by  solving  Eqn  (2.2)  was  plotted 
as  a function  of  time  and  cumulative  evaporation  to  determine  the 
appropriate  model  for  describing  the  data.  from  the  time  series 
analysis  It  was  concluded  that  RV  could  not  be  described  by  time  since 
1t  did  not  account  for  the  Intensity  of  evaporatlvity  as  well  as  soil 
water  status.  Some  of  the  problems  of  the  time  dependent  model  were 


partfellj'  solved  Dy  using  cumulative  evaporation  as  Independent  vari- 
aDle  Instead  of  tine.  However  this  approacn  was  Inadequate,  because 
It  did  not  account  for  variations  1n  irrigation  amounts  [note  diff- 
erences between  TR2  and  TR3  ) as  welt  as  residual  soil  water  from 
previous  Irrigations. 

Sardner  (1970)  employed  a ratio  of  cumulative  evaporation  to 
amount  of  water  applied  as  tbe  Independent  variable  In  modeling  evapo- 
ration. This  approacn  would  work  If  there  1s  always  the  same  amount  of 
water  present  In  the  soil  prior  to  each  irrigation  and/or  that  all 
water  from  previous  application  had  evaporated  or  percolated  deep  Into 
the  soil  making  It  unavailable  for  subsequent  evaporation.  These 
requirenents  may  not  be  always  satisfied  1n  a field  situation  witn 
rainfall  occurring  at  random.  Additionally,  Gardner  and  Gardner 
(19S9)  nad  noted  that  at  low  application  rates  almost  all  applied 
water  as  well  as  some  additional  water  was  lost,  but  as  the  amount  per 
application  was  Increased  , the  fraction  evaporated  decreased  with 
water  moving  deeper  Into  the  soil,  suggesting  tnat  all  the  applied 
water  was  not  lost  during  subsequent  evaporation  periods. 

Therefore  even  for  a Initially  saturated  soil  there  must  be  an 
approximate  upper  limit,  CE^,  on  cumulative  evaporation.  CE„,  would 
depend  upon  physical  soil  properties  and  tne  maximum  depth.  D^,  over 
which  soil  evaporation  may  occur.  C£„  could  be  determined  for  a soil 
either  by  actual  measurements  or  by  simulation.  Some  procedures  have 
been  reported  for  determining  location  of  the  drying  front  but  In 
general,  measurement  of  evaporation  front  has  been  a difficult  task, 
Wiegand  and  Taylor  (1962)  monitored  temperature  profiles  of  a soil 
column  In  a laboratory  and  Inferred  that  the  location  of  the  evapo- 
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ration  frortt  was  associated  witii  the  largest  tenperature  depressiori; 
Bardner  and  Hanks  (1966)  used  a "heat  sink"  method  to  locate  the  eva- 
poration front.  Gnover  (1972)  experimentally  found  that  the  location 
of  the  evaporation  front  corresponded  to  a region  of  fixed  moisture 
content  (0.5t  over  0.01  m depth).  HcGowan  (1973)  Identified  the  eva- 
ponation  front  with  a plane  of  zero  moisture  flux  using  t1me-no1sture 
content  relationships:  HelUng  and  Bel)  (1980)  followed  a similar 
procedure  with  soil  moisture  data.  During  this  study  water  content 
profiles  were  measured  every  alternate  day  at  0.15  m Intervals  and 
temperatures  were  measured  up  to  0.10  m depth.  This  type  of 
resolutloh  seemed  inadequate  to  implement  any  of  the  aSove  procedures. 

Therefore  a numerical  model  was  developed.  Earlier,  similar 
models  have  been  published  by  Hanks  and  Gardner  (1965),  Hlllel  (1975), 
and  van  Keulen  (1976).  This  model  was  similar  to  HlHel's  model  but 
was  written  In  FORTRAN  77  and  used  experimentally  measured  evaporation 
rates  as  the  surface  flux  boundary  conditions.  The  model  soil  profile 
was  divided  into  eleven  layers  with  thickness  varying  from  0.02  m at 
the  top  to  0.175  m at  Ootton  of  profile.  Mater  flux  from  each  layer 
was  found  from  Darcy's  law.  The  results  of  water  content  with  depth 
for  each  seven  day  period  were  used  to  visually  identify  the  location 
of  the  evaporation  front  using  the  criteria  of  McGowan  (1973).  The 
results  of  twenty,  seven  day  periods  showed  that  the  plane  of  zeno 
flux  was  located  approximately  0.3D  m from  the  soil  surface.  It  was 
also  found  that  the  top  0.02  to  0.03  m soil  was  near  0 I moisture 
content  Out  lower  soil  was  relatively  moist  with  average  moisture 
content  of  56.  Using  this  basic  information  CE^  was  computed  (Table 


- 2.6)  for  all  Irrigation 
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Table  2.4  Procedure  for  conoutlng  CL,  nailnum  possible  cusKilatIve 
soil  evaporation  associated  with  an  Irrigation  or 
rainfall  event. 


The  Implementation  of  the  model  required  information  on  CE„,  the 
maKimum  possible  evaporation  associated  with  an  Irrigation  or  rainfall 
event.  The  CE,„  was  less  than  or  equal  CL.  the  potential  evaporation 
associated  with  a soil  Initially  irrigated  to  field  capacity  and  then 
subjected  to  drying.  When  soli  drying  occurs,  evaporation  usually 
stops  well  before  the  soil  Is  completely  dessicated.  If  the  average 
final  water  content  of  cm  of  soil  Is  denoted  as  9^  , then  CEp.soll 
can  be  defined  as 

CEp  . /Te,p  - e,)  da  f Ej  (2.3) 

where  z was  the  soil  depth,  field  capacity  volumetric  water  cont- 
ent . 9j  final  water  content  at  the  end  of  a drying  period,  and  Ej 
was  maximum  observed  first  day  evaporation  amount,  Results  of  numer- 
ical simulation  suggested  that  the  maximum  depth  of  drying,  p,  was  0.3 
m;  9f^  was  equal  to  11. 5X  (Zur  and  Jones.  1982)  , and  and  Ej 
were  measured  to  Oe  51  and  7.23  m respectively.  Substituting  these 
values  in  Eqn  2.3  gave  CEp  a value  of  26.73  im.  The  first  day  evapor- 
ation was  added  to  account  for  evaporation  loss  while  water  was  dis- 
tributed to  field  capacity.  Therefore,  E^-0  if  soil  Is  at  the  field 
capacity  at  the  beginning  of  the  drying.  This  computed  value  cotepared 
well  with  maximum  7 day  cumulative  evaporation  of  26.59  m during 
EXP2-REP2,  To  compute  CE^  associated  with  each  Irrigation  event,  the 
following  procedure  was  used: 


1.  Hessure  total  soil  water  prior  to  irrigation  (TSW|,^) 

Z.  Irrigate,  and  compute  total  soil  water  after  Irrigation, 
TSU,,>  TSW|ji  ♦ Irr 

3,  From  TSU,,  subtract  water  tliat  would  Oe  held  bjr  the  soil 
at  the  end  of  the  drying  period  (7  days  In  this  study), 

TSW  - TSWj,  • 

This  allows  for  drainage  below  the  evaporation  zone  as  well  as 
evaporation  water  loss. 

4.  Determine  CE,,  using  following  criteria: 

C£„  • TSW  If  TSW  <■  CEp 

CE„  ■ CEp  If  TSW  > CEp 

The  values  of  CE„  computed  using  above  four  criterias  are  listed  in 
Table  2.5  for  TR2  and  Table  2.6  for  TR3,  note  that  for  TR3,  CE^  values 
were  usually  at  the  mailmun  value  because  enough  water  was  added  to 
bring  the  evaporation  zone  to  field  capacity  and  allow  some  drainage 
below  that  depth. 


conputeil  for  T1I2 


ploCteP 


function  of 


nomalizM  cuTulatiye  evaporation  CE^,  (CE„  • CE/CE„)  1n  Figs.  2.5- 
2.8.  The  trend  in  progression  of  RV  during  different  experiments  was 
similar,  indicating  that  one  functional  form  may  be  adequate  to 
describe  R»  under  various  eiperimental  conditions  experienced  during 
this  study.  Initially  RV  remained  constant  for  about  7 to  8 hours 
following  an  irrigation  and  then  Increased  exponentially  as  C£. 
increased  beyond  a threshold,  CE^.  The  constant  phase  of  RV  corresp- 
onded to  the  stage  1 of  evaporation;  i.e.  when  RS  remained  zero.  The 
resistance  during  this  phase  was  only  frcm  the  laminar  air  layer 
adjacent  to  the  surface.  Since  the  experiments  in  this  study  were 
conducted  In  unlfonn  wind  velocity,  the  boundary  layer  resistance  RB 
remained  relatively  constant.  As  the  drying  continued,  the  soil  was 
no  longer  able  to  meet  the  atmospheric  demand,  and  to  maintain  the 
balance  between  the  demand  and  supply,  soil  resistance  RS  increased. 

RS  oscillated  showing  maxima  at  noon  and  minima  during  mornings 
and  evenings.  These  oscillations  may  have  partly  been  an  artifact  of 
computational  procedure  caused  by  Eqn.  2.2,  diurnal  oscillations  In 
surface  temperature,  evaporation  £.  and  the  phase  difference  In  the 
two.  The  moisture  Front  migrates  in  response  to  diurnal  variations  in 
the  atmospheric  demand,  the  front  being  farthest  from  the  surface  in 
the  afternoon  and  closer  to  the  surface  during  other  times.  An 
analogy  seems  apparent  between  the  oscillation  In  RS  and  location  of 
the  evaporation  front.  In  case  of  EXP2-REP2,  RV  remained  low  conpared 
to  all  other  cases.  This  situation  arose  because  of  alternating 
cycles  of  cloudy  and  clear  days  chat  caused  moisture  migration  and 
pernaps  rehydratlon  of  dry  upper  soil  layers.  These  dryihg  and 
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wetting  cycles  Taintelned  nigh  prolsnged  evaporation  rates,  Kept 
surrace  temperatures  low,  and  resulted  in  lower  RV  estlinates. 
Therefore  It  was  also  concluded  that  to  completely  describe  the  soil 
evaporation  process  some  provision  for  simulating  varying  sty 
conditions  was  needed.  Employing  a running  average  of  net  radiation 
was  found  to  be  a feasible  solution. 

In  sumiary  four  factors  were  needed  to  completely  describe  soil 
resistance  and  subsequently  the  evaporation  process  between  any  two 
water  application  events. 

1.  Cumulative  evaporation,  CE. 

2.  Amount  of  water  applied. 

3.  Available  soil  water  for  evaporation,  CE„. 

4.  Dally  running  average  of  net  radiation  following  an  Irriga- 
tion or  rainfall  event. 


Model  Develoonemt 

From  the  following  model  for  RV  was  proposed: 


If  CE^  < C£^  (2.4a) 
If  CE^  > CE^  (2.4b) 


where  a,  6 and  n were  constants  to  be  found  from  non-linear  least 
square  fitting  procedure  and  RN,  was  the  running  average  of  dally  net 
radiation.  Net  radiation  was  included  to  account  for  day  to  day  vari- 
ations In  the  sky  conditions  and  allow  simulation  of  renydratlon. 
Additionally  RB  and  CE^  were  also  unknowns  In  the  RV  equation.  In  the 
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literature,  methods  to  compute  RB  Based  on  wind  speed  and  surface 
roughness  are  recommended  {Honteith,  1965;  Goudriaan,  1977).  On  the 
other  hand  CE^  Is  a soil  specific  constant.  Since  the  soil 
Independent  Stage  of  evaporation  lasted  for  a day  in  most  cases,  the 
values  of  RV  during  the  first  day  were  analyzed  to  discern  value  of 
CEc  and  RB  using  following  procedure; 

1.  RV  was  plotted  as  a function  of  CE^  during  first  two  days 
following  irrigation, 

E.  RV  values  oscillated  during  the  day,  but  the  point  in  time 
when  RV  increased  and  continued  to  do  so  thereafter  was 
taken  as  indication  of  the  end  of  the  soil  independent  stage 
of  evaporation.  The  value  of  CE,^  at  this  time  corresponded 
to  CEj. 

3.  RB  was  computed  as  the  arithmetic  average  of  RV  values  up  to 

from  Che  analysis,  RB  had  a mean  of  1S3.B  s/m  with  standard  devi- 
ation of  33.1  s/m  and  80t  of  the  values  were  within  one  standard  devi- 
ation of  the  mean.  The  value  of  RB  compared  well  with  some  puDlished 
values  (Table  2.7).  From  this  analysis  CE^  was  15X  and  It  required 
7.5  hours  (sd  • 0.5  hours)  for  CE,.,  to  equal  this  value.  With  RB  and 
CE,.  determined,  optimal  values  of  o,  S and  n were  predicted  using  a 
non-linear  least  square  procedure; 


1.0942  * 10* 


-3 


(2.5a) 

(2.5b) 

(2.5c) 


8 ■ 2.0197  * 10  ■ ^ 


n ■ 1.36757 


CoTparison  of  boundary  layer  resistance  (s/cn)  obtained 
In  this  study  and  from  some  selected  literature  under 
similar  situations. 


This  study  Bare  soil 

wind  speed  0.085  m/s  153  8 

Soudriann,  1577  soil,  clod  site  0.05  m 138.0 

Eqn.  3.28  wind  speed’O.OBS  n/s 

Conaway  and  van  3ave1,  1969  0.0256  cm.  Z>1  m 

Edn.  3 w^nd  speedwlm/s  928.0 

most  wind  velocities 


Rijtenia.  1973 


choices  of  psraneters 


evaporation  laodel  became: 


•(?  ) - «(5  ) 
XT  1537 


CE  < 15  t (Z.ea) 


Id  <dien  CE  > 15  t E is  given  by; 


(2.6b) 


153.8  e 

The  model  predicted  values  of  RV  were  superimposed  on  Rigs,  2.5- 
2,8  for  direct  comparison  with  the  data.  The  model  seems  to  simulate 
RV  well  early  during  the  drying  cycle,  however  later  it  sisMjIates  some 
average  path  rather  than  the  oscillations  of  RV,  However  at  higher  RV 
values  . E would  be  very  small  and  was  not  expected  to  give  any 
noticeable  errons.  During  EXP  4 the  model  predicted  low  resistance 
perhaps  due  to  the  running  averaging  of  net  radiation  and  was  likely 
to  Introduce  some  errors. 


node!  Verification 

To  determine  the  accuracy  of  the  model  and  establish  confidence 
In  its  predictions.  It  was  used  to  compute  evaporation  in  response  to 
Irrigation  amounts,  soil  and  ain  temperatures,  and  evaporativity.  The 
hdurly  evaporation  rates  as  well  as  cumulative  evaporation  was  compar- 
ed witn  tne  experimental  measurements  in  Figs.  2.9-2.18. 

During  the  design  of  the  experiment  it  was  expected  that  the  soil 
independent  evaporation  stage  (stage  1)  could  be  maintained  for  at 
least  24  hours,  which  was  the  irrigation  frequency.  To  verify  this, 
simulations  were  carried  out  by  assuming  that  RV  • R8  between  two 
Irrigations.  The  model  simulated  evaporation  rates  Fairly  accurately 
up  to  6-7  hours  after  irrigation  and  then  resulted  in  overestimation 
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fig.  2.9  Comparison  of  measured  and  simulated 
evaporation,  EXPl  - REPl 
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2.13  Conparison  of  measured  and 
evaporation,  EXP2  - RCPl 
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2.14  Cooipdrlson  of  inaosi 
evaporation,  CXP2  - REP3 
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F1g.  2.16  Comparison  of  neasuretl  and  simulated 
evaporation.  EXP3  - REP2 
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F1g.  2.17  Codpdri son 
evaporation. 


measured  and 
REPS 


simulated 


rSe 

-40 


61 


2.18  Comparison  of  measi 
evaporation,  EXP4  . REPl 
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suggesting  chsc  the  stage  I duratign  was  smaller  than  expected 
initially.  It  was  found  that  stage  I lasted  an  average  of  7.5  hours. 

Perhaps  the  true  test  of  the  model  was  provided  by  TR2  and  TI13 
which  were  subjected  to  a week  of  drying.  The  model  accuracy  was  good 
early  in  the  drying  cycle  when  evaporation  rates  were  high.  The  model 
predictlohs  were  on  average  within  3*  of  the  measured  value  Oh  the 
first  day.  101  on  the  second  day.  and  151  on  the  third  day.  However 
the  model  was  less  sensitive  at  low  evaporation  rates.  An  equation 
relating  cunolatlve  evaporation  to  time,  Eqn.  2.7,  was  developed  to 
determine  if  the  soil  evaporation  model  developed  in  this  study  provi- 
ded improvement  over  currently  used  approache.  The  equation  was  fit 
tp  combined  data  of  TR2  and  TR3  using  hon-11near  least  square 
procedure.  The  data  for  TRI  were  not  Included  since  that  treatment 
had  received  daily  irrigation.  The  resultant  model  had  a R^  of  0.92 
and  the  exponent  of  time  was  0.325  which  was  different  than  0,5  as 
derived  by  Sardner  (1959)  and  Philips  (1957)  from  their  solutions  of 
the  flow  equation: 

CE.A.16t“-“=  (2.7) 

Kijne  (1975)  had  recommended  two  different  exponents.  0.67  for  the 
first  four  days  of  drying  and  O.A  for  15  to  20  days  there  after. 
Relative  errors  (Eqn.  2.8)  between  measured  dally  evaporation  rates 
and  predicted  evaporation  rates  using  Eqn,  2.6  and  2.7  were  compared 
In  Table  2.8  for  TR2  and  Table  2.9  for  TR3. 

Relative  Error  . ‘^Ai"hjlated  ‘^measured’  . (j.j, 

measured 

It  can  be  seen  from  these  tables  chat  the  Individual  as  well  as 
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medel  when  compared  with  the  time  dependent  model.  The  smaller 
standard  deviations  among  all  experiments  suggest  that  the  dally  eva* 
poratlon  rates  simulated  with  the  new  model  were  more  close  to  the 
measured  values  on  any  day  than  Che  time-model. 

Late  In  the  drying  cycle  the  relative  errors  and  standard  devia- 
tions were  high  with  the  new  model,  but  still  much  smaller  than  the 
time  dependent  model.  In  this  study  this  phenomena  was  observed  app- 
roximately three  days  after  drying,  when  approximately  SDt  of  evapora- 
tion had  been  completed  and  the  dally  evaporation  rates  were  much  less 
than  1 nm/day  and  therefore  the  errors  were  insignificant  In  terms  of 
total  water  lost.  The  RV  computed  from  the  data  showed  diurnal  osci- 
llations, whereas  the  model  predicts  an  average  RV  for  each  day. 
Since  the  evaporation  rates  were  small  during  these  periods,  this 
approximation  was  not  expected  to  give  serious  errors. 

The  maximum  errors  In  cumulative  evaporation  were  recorded  during 
EXP2-REP2  which  was  characterized  by  variable  cloudy  conditions.  The 
errors  were  heavily  biased  by  wrong  calculations  on  the  second  and  the 
third  day  in  both  the  cases.  Extreme  variations  in  the  above  ground 
conditions  may  have  produced  oscillations  in  the  evaporation  front, 
resulting  to  cycles  of  drying  and  rehydration  of  dried  zones  that 
resulted  In  sustained  evaporation  rates.  This  model  can  adjust  the 
rate  of  upward  trend  of  RV  under  low  radiation  conditions  but  does  not 
nave  any  mechanism  to  reduce  the  resistance  as  It  may  occur  due  to 
rehydration  of  soil.  To  Incorporate  such  a mechanism  in  the  model,  a 
combined  heat  and  moisture  balance  model  would  be  required. 

The  errors  (differences  between  measured  and  simulated  daily 
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evaporatlsn)  were  normally  distributed  suggesting  that  the  errors  were 
Indegendent.  The  data  base  developed  bad  covered  a wide  variety  of 
situations  and  performed  well  in  all  cases  except  under  extreme  situa- 
tions. The  periodicity  in  evaporation,  though  small  in  amplitude, 
prevailed  well  after  the  end  of  stage  I.  This  was  contrary  to  the 
classical  assumption  that  evaporation  following  stage  I was  environ- 
ment independent.  It  appeared  that  there  was  a transition  stage 
between  stage  I and  the  environment  Independent  stages  II  and  III,  and 
the  duration  of  this  stage  depended  not  only  on  the  atmospheric  demand 
but  also  on  the  strength  of  the  evaporation  source,  mainly  the  water 
content  of  the  evaporation  zone. 

Sensitivity  Unalvsis 

In  the  preceding  discussion  It  was  shown  that  the  new  model  was 
superior  compared  to  the  strictly  tine  dependent  model.  This  model 
allows  diurnal  simulation  of  evaporation  and  will  be  simple  to  inte- 
grate with  the  canopy  models  to  develop  an  Integrated  soil-plant 
-atmospheric  model  to  predict  E,  T and  ET  In  response  to  soil,  plant 
and  atmospheric  factors.  The  new  model  required  more  information  and 
the  additional  effort  in  obtaining  them  seemed  Justified.  The  success 
of  the  soil  evaporation  model  depends  upon  the  accuracy  with  which  the 
soli  resistance  model  parameters  RB.  CE^,  a,  a and  n are 
specified.  To  Illustrate  model  sensitivity  to  these  variables, 
seventh  day  CE  was  computed  by  varying  one  variable  at  a time.  For 
this  computation  three  one  week  periods  were  selected  based  on  the 
accuracy  with  idilch  they  were  simulated.  The  selections  were  the  best 
simulated  week  (EXP3-REPE,  TR3).  the  worst  simulated  week  (EXPE-REPE, 
TRE)  and  an  average  week  (EXFl-REPl,  TR3].  Computations  were  done  at 


the  model  assigned  values  for  these  parameters,  as  well  as  80t.  and 
130X  of  that  value.  The  relative  errors  In  seventh  day  CE  were 
compared  In  Table  2.10  The  negative  errors  In  the  table  indicate  tnat 
simulated  values  were  less  than  measured  values  and  positive  values 
indicate  over  estimation.  The  model  was  most  sensitive  to  n , and 
least  sensitive  to  s . The  Importance  of  RS  and  CE^  was  only  for  a 
short  period  In  drying  cycle,  and  since  the  predictions  were  satis- 
factory ( see  Figs.  2.9  - 2.18)  during  this  time,  the  errors  assoc- 
iated with  them  were  small.  The  effect  of  20  S variation  in  s 
resulted  In  only  n variation  In  CE,  which  also  suggest  that  Its  accu- 
racy was  not  that  critical.  The  most  errors  were  due  to  variations 
In  n , therefore  care  must  be  taken  in  determining  its  value.  In 
applying  this  model  to  other  soli  types  these  parameters  may  be  obta- 
ined either  by  direct  measurements  and  curve  fitting  as  done  in  this 
study,  or  a combined  heat  and  moisture  model  may  be  used  to  derive 
values  of  CE^,  a,  3 and  n using  model  fitting  prdcedores  similar  to 
the  ones  outlined  In  this  chapter. 

Direct  measurements  of  RB  has  been  difficult  and  the  disagreement 
between  various  measurements  Is  evident  from  Table  2.7,  where  R8  comp- 
uted in  this  study  was  compared  with  some  published  values.  Therefore 
perhaps  the  best  way  of  obtaining  it  Is  by  calibration  of  the  observed 
data.  The  model  Is  sensitive  to  RB  only  during  stage  I of  evaporation 
(RS*0),  following  this  stage,  RS  becomes  several  orders  of  magnitude 
higher  than  RB  and  RB  will  not  have  any  noticeable  impact  on  the  pred- 
ictions. In  sandy  soils  stage  I duration  is  small,  therefore  small 

errors  In  RB  will  not  have  much  Impact  on  the  long  tine  predictions. 
But  this  may  cause  appreciable  errors  in  heavy  soils  where  the  stage  I 
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Tabli  2,10  Results  of  sensitivity  analysis.  Percent  change  in  CE 
on  day  7 caused  by  t change  in  Che  best  oicked  par- 
ameters of  the  soil  evaporation  model  Egn.  2.6. 


CASE 


RS  CE,  a 8 

0.8  1.2  0.3  1.2  0.8  1.2  0.8  1.2 


0.8  1.2 


8£ST  5 -4  -4  4 0 0 7 -6  48  -26 

MORST  4 -4  -3  3 0 0 7 -5  55  -24 

AVERAGE  6 -4  -4  4 0 0 7 -G  49  -25 
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duration  is  likely  to  be  longer.  Therefore  the  procedure  used  to 
compute  RB  must  be  decided  accordingly. 

measures  the  duration  of  stage  I following  an  irrigation  or 
rainfall  event.  It  can  also  be  used  to  calibrate  a model  as  was  done 
in  this  study  so  that  the  simulated  and  measured  values  were 
similar.  The  calibration  procedure  predicted  CBc'lBI  in  this  study. 
However  this  choice  of  CE^  results  in  underprediction  of  evaporation 
from  EXP2-REP2,  probably  because  of  moisture  migration  that  sustained 
evaporation  rates. 

Summary  fnd  Conclusions 

The  need  to  predict  evaporation  and  transpiration  components  of  a 
developing  crop  canopy  was  identified  as  an  important  Issue  in  determ' 
Ining  appropriate  water  saving  strategies.  Various  technigues  differ- 
ing in  details  and  complealty  have  been  developed  td  tackle  this 
problem.  Evaporation  and  transpiration  can  be  computed  as  Independent 
components  and  added  together  to  give  ET  from  the  crop  (Ritchie,  1972) 
or  they  can  be  treated  simultaneously  es  interdependent  processes 
(Shawcroft  et  al.,  1970;  Goudriaan,  1977).  The  later  approach  seemed 
more  promising  because  It  treats  soil-pi  ant-atmosphere  as  an  energy 
exchange  system  and  has  potential  of  simplatlng  even  the  microclimate 
of  the  canopy.  To  Implement  this  later  approach,  a soil  evaporation 
model  based  on  vapor  pressure  difference  and  soil  resistance  was  des- 
ired. The  principles  applied  in  measuring  plant  resistance  cannot  be 
applied  to  a soil  due  to  the  lack  of  Information  of  the  precise  loca- 
tion of  the  drying  front  and  also  due  to  the  flhite  thickness  over 
which  evaporation  occurs.  Therefore  soli  resistance.  RV,  was 
ccvnputed  from  measured  evaporation  rates,  soil  temperature  and  a1r 
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tumidity.  By  as$iming  that  evaporation  was  propoptional  to  the  vapor 
pressure  difference  between  the  soil  surface  and  atmosphere.  A ten 
week  evaporation  study  was  completed  in  a cootrolled  environment  where 
son  was  subjected  to  natural  radiation  but  temperature  and  surround- 
ing humidity  were  controlled  to  give  three  evaporativlty  levels. 
During  the  study  the  soil  was  Irrigated  daily,  and  once  every  week. 
There  were  two  once  a week  Irrigated  treatments,  with  the  second  rece- 
iving twice  the  amount  of  the  first.  The  resulting  time  series  in  RV 
was  used  to  develop  a semi-empirical  model  of  RV,  the  soil 
resistance.  The  model  perrormed  well  and  the  errors  were  estimated  to 
be  much  smaller  than  a strictly  time  dependent  model  suggested  by 
fiardner  (1959). 

In  this  model  many  of  the  possible  factors  controlling  evapora- 
tion were  accounted  for.  This  model  has  the  characteristics  that, 
following  an  Irrigation  or  rainfall  event,  the  soil  will  start  out  In 
stage  I and  will  remain  in  stage  1 a variable  length  of  time  depending 
on  the  amount  of  water  added.  To  further  improve  the  accuracy  of  this 
model,  one  may  need  a comhined  heat  and  mass  transfer  model. 


CHAPTER 


HOBEL  FOR  PREDICTING  MATER  USE  AND  HICROCLIHATE 
OF  A DEVELOPING  CROP 

Introduction 

Crop  ;rle1d  is  set  Oy  the  genetic  potential  oF  the  crop,  ecologi- 
cal conditions,  and  by  proper  and  Cimely  supply  oF  water,  Fertilizer 
and  pest  control.  Water  Is  perhaps  the  most  important  oF  all  manage- 
ment variables  and  must  be  supplied  in  ouanclty  and  time  in  accordance 
with  the  crop  water  needs  For  optimum  growth  and  yields.  The  timing 
and  amount  of  water  applied  is  governed  by  prevailing  climatic  condi- 
tions, stage  of  crop  growtn,  soil  water  nolding  capacity  and  the 
extent  of  root  develotment.  Estimating  ET  for  a specific  crop  can  Pe 
a very  complex  issue  depending  upon  the  degree  of  accuracy  desired. 
To  obtain  most  accurate  estimates,  a11  of  the  above  crop  and  environ- 
mental conditions  need  to  be  taken  into  account. 

Since  Dalton  first  introduced  the  mass  transport  equation,  there 
have  been  numerous  models  and  methods  to  evaluate  crop  water  needs. 
Tnese  methods  Include  various  types  of  lysimeters,  portable  chambers, 
sap  Flow  techniques,  atmometers,  and  pan  evaporimeters  (Glad,  1082). 
Precision  weighing  lysimeters  provide  the  only  direct  way  of  accura- 
tely measuring  ET.  However  due  to  the  expense  and  time  involved  in 
using  field  measurement  methods,  predictive  means  of  estimating  ET 
with  varying  level  of  accuracy  and  complexity  have  evolved. 

Tne  choice  of  a ET  estimation  method  depends  upon  the  Intended 
use  and  tne  time  scale  of  application.  For  example  crop  irrigation 
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managenent  requires  only  daily  or  maybe  hourly  estimates  of  ET  to 
allow  timely  Irrigation  decisions  to  avoid  water  stress  and  yield 
reduction.  On  the  otherhand,  monthly  estimates  may  he  adequate  for  a 
regional  water  management  proOlera,  In  this  paper  short  term  methods 
are  considered. 

Thornthwaite  (19441  and  Penman  (1948)  introduced  the  concept  of 
potential  ET,  to  characterize  the  effect  of  climate  on  crop  water 
requirements.  Potential  ET  was  defined  as  the  rate  of  evapotranspira- 
tion  from  an  entensive  surface  8 to  15  cm  tall,  green  grass  cover  of 
uniform  height,  actively  growing,  completely  shading  the  ground  and 
not  short  of  water  (Jones  et  al.,  1984a).  Several  ET  methods  requir- 
ing meteorological  data  have  been  proposed  and  are  used  to  estimate 
potential  ( or  now  more  commonly  referred  to  as  the  reference  ) crop 
evapotnanpiration  (Ooorenbos  and  Pruitt,  1977).  The  simplest  among 
these  models  require  only  the  average  air  temperature,  day  length,  and 
crop  coefficients.  The  generally  better  performing  methods  need  daily 
radiation,  temperature,  vapor  pressure  and  wind  speed  data  (Jones  et 
al.,  1984).  The  actual  crop  water  requirement  is  computed  from  refer- 
ence Bvapotranspiration  using  a crop  coefficient,  which  is  characteri- 
stic of  the  crop,  and  its  stage  of  growth.  These  methods  have  been 
reviewed  along  with  their  potential  applications  by  Jensen  et  al, 
(197D),  Doorenbos  and  Pruitt  (1977),  Burman  et  al.  (1981),  and  by  81ad 
(1982). 

The  concept  of  potential  ET  is  useful  for  many  purposes  and  the 
crop  coefficient  method  has  gained  widespread  acceptance.  Potential 
ET  requires  that  plants  behave  passively  as  wicks  and  that  plants 
evercise  no  control  over  transpiration  (T).  In  reality  ET  Is  Influ- 
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enced  soil  and  plant  factors  and  consaquently  1t  is  loaar  than 
potential  ET.  Honteith  (196$)  combined  energy  balance,  aerodynamic 
parameters,  and  surface  parameters  in  an  eouation  valid  for  vegetation 
of  any  type  In  any  state  of  Mater  stress.  The  development  of  the 
Honteith  Tfiodel  Mas  a historic  event  and  the  model  has  become  the 
building  block  of  many  ET  models  for  closed  canopies  ( Zur  and  Jones, 
19B2;  Jagtap  et  al.,  1984;  Sinclair  et  al.,  1976)  and  has  been  widely 
quoted  in  the  literature.  However,  their  approach  is  vague  about  soil 
evaporation  and  as  a result  the  Honteith  model  has  not  been  applied  to 
the  problem  of  soil  evaporation  and  Incomplete  canopy  cover. 

Ground  cover  affects  ET.  Hhen  plants  cover  only  a small  portion 
of  the  soil  surface,  ET  is  dominated  by  E.  As  the  plant  cover  Increa- 
ses, evaporation  from  the  soli  surface  is  decreased  and  the  relative 
importance  of  transpiration  Increases.  Xlocke  et  al.  (1982)  in  a 
field  study  showed  that  soil  evaporation  was  reduced  with  canopy  deve- 
lopment but  it  did  not  cease  even  at  full  crop  cover.  Phene  at  al. 
(19SS)  reported  that  with  dry  soil  surface  and  constant  soil  water 
content,  the  bare  soli  evaporation  rate  was  14t  of  potential  ET  for 
alfalfa.  The  Importance  of  soli  evaporation  at  small  LAI  cannot  be 
Ignored  and  even  at  full  ground  cover,  soil  evaporation,  though  small. 
Influences  canopy  microclimate  (humidity  and  temperature). 

Recently  Black  et  al.  (1969),  Ritchie  (1972),  and  Tanner  and  Jury 
(1975)  developed  empirical  models  to  estimate  soil  evaporation  and 
transpiration  from  well  watered  crops  where  soil  goes  through  cycles 
of  drying  and  wetting.  In  these  models  transpiration  was  related  to 
leaf  area  development  and  was  computed  by  multiplying  potential  ET 
with  an  empirical  function  of  leaf  area  to  account  for  the  degree  of 


crop  cover.  The  soil  evaporation  was  computed  in  two  stages.  Follow- 
ing a rain  or  Irrigation  event  and  until  a certain  cumulative  evapora- 
tion occured,  E was  calculated  based  on  the  net  radiation  reaching  the 
soil  surface.  Frail  this  tine  on,  E was  determined  using  the  eguation 
developed  by  Black  et  a1.  (ig69)  relating  cumulative  evaporation  (CE) 
to  the  square  root  of  time.  The  transpiration  and  soil  evaporation 
coieponents  were  then  added  together  to  give  ET. 

The  separation  of  ET  into  E and  T and  their  mutual  interrelation- 
ship forms  an  important  problem  in  water  management.  Hodern  precision 
techniques  of  Irrigation  application  and  methods  of  evaporation  and 
transpiration  control  now  being  developed,  require  this  information  to 
allow  their  effect  on  evapotranspiratlon  to  be  evaluated.  One  of  the 
promising  approach  to  achieving  this  separation  Is  to  treat  the  soil- 
plant-atmosphere  continuum  as  an  energy  and  mass  exchange  system. 
Therefore  the  objectives  of  this  paper  were  to; 

1.  Develop  a mechanistic  model,  that  treats  soil  evaporation  and 
transpiration  as  separate  processes,  for  predicting  the  water 
use  and  microclimate  of  a well  watered  developing  crop,  where 
the  soil  surface  may  go  through  cycles  of  wetting  and  drying. 

2.  Develop  solution  for  the  mechanistic  model  that  will  allow 
its  use  with  normally  measured  climate  variables  so  that  1t 
can  be  used  for  practical  applications  and  evaluations. 

3.  Illustrate  with  an  example  the  various  prediction 
capabilities  of  the  model 

4.  Conduct  a sensitivity  analysis  of  the  model  parameters. 
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Woael  Oavelopmerit 

In  this  section  the  funnamental  equations  governing  the  energy 
exchange  aet*een  an  evaporating  body  and  its  surroundings  are 
presented.  The  discussion  will  explain  how  the  energy  exchanges  In 
the  canopy  and  soils  are  integrated,  and  assumptions  Involyed.  The 
proposed  model  includes  a number  of  rarely  measured  canopy  and  soli 
attributes,  so  algorithms  and  theory  required  to  obtain  them  from  some 
cotimonly  measured  inputs  will  Qe  presented. 

Fundamental  Equations 

The  evaporation  and  temperature  of  a single  leaf  can  be  predicted 
from  its  energy  balance  if  the  resistances  to  heat  and  vapor  transfer 
from  the  leaf  are  known.  Since  the  photosynthetic  heat  flux  and  time 
rate  of  canopy  heat  storage  are  usually  small  compared  to  other  energy 
balance  terms,  the  principal  of  energy  balance  states  that  the  net 
radiation  (RN)  absorbed  goes  intn  latent  (L)  heat  {evaporation  on 
dew),  and  the  baince  goes  into  sensible  (S)  heat  (neating  or  coaling 
of  air).  The  latent  neat  flux  density  from  a leaf  Is  proportional  to 
the  water  vapor  concentraion  gradient  between  the  inside  of  the  leaf 
and  turbulent  air  just  outside  the  boundary  layer  and  is  inversely 
proportional  to  the  simi  of  the  stomatal  and  boundary  layer 
resistance.  Similarly,  sensible  neat  flux  density  is  given  by  tempe- 
rature gradient  ecross  the  boundary  layer  resistance,  hathematical ly. 


(3.1) 


(3.3) 


WBere  c,  is  the  volunetric  heat  capacity  of  the  air,  y is  the  psychr- 
ometrlc  constant,  e(T^)  is  saturated  vapor  pressure  at  leaf  surface 
temperature  Tj.  e^  and  are  vapor  pressure  and  temperature  of  the 
air,  ano  r,  and  r|,  the  stomata!  and  boundary  layer  resistances  respec- 
tively. To  compute  the  evaporation  loss  from  the  leaf,  the  sensible 
heat  flux  density  S.  is  replaced  by  mathematical  substitution  using 
the  saturation  vapor  pressure  vs.  temperature  relationship, 
e(T,)  - e,  . i (T,  - T,)  * (e(T,  - e,) 

■ 4 (Tj  - T,)  ♦ 0 (3.4) 


eherw  D is  vapor  pressure  deficit  or  dryness  of  the  air.  Monteith 
(1965)  for  Che  benefit  of  simplicity  ignored  the  vertical  dimension  of 
the  canopy  and  used  the  single  leaf  model  to  estimate  latent  heat  loss 
from  a canopy  at  full  cover  by  replacing  single  leaf  resistances  eith 
corresponding  canopy  resistances.  Also  the  net  radiation  absorbed  by 
the  canopy  was  computed  by  subtracting  soil  heat  flux,  6,  from  the  net 
radiation  above  the  canopy.  Accordingly,  the  latent  heat  loss  from  a 
canopy  was  given  by 

a ( RN  - S ) ♦ C 0 /SB 

'•  ■ TTTTTnfSTHBT 


(3.6) 
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Hn«re  RB  and  RS  are  the  canopy  aerodynamic  and  surface  resistances  to 
heat  and  vapor  flo«  respectively.  Honteitli  computed  the  canopy 
surface  resistance  to  vapor  transport.  RS  as  Che  mean  resistances  of 
various  layers  of  leaves  acting  In  parallel.  Thus  If  r,  is  the  mean 
stomatal  resistance  (considering  top  and  hottom  surfaces  acting  in 
parallel),  and  the  canopy  with  leaf  area  Indeii  of  LAI,  RS  becomes, 

RS  - rj/LAI  (3.7) 

Eqn.  (3.7)  relates  changes  in  the  canopy  resistance  with  the  crop 
development  and  consequently  the  water  loss  through  Eqn.  (3.S).  There- 
fore Eqn.  (3.6)  could  be  applied  over  the  entire  growing  season,  pro- 
vided that  (RN-S)  is  replaced  by  energy  absorbed  by  the  canopy  only. 
However,  a similar  model  for  soil  evaporation  and  soil  resistance  was 
lacking  which  made  it  difficult  to  apply  the  enegry  balance  equations 
Co  soils  except  when  the  soil  was  saturated  or  dry,  A model  of  soil 
evaporation  using  expressions  to  characterize  soil  resistance  was 
developed  and  tested  In  cnapter  II.  Therefore  it  is  now  possible  to 
combine  lAbnteith’s  canopy  model  with  Che  soil  evaporation  model  along 
with  an  appropriate  energy  distribution  model  to  develop  an 
Integrated  soil-canopy  model  applicable  over  the  entire  range  of 
growing  season. 

Implicit  In  Eqn  (3.6)  Is  the  assumption  that  Che  entire  canopy 
acts  as  a single  isothermal  sink  for  total  absorbed  energy,  Rn-G.  On 
the  contrary  this  assianption  Is  rarely  true  except  may  be  for  canopies 
with  low  LAI,  In  reality  Che  canopy  behaves  as  multiple  sources  and 
sinks  distributed  throughout  tne  canopy  depending  upon  distribution  of 
radiation  absorption  and  stomatal  resistances.  Roughly,  tne  state  of 
leaves  In  the  canopy  can  be  described  as  either  sunlit  or  shaded  Chat 
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change  dynapnlcally  as  a function  of  sun  elevation,  row  geonetry, 
canopy  extinction  cnaracterlstics  ano  LAI,  As  a result,  only  a frac- 
tion of  LAI  is  exposed  to  direct  radiation,  while  remaining  LAI  expe- 
riences only  the  diffuse  radiation.  Since  the  anount  of  light 
Incident  on  a leaf  affects  Its  stonatal  resistance,  the  sunlit  and 
shaded  leaves  should  be  distinguished  to  obtain  more  accurate 
estimates  of  transpiration.  Sinclair  et  al,  (1976)  developed  a "big 
leaf"  model  where  the  canopy  was  mathematically  condensed  Into  a 
single  plane,  with  all  leaves  exposed  to  same  microenvironment.  They 
calculated  the  stomatal  resistance  of  the  sunlit  and  shaded  leaves 
separately  and  summed  It  with  the  boundary  layer  resistance  In 
series.  8S  was  obtained  by  sunning  In  parallel  the  total  resistances 
of  sunlit  and  shaded  leaves  weighted  by  their  respective  leaf  areas, 
More  recently  Shuttleworth  and  Wallace  (1985)  have  theoretically  shown 
how  to  Integrate  one  zone  canopy  model  of  transpiration  with  soil 
evaporation  to  develop  an  Integrated  ET  model  applicable  over  entire 
growing  period. 

The  Combined  ET  Model 

In  the  approach  presented  here,  the  canopy  was  divided  Into  sun- 
lit and  shaded  leaves  with  leaves  In  each  zone  at  different  leaf  temp- 
eratures. The  two  zone  canopy  structure  Integrated  with  the  soil 
substrate  is  presented  as  an  electrical  analog  in  Fig.  3.1.  The  model 
is  valid  for  a non-stressed  crop  idiere  the  soil  may  go  through  cycles 
of  drying  and  wetting  due  to  irrigation  or  precipitation  events.  This 
representation,  consisting  of  potential  and  curnent  sources  and 
resistors  greately  helps  setting  the  system  of  equations  and  solving 
it  for  various  energy  fluxes.  There  are  three  zones  to  be  considered: 
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txe  soil  zone,  tne  sonlft  conopy  zone  and  the  shaded  canopy  zone.  In 
each  zone  tnere  Is  a current  source  of  net  radiation  which  is  balanced 
by  the  sensible  (S)  and  latent  neat  (L)  fluxes.  The  driving  forces 
(potential  sources]  for  the  sensible  heat  flux  and  the  latent  heat 
flux  away  from  the  canopy  are  shown  as  CpT^,  and  CpCp/  y wnene  and 
e,  are  the  temperature  and  vapor  pressure  at  the  reference  height 
above  the  canopy  (usually  2m  for  agrononic  crops). 

The  microclimate  In  the  canopy  was  assumed  to  differ  from  the 
ambient  conditions  above  the  canopy.  LacX  of  this  consideration  leads 
to  underestimation  of  fluxes  as  found  by  Sinclair  et  al.  (19?6).  It 
was  assumed  that  the  aerodynamic  mixing  within  the  canopy  can  be 
sufficiently  good  so  tnat  the  microclimate  of  the  canopy  was  described 
by  one  value  of  temperature  and  humidity  and  its  location  corresponded 
to  the  mean  canopy  airstream  height  (Lemon.  1965;  Thom,  1972).  Models 
based  on  this  asstanptlon  have  been  previously  tested  (Slack  et  al., 
1970;  Sinclair  et  al..  1976).  Witnin  each  zone  the  exchange  of  sensi* 
ble  and  latent  heat  occurs  1h  response  to  the  temperature  and  vapor 
pressure  gradients  between  the  mean  canopy  air  flow  height  and  the 
mean  surface  properties  of  the  three  zones  respectively. 

To  develop  the  energy  balance  equation  It  1s  generally  assumed 
that  the  canopy  1s  homogeneous  1n  the  horizontal  direction  and  horizo- 
ntal fluxes  are  much  smaller  relative  to  the  vertical  fluxes  of  energy 
exchange  so  that  they  can  be  Ignored.  The  vertical  transfer  of  vapor 
and  heat  between  the  canopy  air  stream  and  the  air  was  proportional  to 
the  vertical  differences  in  humidity  and  temperature  and  It  was  cont- 
rolled by  a bulk  aerodynamic  resistance,  R.  Based  on  the  electrical 
analog,  and  the  principle  of  energy  balance,  the  model  shown  in  Fig. 
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3.1  can  be  nathemaffcal ly  deicribed  by  the  fundamental  equations  as 


Zone  of  sunlit  leaves. 

RNj,  . Lj,  ♦ Sj, 


Zone  of  shaded  leaves. 


Rtlq,  . L^,  ♦ S^I 


Soil  surface. 

RNjS  - S . L,s  * S55 

'=0  '''^s>■^ 

'■S5  • T '"^s  * "^Ss’ 


(3.8a) 


(3.80) 


(3.8c) 


(3.9a) 


(3.90) 


(3.9c) 


(3.10a) 

(3.100) 

(3.10c) 


Similarly  exchange  between 


canopy  airstream  beigOt 


reference  Height  1$  given  by; 


•■si  * '■dl 


(3.11a) 


(3.11b) 


Hhere  Ti  the  temperature  ( subscript  a-a1r,  c«canopy,  s1-sunl1t 
leaves.  dl*dark  leaves,  ss-soll  surface).  SB,  and  RS,  the  sunlit  and 
shaded  leaf  ensemble  boundary  layer  and  surface  resistances  to  heat 
and  vapor  transfer  (the  sum  of  RB  + RS  mil  be  denoted  by  RV,  total 
resistance  to  vapor  transfer),  R the  turbulent  eddy  diffusive  resist- 
ance to  vertical  transfer,  e(T^)  are  the  saturated  vapor  pressure  at 
temperature  Tj  and  RN,  are  the  net  radiation  absorbed  at  each  zone. 

The  soil  resistance  RSjj  in  s/t  «as  defined  in  Chapter  II  as: 


(.  r 8 RNJ  { CE  - CEJ" 


RSjj  . 154  ( e ° ’■  '■  . 1 ) (3.12) 

o ■ 1.0242  ' 10'^.  B ■ 2.0197  * 10"^  and  n ■ 1.36757  were  empirically 
determined  constants  for  the  sandy  soil  In  this  study.  RN,  was  the 
average  dally  net  radiation  (mm  of  water)  since  beginning  of  the  dry- 
ing period  but  RNj  should  not  go  below  2 mm.  CE„  was  normalized  cumu- 
lative evaporation  (in  3)  obtained  by  dividing  cianulative  evaporation 
In  nn  by  potential  cumulative  evaporation  (CEg'lR-B  nn  for  fine  sand 
Initially  at  the  field  capacity  water  content)  associated  with  last 
water  application  event,  and  CEj  was  the  critical  value  of  CE„  when 
stage  one  of  soil  evaporation  en0s(CEj.i5%  of  CEp  for  Arredando  fine 
sand).  Procedures  on  how  to  obtain  these  parameters  for  specific 


discussed  in  Chapter  II, 


D1St1nqul5h1nQ  Bttiigeii  Sunlit 


Shadtd  I.B«v85 


The  relaffve  csntrlbutlon  of  sunlit  and  shaded  leaves  to  canopy 
transpiration  depends  upon  the  LAI;  at  small  LAI  most  eater  loss  eoulO 
be  expected  from  the  sunlit  leaves  eh1le  at  full  canopy  the  shaded 
leaves  may  make  substantial  contribution.  In  practice  only  the  net 
radiation  above  the  canopy  (measured  or  computed)  and  leaf  area  index 
could  be  easily  supplied,  while  the  radiation  regime  within  the 
canopy,  Che  resistances,  as  we1l  as  sunlit  and  shaded  LAI's  are  rarely 
measured.  The  following  sections  describe  procedures  to  determine 
these  unknowns. 

The  total  LAI.  LAIj,  can  be  apportioned  Into  sunlit  LAI  (LAIj,) 
and  shaded  LAI  (LAI^^]  using  Beer's  law.  Accordingly,  the  net  radia- 
tion falling  on  a horizontal  surface  below  leaf  area  Index  of 

LAI  Is  related  to  the  net  radiation  above  the  canopy  (Rh).  and  Che 
canopy  extinction  coefficient  K; 

RNm  . RN  • exp(-  K*LAI  ) (3.13) 

The  ratio  (RNui/RH)  is  also  equal  to  the  fractional  area  of  sunlit 
leaves  on  a horizontal  surface  below  the  LAI.  Therefore  the  total 
area  of  sunlit  leaves  1n  the  canopy  Is  obtained  by  Integrating  Eqn. 
(3.14)  from  top  (LAI-0)  to  bottom  (LAI-LAI^). 

RNu^,/RN  . exp  (-K*LAI)  (3.14) 

Ulj 

LAIj,  ■ / exp  ( - K * LAI  ) d (LAI)  (3.15) 


{-K’UD) 
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= (1/K)*  (1.  - exp 

• L4It  - UIsi  (3.16b) 

Sinflarly,  tbe  net  radUttpn  captured  by  tne  sunlit  canopy  (RNj]) 
and  that  reaching  the  soil  surface  (Rhjj]  is  radially  obtained  from 
Eqn.  (3.13): 

RN,,  . RN  • (1  - exp(-K’LAl),  A (3.17) 

RN55  • RN  • exp  (-K'LAl)  - RNji  (3.18) 

where  is  the  net  radiation  received  by  the  shaded  leaves.  An 

empirical  expression  relating  RN  and  was  developed  using  neasure- 
oients  in  a fully  developed  well-watered  soybean  canopy.  Accordingly. 
RN^I  was  computed  as 

RHjI  ■ 0.05  • RN  • R • LAIfl]  (3.19) 

Leaf  Resistance 

To  apply  the  fundamental  Eqns.  (3. 1-3. 7)  to  the  sunlit  and  shaded 
regions  of  the  canopy,  individual  leaf  stomatal  and  boundary  layer 
resistances  are  required.  The  stonatal  diffusive  conductance  was 
expressed  as  a function  of  photosynthetic  photon  flux  density  (PPFO): 

C ■ C„ax  [1  - exp  ( - 5 • PPFD  /Z„„  ) (3.20) 

and  rj  • 1/C  (3.20b) 

Uhene  was  the  maximum  conductance,  and  5 the  initial  slope  of  C 
vs  PPFO.  All  the  parameters  in  Eqn.  (3.20)  are  inputs  to  the  model. 
However,  the  PPFO  on  the  shaded  leaves  is  usually  not  available.  Due 
to  1acX  of  better  Information,  similar  to  RNei.  a locally  developed 
empirical  function  to  conaiute  the  PPFD  on  shaded  leaves  was  used. 
Consequently.  PPFO  on  shaded  leaves  was  computed  as  51  of  the  above 

ifi^ 


LAIsi 

LAIdl 
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canopy  PPFD.  The  mean  sto^atal  resistance  of  the  sunlit  anO  shaded 
leaves  «as  obtained  by  substituting  their  respective  PPFO  in  Egn. 
(3.20). 

Tne  individual  leaf  boundary  layer  resistance  was  assuned  cons- 
tant at  10  s/m  in  this  model.  Under  field  conditions  it  may  vary  with 
leaf  geometry,  but  for  a plant  coimunity  as  a whole,  it  will  be  always 
under  100  s/m  and  will  probably  range  from  6 to  60  s/m  (Jackson, 
1902;  Monteith.  1973).  It  will  be  shown  later  that  the  model  Is  not 
very  sensitive  to  the  leaf  boundary  layer  resistance  and  the  value 
selected  is  an  appropriate  approximation. 

The  leaf  ensemble  resistances  of  sunlit  and  shaded  leaves  were 
obtained  using  U!,,  and  tai^,  in  Egn.  (3.7) 

Diffusion  Resistance 

In  the  model  two  additional  diffusive  resistances,  Rbss  from  the 
soil  to  the  mean  canopy  airstrean  and  R from  there  on  up  to  the  refer- 
ence height  were  introduced.  The  diffusion  resistances  depend  upon 
the  wind  speed  distribution  within  and  above  the  canopy.  According  to 
Them  (1972)  turbulent  eddy  diffusive  resistance  between  two  heights  tj 
and  i2  is  obtained  by  following  equation: 


R ■ /(1/6)  dh  (3.21) 

where  5 Is  the  eddy  diffusivlty  coefficient.  i above  a folly  deve- 
loped crop  of  height  H is  given  by, 

a(H).  U ( H - d )/  In 


[ (a  - d)/a^] 


(3.22) 


(0.4). 


d 1s  the  zero  plane  displ- 
acement height,  z Is  the  reference  height  (2  m),  and  z„  Is  the  roughn- 
ess characteristic  of  the  vegetation.  The  size  of  Zg  specifies  the 
buU  effectiveness  of  vegetation  as  a monentun  aosorber,  while  d can 
be  considered  to  Indicate  the  mean  level  at  which  momentum  is  absorbed 
by  the  plant  community.  Thus  In  crops  of  uniform  height,  (a  * Zij) 
defines  the  height  where  mean  canopy  flow  conditions  apply.  In  pract- 
ice d Is  found  by  plotting  U{z)  against  in  (z-d),  for  various  d values 
until  the  value  of  d for  which  a straight  line  plot  results  Is  the 
accepted  value  of  a.  The  can  be  then  found  from  the  Intercept  of 
this  straight  line  on  the  In(z-d)  scale.  Stanhlll  (igeg)  found  a 
relationship  between  d and  N and  Szeicz  et  a1,  (19(9)  found  a relation 
between  z^  and  M valid  for  a number  of  crops. 

4 ■ O.T  « (3.23a) 

*0  ■ 0-13  H (3.23b) 

The  wind  speed  and  the  diffusivlty  coefficient  are  assumed  to  extingu- 
ish exponentially  with  height  (h)  In  the  canopy  (Lemon,  1965;  Thom, 
1972)  and  it  Is  given  In  terms  of  the  velocity  attenuation  coefficient 

4(h)  . 4 (H)  exp  ( -n  (H  - h)  /H  ) (3.24) 

for  most  agricultural  crops  n varies  between  2 to  4 (Lemon,  1965).  9 
can  be  computed  by  Integrating  Eqn.  (3.21)  first  from  dtz,  to  H and 
then  from  H to  a (see  Fig.  3.1  for  details)  with  diffusivlty  coeffic- 
ient given  by  Eqn.  (3.22)  and  (3.24)  and  d and  Zg  by  Egn.  (3.23). 

D : dh  , dh 

“ ■ ) "aTh)  ‘ “4THT 


(3.25a) 


[ (i-a)/*,] 


[ in  [(2-d)/(H-l)  ] * -f^(C5  - 1)1  (3.25b) 


Sinilbriy,  RBjj  can  b:  obtained  as: 


R0„ 


(3.26a) 


In  [ ( 2-d)/2  1]  „ 

Mss  ■ fi  ( H - d ) C '^7  - 

-here  Cg  . exp  C n(H  - d - 2g)/H]  ana  C7  . exp  [ n(H  - 2j,)]. 

Fon  bane  soil  R and  RBjj  are  given  by  Eqns.  (3.27)  ana  (3.28). 

R "In  (2/2([)  In  (2/f*2o)/  k*k*U  (3.27) 

RBjs  ■ ln(2/2oi)  ln:(d«o)/2a.)J/  k*kU  (3.28) 

«here  2g.  «as  tbe  nougbness  height  of  tne  so11,  nomally  taken  as  the 
clod  height  on  In  flat  fine  sand  we  assumed  It  0.01  m.  In  deriving 
Eqns.  3.27  and  3.28  It  was  assuned  that  the  height  of  mean  canopy 
alrstream  remains  at  [0^2^*0.8311}  for  both  the  bare  soli  and  fully 
developed  canopy. 

The  effect  of  crop  height  at  full  canopy  on  water  balance  will  be 
further  discussed  In  the  sensitivity  analysis  section.  For  the  tine 
being  it  can  be  stated  that  for  short  crops  this  assumption  will  not 
have  much  quantitative  effect  on  the  water  balance  components.  The 
validity  of  this  assumption  especially  for  bare  soil  and  tall  crops 
may  be  questionable.  Sinclair  et  al.  (197S)  used  different  estimates 
of  mean  source/sInk  heights  for 


sensible  (d+z^wO.! 


sa 

(dtZg'fl.SH)  neat  respectively.  There  Is  no  douht  Chat  the  crop  nelght 
and  planting  density  w11t  have  some  effect  on  and  R;  hoMever  at 
this  time  there  Is  no  quantitative  information  available  Co  character- 
lae  tnis  effect.  However  if  the  Information  on  the  variation  of  mean 
source/sink  height  with  canopy  height  Is  available  one  can  use  Eqns. 
(3.aa)  and  (3.26)  to  compute  R and  R8ss.  In  the  model  presented,  R 
and  R8jj  were  computed  at  full  canopy  height  ana  bare  soil  using  user- 
specified  mean  wind  speed  and  their  values  at  Intermediate  stages  were 
computed  using  linear  Interpolation  between  the  two  limits.  Following 
the  stage  of  full  canopy,  the  source  sink  height  remained  constant  at 
0.33H,  and  Eqns.  (3.2S)  and  [3,26]  were  used  to  compute  R and  RB^j  1n 
response  to  varying  wind  speed. 

Model  Solution 

The  model  was  described  by  eleven  simultaneous  equations  [Eqns. 
3.B-3.11).  The  model  has  a total  of  twelve  unknowns:  four  tempera- 
tures [Tg],  Tjj,  Tgj,  and  Tg),  vapor  pressure  e^,  three  latent  and  3 
sensible  neat  flukes,  and  soil  heat  flux  G,  This  system  of  equations 
can  be  solved  only  If  one  of  the  twelve  unknowns  Is  known  or  using  an 
Iterative  solution.  Here  both  the  solution  scenarios  are  presented. 

First  we  win  show  how  to  eliminate  T^  and  e^  by  mathematical 
substitution  using  Che  saCuhatlon  vaoor  pressure  vs.  temperature 
relationship  [Eqn,  3.t)  to  obtain  three  simultaneous  equations  depend- 
ent on  surface  temperatures.  Substituting  for  e^  using  Eqn.  (3.1U) 
In  Eqns.  3.8b,  3.9b,  and  3.10b  and  then  using  the  definition  of  L 
from  Eqn.  3.4,  these  three  equations  can  be  written  In  matrix  form. 


S+RV, 


] CLsO  [T.,-T.] 

"[;■!:$ -'’Bj 


C RL  ][Lj]  . Cj  [T,-T,]  ♦ C3O 

Sfnce  all  the  elements  of  [ RL  ] are  known,  both  slOes  can  be  miltlol- 
led  by  the  Inverse  [IRL]  of  [RL]  to  obtain  an  eouatlon  In  L’s.  There- 


CL,]  ■ C2  CIRL]  [Ti-Tj]  ♦ Cj'O  [IRL]  (3.29a) 
Similarly  can  be  replaced  from  Epns.  3.8c,  3.9c,  ana  3.10c  using 
Eqn.  (3.11b).  to  obtain  a natrlv  for  sensible  heat. 


] CSjiJ  [T.i-T,] 

] [Sa  ] - Cl  [t!  -T*] 

] [s„]  ^ CTsl-i] 


or  C RS  )[$,]  -CiCTi-T,] 

and  CS,]  - Cj  [IRS}  [Tj-T,]  (3.29b) 
Where  Cj-Cp,  Z^K^/r  and  aCj,  [RL]  and  [RS]  are  the  matrices  of 
latent  heat  and  sensible  heat  associated  resistances,  with  Inverse 
[IRL]  and  [IRS],  [Ti*Tg]  Is  a matrix  of  temperature  differences 
between  each  of  the  three  surfaces  and  air  at  the  reference  height. 
Finally  substituting  Eqn.  (3.29a)  and  (3.29b)  In  the  energy  balance 
Eqns.  (3.8a,  3.9a.  and  3.10a,  gives 

[RN,]  . [L,]  ♦ [5,]  • (3.29c) 
Rearranging  the  terms  after  substitution  gives  the  deslrad  set  of 
three  simultaneous  equations. 


C"^s1*^a^  [RRsl  ^ 

C^dl-T,]  ■ [IC4]  [RNj,  ] - [Cg]  (3.30) 

[Tss-T,]  [RHss  - 
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Where  [lOj]  . Cj  [iRS]  ♦ Cj  C'RL]  and  [C5]  vCj'D  [IRL],  Eqn.  (3.30) 
is  a set  of  three  slnultaneoiis  equations  with  four  unknowns  S,  Tjj, 
T51  and  T^|. 

Case  1:  iCnown  Soli  Surface  Temperature 

Among  the  four  unknowns,  soil  surface  temperature  T55,  is  perhaps 
the  easiest  to  measure.  With  the  soil  temperature  known,  first  fi  can 
be  obtained  from  Eqn.  (3.30)  after  rearranging: 

■(tC4_3i*RN,,  + 1C4_32*RNj,  - Cj^jj  ♦ T,  - 

where  ICgji,  stands  for  element  in  the  j'th  row  and  k‘th  column  of 
inverted  matrix  CCj7.  A similar  explanation  follows  for  Cj  33.  Know- 
ing G,  first  the  sunlit  canopy  temperature  (Tj^)  and  Che  shaded  canopy 
temperature  (T^j)  are  obtained  from  Eqn,  (3,30),  These  temperatures 
when  substituted  in  Eqns.  (3.29a)  and  (3.29b)  give  the  latent  and 
sensible  heat  fluxes.  Finally  to  complete  the  solution,  Eqns. 
(3.20a)  and  (3.290)  are  used  with  known  fluxes  of  L^'s  and  Sj's  to 
obtain  temperature  and  vapor  pressure  e^  describing  the  mean  canopy 
microclimate.  Thus  E,  T and  ET  at  any  time  Is  obtained  by  dividing 
the  respective  latent  heat  flux  density  with  latent  neat  of  vaporlza- 

£T  . (Lj,  4 L,,  * Lj,)/  i (3.31a) 

E ■ LSS/  \ (3.31b) 

T . (LSI  4 Ldl)/  \ (3.32c) 

Cate  11:  Unknoiei  Soil  Surface  Temperature 

In  many  situations  soil  temperature  Is  not  known  and  an  Iterative 
scheme  is  needed  to  solve  Eqn.  (3.30).  In  the  iterative  scheme  first 
G would  be  guessed  to  obtain  surface  temperature  from  above  ground 
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Enn.  (3.30),  say  Since  the  above  ground  ana  below  ground  proc- 
esses are  coupled,  soil  surface  temperature  can  also  be  obtained  fron 
a heat  balance  of  soil  surface  using  a numerical  model  of  layered 
soil,  say  Hhen  these  two  temperatures  agree  with  each  other 
within  a certain  error  bound  the  convergence  criteria  Is  met. 

To  develop  the  general  purpose  numerical  sub-model  of  heat  flow, 
the  soil  was  represented  by  a nine  layers  varying  in  thickness  from 
0.1  to  0.4  m (Fig.  3.2).  The  heat  flow  rate  in  and  out  of  a layer  Is 
calculated  from  the  temperature  difference  between  adjacent  layers  and 
the  the  thermal  conductivity.  The  instantaneous  heat  flow  rates  are 
realized  over  a small  time  Interval,  and  based  on  the  principle  of 
conservation  of  energy,  the  net  change  (input-output)  in  heat  is  used 
to  obtain  the  volumetric  heat  content  of  each  layer  a finite  instant 
of  time  later.  The  calculations  are  then  repeated  for  the  next  time 
step  (De  Uit  and  van  Keulen,  1972),  It  was  assianeC  that  the  soil 
surface  temperature  varies  In  response  to  the  energy  balance  of  the 
top  soil  layer,  and  there  is  no  heat  flow  across  the  bottom  at  2.1  m 

Therefore,  If  the  soil  surface  temperature  is  T^j,  the  heat  flow 
S Into  the  top  soil  layer  of  thickness  T^,  temperature  T55  ^ and  the- 
rmal conductivity  is  given  by: 

6 ■ (^ss  - ^s,l)/Ti  ■ . 

^ss  ■ ^ss.l  f t > (3.33a) 
Initially  (at  tlme>0)  Tj,  ^ is  known  fron  the  initial  condition  and  at 
times  greater  that  zero,  it  is  again  known  from  the  previous  step 
calculations.  Therefore  from  known  8 and  soil  surface  tempera- 
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ture  can  Oe  obtained  from  below  ground  calculation  using  Eqn. 
(3.33a).  If  Tj5_3  and  are  equal  within  some  tolerance  limit 
we  have  the  right  magnitude  of  G and  all  remaining  unknowns  can  be 
obtained  as  outlined  in  the  eanlien  section  above.  Usually,  they  wi11 
be  diffenent  and  few  iterations  would  be  necessary.  The  aim  of  the 
iteration  scheme  Is  to  find  an  Increment  of  G that  will  make  above 
ground  computed  soil  temperature  ^ ^rid  T55  jj  equal.  This  can  be 
done  with  the  idea  suggested  by  Chan  (1984).  Let  Che  required  incre- 
ment of  G be  denoted  as  G*  and  the  corresponding  surface  temperature 
Tjj' . then  we  can  write; 

^ss'  ■ ^ss.a  ■*  (“'’’ss.a/'*®'  (3.33b) 

and  T„'  . r (dT^s^tj/dS)  G'  (3.33c) 

where  dT^^  ^/dG  and  dT^^  ^/dG  denotes  the  derivative  of  ^ and 
Tjj  j,  with  respect  to  G respectively.  Eliminating  T55’  from  Eqns. 

(3.33b)  and  (3.33c)  the  Increment  for  G car  be  found  as; 

i(«. 

G'  • - Tj5_^)  /(dTjj_|j/0G  -dT5,_j/d6)  (3.33d) 

Expression  for  dT^^  jj/dG  can  be  readily  obtained  from  Eqns.  (3. 333).^"^'^^ 
and  dTjj^j/dS  can  be  obtained  by  first  combining  Eqns.  (3.10a-c)  and 
then  differentiating  with  respect  to  G: 


or  dTss_,/aa  - (-  T /Cp)  R055*RV5s/(  1 1)855  + '<  "''ss'  33e) 

thus  6'  cdn  be  computed  ec  each  step  to  obtain  the  next  value  of  8 and 
continue  this  Iteration  procedure  until  the  error  between  T5J  5 and 
T55  p is  within  certain  Unit.  When  the  convergence  criterion  is  met. 
Che  soil  temperatures  of  all  layers  are  conputed  using  the  numerical 

The  iteration  scheme  outlined  above  converges  quickly.  Calcula- 
tions have  shown  that  with  the  top  layer  of  O.QI  n thickness  and  30 
minute  tine  step,  two  iterations  were  needed  to  get  convergence  within 
3%.  The  solution  nethod  presented  above  requires  repeated  matrix 
Inversion  operations  which  nay  also  be  feared  to  slow  the 
computations.  However  since  we  are  dealing  with  only  a 3x3  matrix 
this  is  not  expected  to  be  a problem  if  a proper  matrix  inversion 
method  is  used.  The  Gauss-Jordan  elimination  procedure  with  pivoting 
strategy  was  used  in  the  present  model. 

hodel  Inouts/Outputs 

The  model  inputs  were  grouped  into  three  groups;  the  weather 
conditions,  plant  specific  variaPles,  and  soil  specific  variables. 
Tne  parameters  which  remain  constant  are  indicated  by  a star  In  the 

Weather  Variables 

Weather  parameters  needed  are 
1.  Net  radiation  (RN)  in  u/a? 

Z,  Photosynthetic  photon  flux  density  (PPfO)  In  u E/m^-sec 
3.  *1r  temperature  (Tj)  1"  °C 
a.  Dew  point  temperature  In  °C 


S.  Mind 


(U]  In  n/sec 


PTint  VariaOles 

Plant  characteristics  required  are: 

1.  Lear  area  Index  (UI) 

2.  • Canopy  light  extinction  coefficient  [K) 

3.  • LAI,  at  full  canopy 

4.  • Height  at  full  canopy  (H)  in  n 

5.  * Wind  attenuation  factor  (n) 

6.  ' Haxlmum  stomatal  conductance  In  m/s 

7.  • Initial  slope  of  conductance  vs.  PPPO  {$)  In  mV  u E 
Soil  tfarlables 

Soil  characteristics: 

1.  * Constants  a,  8 and  n In  the  model 

2. *  Cnltlcal  value  of  normalized  cumulative  evaporation  (CE^) 


3.  * Potential  soil  evaporation  (CEg)  In  nm 

4.  • Thicitness  of  evaporation  zone  In  m 

5.  * Soil  roughness  height  (z„')  In  m 

6.  Thermal  conductivity  of  soil  as  a function  of  water  content  In 

J/m-sec-“C 

7.  Heat  capacity  of  Che  soli  as  a function  of  water  content  in 
J/m^-  ‘C. 

0.  •*  Number  of  soil  layers  and  their  thicknesses  In  m 

9.  * Error  bound  it  iterative  solution  is  desired 

10.  Initial  soil  temperatures  in  °C 

11.  Initial  soil  water  content  by  layer  on  volumetric  basis 

12  • Field  capacity  water  content  by  layer  on  volumetric  basis 


13.  * Lower  limit  of  soil  water  extraction  by  layer  on  volometric 

The  main  output  variables  ere  listed  here.  Other  variables  may 
be  outputed  as  per  objectives. 

1.  E,  T and  ET 

2.  Temperatures  of  sunlit,  shaded  canopy  zones 

3.  Soil  surface  temperature 

4.  hean  canopy  temperature  and  vapor  pressure 
OroanUatlon  of  The  >loee1 

The  model  was  written  In  FOBTRAN  77  and  was  organized  in  modular 
form  so  that  It  would  be  easier  to  integrate  it  with  other  models  or 
into  other  models.  Program  MAIN  is  the  driver  and  calls  other  subrou- 
tines. Convnon  blocKs  were  used  to  share  Infomation  between  various 
subroutines.  Two  output  routines,  one  daily  and  one  seasonal  are 
included  to  obtain  a variety  of  output  variables.  The  flow  chart  of 
the  model  is  presented  in  Fig,  3.3. 

The  model  begins  with  Initializing  certain  Xey  variables  such  as 
seasonal  evapotranspiration , evaporation  and  transpi ration.  The  time 
Invariant  variables  are  read  next  in  the  subroutine  FIXEOI  (FIXEO 
Input)  from  the  input  file.  One  of  the  constants  read  in  this  routine 
is  a key  for  selection  of  the  solution  method  outlined  earlier  (KET-l 
solution  for  known  soil  Surface  temperature,  KET«2  for  unknown  soil 
surface  temperature). 

At  the  beginning  of  each  day,  half  hourly  weather  data  beginning 
at  0.5  hours  for  air  temperature,  dew  point  temperature,  net 
radiation,  PPFO,  wind  speed,  soil  surface  temperature  (if  KEt*l)  is 
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read  in  the  subroutine  HEATHER.  Next  subroutine  PUNT  reads  UI  once 
euery  daji.  At  tbe  sresent  tine  subroutine  IRRIGATE  Is  called  at  user 
specified  irrigation  frequency.  It  is  assined  that  irrigation  mount 
is  adequate  to  refill  the  soil  orofile  to  its  field  capacity.  Soil 
tenperature  in  each  layer  [if  K£f»2)  is  updated  using  weighted  soil 
and  irrigation  water  temperature  by  assuming  that  irrigation  water  was 
at  air  temperature.  Cumulative  soil  evaporation  (CE)  and  cumulative 
running  average  radiation  on  soil  (RNa)  are  Initiallted  to  zero. 

During  the  day  at  each  time  step,  RN^,,  RH^,.  RN^^.  UIj,,  LAl^,. 
RVji.  RVj^,  RV^j  are  computed  in  the  MAIN  program.  Subroutine 
SOLUTION  includes  ail  equations  and  steps  described  in  the  solution 
section.  It  calls  three  other  subroutines  to  perform  Intermediate 
calculations.  Subroutine  GAUSSJ  is  a matrix  inversion  algorithm  based 
upon  the  Gauss-Jprdan  elimination  method  using  pivoting  strategy. 
Subroutines  tiATRRO,  and  HATADO  perform  matrix  multiplications  and 
additions  respectively.  Subroutine  SOLUTION  returns  L$^,  Lgg, 

^sl-  Sd|.  Sjs-  8-  ^si’  ^d1-  ^ss  T,.  and  e,.  values  to  the 

MAIN  program.  If  KETwR  soil  temperatures  are  computed  in  Che  subrou- 
tine SOILTEHP. 

The  diurnal  energy  balance  components  are  integrated  in  Che  subr- 
outine OUTPUT  (OUTPUT).  Additional  calculations  performed  by  tbe 
routine  Include: 

- Difference  between  canopy  mean  temperature  and  air  at  lOOO  n 

- Maximum  and  minimum  mean  canopy  temperature 

- Maximum  and  minimum  soil  surface  temperature 

The  user  may  need  to  modify  this  routine  to  obtain  soecific  outputs. 

Mater  lost  to  evapotranspi ration  is  used  to  update  the  soil  water 


evapo^^E1on 


assune<l  xell  «ater«d. 


water  coritent  below  the  evaporation  zone  renialns  at  the  field 
capaciti'.  Soil  evaporation  occurs  entirely  from  the  evaporation 
tone.  Transpiration  loss  nay  also  occur  from  the  evaporation  tone. 
If  so.  Che  user  Oeflnes  a table  look-up-functlon  which  includes  frac- 
tion of  transpiration  coning  from  the  evaporation  tone  as  a function 
of  LAI.  It  Is  possible  to  change  this  functional  form,  say.  to  days 
after  planting  Instead  of  LAI.  This  feature  is  Included  so  that  root- 
ing patterns  specific  to  a crop,  soil  and  management  practice  could  be 
Incorporated  In  the  model. 

At  the  end  of  Che  day,  the  output  routine  OUTPUT  is  called  to 
generate  daily  sumnary.  Then  the  nett  day's  weather  and  LAI  Is  read 
and  alt  the  calculations  are  repeated.  Ibis  process  continues  until 
Che  end  of  the  growing  season. 

%Oel  Applications 

This  model  has  a variety  cf  potential  applications  to  water 
management.  In  this  section  a brief  highlight  about  short  and  long 
term  model  predictions  Is  presented  for  an  example  of  soybean  crop 
The  soybean  crop  was  assumed  to  be  grown  In  the  constant  weather  cond- 
itions shown  in  Fig.  3. A The  model  SOYGKO  (Mllkerson  ec  al.  I9B3)  was 
used  to  generate  the  leaf  area  Index  for  Che  variety  'Bragg*. 
Irrigated  throughout  Che  season.  The  crop  reached  full  cover 
(LAI*3.S)  after  A9  days  and  attained  maximum  LAI  of  B.56  after  66  days 
In  the  growing  season.  At  full  canopy  the  height  of  the  crop  was 
assumed  D.5  m.  had  light  extinction  character! sties  of  a spherical 
leaf  angle  distribution  (K-O.B).  with  velocity  attenuation  coefficient 
n*3.  The  stomatal  conductance  data 


reported  In  Chapter  IV 
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ttie  canopy  leaf  ensemBle  resistance  RSj^  and  Using  the 
«ina  velocity  of  l.S  n/s,  R and  R8jj  were  calculated  to  be  34, S and 
82. S S/H)  for  bare  soil  and  42  s/n  and  286  s/m  at  full  canopy:  r|,.  the 
leaf  boundary  layer  resistance  was  assumed  10  s/m  and  held  constant 
throughout  the  simulation. 

Diurnal  Hodel  Predictions 

Two  diurnal  simulations  were  performed.  The  canopy  Ut[  was  2.0 
and  in  the  first  case  the  soil  surface  was  saturated  (dally 
irrigation)  whereas  in  the  second  case  it  was  dry  (no  Irrigation).  It 
was  expected  that  with  dry  soil  surface  and  no  soil  evaporation,  the 
humidity  in  the  canopy  would  be  lower  and  transpiration  should 
Increase.  The  changes  in  the  overall  energy  balance  above  the  canopy, 
in  the  canopy  and  at  the  soil  surface,  and  canopy  temperature  and 
vapor  pressure  were  compared  (Fig.  3.5)  to  show  model  predicted  inter- 
actions between  the  canopy  microclimate  and  the  canopy  and  soil. 

Energy  balance  components  when  the  soil  was  wet  indicate  that, 
the  T from  the  canopy  was  negligible  Because  of  high  leaf  resistance, 
but  soil  evaporation  occurred  due  positive  vapor  pressure  gradient. 
The  pre-dawn  soil  evaporation  cooled  the  soil  surface,  lowered  vapor 
pressure  gradient,  and  soil  evaporation  got  progressively  smaller. 
Uith  a rise  in  the  net  radiation,  and  air  temperature,  both  soil  eva- 
poration and  transpiration  increased  simultaneously.  For  most  of  the 
day,  the  net  radiation  Intercepted  by  the  canopy  exceeded  energy  diss- 
ipated by  transpiring  water.  Due  to  the  wet  surface,  soil  surface 
temperature  remained  cooler  chan  the  mean  canopy  temperature  ( 
negative  sensible  heat  transfer  during  the  day  time),  and  soil  evapo- 
ration was  low.  The  evaporation  occurred  mainly  in  response  to  the 
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vjpnr  pressure  gradient.  Due  to  low  soil  temperature  early  in  the  day, 
the  excess  ret  radiation  was  stored  as  soil  heat.  As  the  net  radia- 
tion and  air  temperature  decreased,  there  was  a insersion  In  the  soil 
neat  flux.  The  transpiration  dropped  rapidly  witn  net  radiation. 
However  due  to  slow  cooling  of  tne  soil,  soil  evaporation  decreased 
very  gradually.  For  the  24  hours  of  simulation.  In  terms  of  amount  of 
evaporated  water,  the  canopy  intercepted  an  equivalent  of  2,74  nm  of 
net  radiation,  had  transpiration  of  2, Hi  ran  and  sensinie  heat  of 
-O.lmn.  Similarly,  at  the  soil  surface  1.31  ran  of  net  radiation  was 
distributed  among  evaporation  1.20  nun,  sensible  heat  0.16  nm,  and  soil 
neat  flux  -0.05  rail. 

The  energy  balance  above  the  canopy  suimiarlaes  energy  distribu- 
tion within  the  canopy  and  at  the  soil  surface.  In  sumary,  most  of 
the  energy  available  (4.05  ran)  was  expended  for  ET  (4.04  ran)  and  a 
negligible  fraction  went  into  the  sensible  neat.  Therefore  over  a 
day,  the  heat  flux  into  the  soil  was  practically  zero. 

The  energy  balance  with  dry  soil  surface  was  similar,  Oue  to 
zero  soil  evaporation,  soil  surface  was  wanner  than  the  canopy  air  ( 
positive  sensible  heat  flux  ].  Also  the  heat  flux  into  the  soil  was 
higher;  0.27  ram  compared  to  -O.OS  ran  with  the  wet  surface.  The  ET  was 
low  and  concurrently  the  humidity  was  low  (Fig  3.6b)  while  canopy 
temperature  was  higher  [Fig  3.6a).  The  ET  was  reduced  from  4.04  ran  to 
3.11  (23%  lower),  however  the  transpiration  was  higher  by  about  ID  % 
(3.11  ran  vs.  2.84  ran)  due  to  lower  hranidity.  The  combination  of  high 
canopy  temperatures  and  low  humidity  caused  slight  increase  in  the 
sunlit  leaf  temperatures  (Fig  3.7a).  The  main  effect  of  wet  versus 
dry  Surface  was  visible  on  tne  daily  course  of  soil  surface 


Coriparison  of  sinulated  canopy  air  canperature 
and  canopy  vapor  prasaure  for  wet  (solid  line) 
and  dry  (doccad  line)  soil  surface  conditions. 
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Fig. 


ConparUon  of  sun-lit  leaf  tempenature  and  soil 
surface  temperature  for  net  (solid  line)  and  dry 
(dotted  Hne)  soil  surface  conditions 
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temperature.  In  general  with  ary  sail,  temperatures  were  5 to  6 
degrees  higher  than  with  wet  soil  surface. 

Long  Term  ET  Predictions 

Using  the  leaf  area  development  data  presented  in  Fig.  3.4c.  the 
mode!  was  used  to  simulate  one  growing  season  of  soybean.  In  the 
first  case  soybeans  were  Irrigated  dally  and  In  the  second  case 
weekly.  During  this  simulation  it  was  assumed  that  there  were  no 
roots  In  the  evaporation  zone  so  that  it  was  depleted  by  soil  evapora- 
tion alone.  The  consequence  of  root  water  extraction  from  the  evapo- 
ration zone  1s  discussec  later. 

In  the  dally  irrigated  experiment,  the  dally  ET  rates  were  almost 
constant  until  119  days  after  planting  fOAPL),  tdiich  was  the  day  of 
crop  maturity  (Fig.  3.3a).  Irrigation  was  terminated  at  crop 
maturity.  Soil  evaporation  decreased  gradually  as  the  canopy 
developed  and  less  energy  was  available  at  the  soil  surface.  ET 
during  the  first  3?  days  was  dominated  by  soil  evaporation. 
Thereafter.  ET  continued  to  decrease  and  finally  became  constant  when 
LAI  exceeded  5.  Following  crop  maturity  soil  evaporation  increased  as 
LAI  dropped  and  more  radiation  reached  the  soil  surface,  hut  soon 
this  trend  was  reversed  as  the  soil  became  dry. 

In  the  weekly  irrigation  cycle,  the  soil  surface  went  through 
cycles  of  welting  and  drying,  as  Illustrated  hy  the  varying  amplitudes 
of  C and  ET  In  Fig.  3.8b.  It  was  interesting  to  see  that  initially, 
ET  on  the  irrigated  days  was  higher  compared  to  the  daily  irrigation 
cases.  This  can  be  eiplalned  by  ekamining  the  variations  In  the  soil 
surface  temperature  shown  in  Fig.  3.9.  With  progression  1n  drying, 
the  soil  surface  temperature  Increased  much  more  than  the  daily  Irrig- 


Fig.  3.8 


Slnulatefl  rially  ET,  T arid  E from  planting  to 
harvest  for  the  dally  and  weekly  Irrigation 


Daily  simulated  Taxinum  and  rninlnun  soil  surface 
temperature  under  the  soybean  canopy  with  daily 
and  »ee«ly  irrigation  cycles 


ited  case.  Following  Irrigation,  t«e  soil  surface  was  cooled,  nut 
remained  Mgner  than  the  dallji  Irrigation  case.  The  higher  soil  t«np- 
erature  coupled  with  saturated  soil  resulted  in  higher  vapor  pressure 
gradient  that  caused  more  evaporation.  It  was  also  interesting  to 
note  the  lowering  of  transpiration  rate  following  an  irrigation 
("dips"  In  T curve),  due  to  higher  soil  evaporation  and  humidity.  Bie 
Interaction  between  E and  T had  mutually  ctsmpensating  effect  on  ET. 
The  feedback  of  E on  T became  less  significant  at  lower  E magnitudes. 

The  seasonal  water  use  with  the  weekly  irrigation  cycle  was  7i 
lower  than  the  daily  irrigation  practice.  For  e«ample.  seasonal  E 
with  weekly  Irrigation  was  120.4  mtn  compared  to  163.5  nin  with  the 
dally  cycle,  <dille  the  seasonal  transpiration  loss  382, a nrni  was 
slightly  higher  than  377.6  im  in  case  of  the  dally  Irrigation.  For 
the  daily  irrigation  practice  It  took  65  days  for  the  cumulative  T to 
equal  the  cumulative  soil  evaporation  in  contrast  to  51  days  with 
weekly  Irrigation. 

Soil  Temperature 

Simulated  maximum  and  minimum  soil  surface  temperatures  under  the 
two  irrigation  frequencies  are  compared  in  Fig.  3.9.  For  the  hare 
soil  the  maklraum  temperature  varied  from  30  to  45  degrees.  Similar 
temperature  variation  was  observed  in  the  bare  soil  evaporation 
studies.  With  canopy  develOfmient,  the  amplitude  reduced  progressively 
and  at  LAI  exceeding  2.3  it  remained  practically  constant.  On  the 
ocher  hand  the  minimum  surface  soil  temperature  varied  by  less  than  5 
degrees  during  the  incomplete  canopy  stage.  Ac  nigh  LAI  following  an 
irrigation,  the  soil  was  cooled  temporarily  ("dips"  in  the  minimum 
temperature).  The  cooler  soil  produced  low  soil  evaporation  as  seen 
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1n  the  Fig  3.8b. 

Canopy  • Mr  Twipenture 

The  temperature  at  the  mean  canopy  source/sink  height  «as  assumed 
to  be  representative  of  the  mean  canopy  temperature.  The  difference 
between  canopy  and  air  temperatures  computed  at  1000  and  1200  h for 
the  two  Irrigation  are  shown  1n  F1g.  3.10.  Since  the  crop  was  well- 
watered  the  differences  in  the  temperature  at  two  tines  was  mainly  due 
to  the  change  in  the  air  temperature  relative  to  the  canopy  tempera- 


the  crops  (Ehrler,  1S73).  These  figures  reveal  that  diurnal  patterns 
may  occur  even  in  well  watered  crops. 

In  the  earlier  discussion  it  was  assumed  that  there  were  no  roots 
in  the  evaporation  zone.  This  may  not  He  true  in  reality.  Simula- 
tions were  performed  by  assuming  that  0 , 2S,  50  and  100  i of  transpi- 
ration was  derived  from  the  evaporation  zone.  The  results  of  seasonal 
water  demand  under  these  situations  are  sumartzed  below  in  Table 
3.1.  hs  the  water  was  applied  every  day,  the  extraction  fraction  had 
little  effect  on  the  water  demand  of  the  daily  irrigated  treatment. 
With  weekly  irrigation  cycle,  increasing  extraction  reduced  soil  evap- 
oration due  to  quicker  drying.  On  Che  whole,  as  long  as  the  crop  is 
able  to  extract  adequate  amount  of  water  from  deeper  depths  (except 
the  loot  case  above),  these  results  suggest  that  computed  water  use  of 
the  crop  would  be  little  affected  Oy  this  assunption. 

Sensitivity  Analysis 

In  this  section  model  sensitivity  to  various  input  parameters 
listed  in  the  Teble  3.2  is  exsmined.  For  the  sensitivity  analysts 


Comparison  of  canopy  to  air  tamperatore 
dlffarenca  racorded  at  1000  and  1?00  n 
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TAB1.E  3.1 


Influence  of  amount  of  transpiration  occurring  from  the 
evaporation  zone  on  the  seasonal  water  use  under  daily  and 
weekly  irrigation  cycle. 


0 163.5  377.6 
as  160.9  373.5 
50  160.3  379.1 
100  161.6  379.6 


541.1  120.4  332.8  503.3 
539.4  111.6  335.1  496.7 
539.4  104.9  386.4  491.3 

541.2  95. 5 388.1  483.6 
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daily  ET  rates  were  first  simulated  at  UU  from  1 to  5 with  the 
standard  values  listed  in  Tahle  3.2  and  then  one  parameter  at  a time 
was  varied  Dy  t SO  1 of  its  original  value  to  compute  new  ET,  The 
sensitivity  was  calculated  as  fraction  change  in  ET  with  respect  to 
the  standard  ET  estimate  and  the  values  in  percentage  are  sunnarlted 
in  the  Tahle  3.3  for  the  daily  irrigated  treatment  and  in  Table  3.4 
ton  the  weekly  irrigated  treatment.  From  the  tables  it  can  be  seen 
that  the  effect  of  rb  and  n.  the  two  most  uncertain  variables,  was 
less  than  lOS.  The  model  is  more  sensitive  to  H end  K which  can 
be  measured  or  obtained  to  reduce  the  model  uncertainty. 

Sunsnary  And  Conclusions 

The  development  and  analysis  of  a mechanistic  three  aone  soil- 
canopy-atmosphere  interaction  model  was  presented.  This  model  provi- 
des a useful  tool  to  study  the  water  use  of  a crop  under  the  Influence 
of  climatic,  agronomic  and  management  practices.  Evaporation  from  tne 
soil  may  constitute  a significant  fraction  of  total  ET  over  a 
season.  The  amount  of  soil  evaporation  varies  with  crop  growth  stage, 
available  soil  water  as  welt  as  eet  radiation  input.  It  was  observed 
that  high  soil  evapbratioh  would  reduce  transpiration.  The  model 
predicts  canopy  and  soil  temperatures  and  the  simulation  showed  that 
they  could  differ  signi ficently  from  the  air  temperature  depending 
upon  the  soil  surface  water  content,  weather  and  stage  of  crop 
growth.  Knowledge  of  these  temperatures  could  be  useful  in  Irrigation 
scheduling  and  predicting  possible  pest  attacks. 

The  model  presented  here  differed  from  the  works  of  Waggoner  and 
heifsnyder  (1968),  Stewart  (1970).  Shawcroft  et  a1.  (1974),  and 
Goudriaan  (1977)  where  the  soil  and  the  canopy  was  stratified  into  a 
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nimSer  of  horizontal  layers.  In  the  oresent  moOel  the  canopy  was 
conceiveO  Co  consist  of  two  canopy  zones,  the  siinitt  leaves  ami  the 
shaded  leaves,  regardless  of  tnelr  position  in  the  canopy.  The  soil 
evaporation  was  treated  using  the  model  developed  in  Chapter  II.  Two 
types  of  solution  algorithms  were  presented  based  on  whether  soil 
temperature  was  a known  input  or  a desired  output.  With  soil 
temperature  known,  an  analytical  solution  was  developed  to  predict 
less  coumonly  measured  canopy  temperatures,  humidity  and  energy 
fluxes.  However,  a rapidly  converging  iterative  solution  scheme  was 
developed  for  situations  when  soil  temperature  was  a unknown. 

The  model  has  a practical  utility  in  Irrigation  scheduling,  pred- 
icting ET.  E and  T under  various  types  of  environments  as  well  as 
estimating  crop  microclimate.  With  advances  in  computer  technology  a 
farmer  can  install  an  automatic  weather  reporting  system  for  less  than 
the  cost  of  a garden  tractor.  Such  self-sustained  systems  are  already 
providing  weather  and  soil  information  necessary  for  the  model.  Crop 
growth  models  such  as  SOYGRO  are  becoming  widely  available  to  provide 
crop  specific  information.  Therefore  this  model  could  help  provide 
up  to  date  information  relative  to  crop  water  use  and  its  management. 

The  model  can  bo  easily  Integrated  into  crop  growth  models  like 
SDYGRO  or  can  be  used  as  stand  alone  version.  Further  work  is 
required  to  validate  the  model  under  field  conditions.  For  general 
application  of  this  model,  additional  work  is  required  to  integrate  a 
companion  soil  water  and  a crop  growth  yield  prediction  model  to 
assess  water  use  and  efficiencies  under  conditions  of  limited  water 


supply. 


CHAPTER  IV 

EFFECT  OF  ATHOSPHERIC  COo  C0NCEHTRAT10N5  OH  SOT0EAN  WATER  USE 
SIHULATIOH  ANir  ESPERIMENTAL  VERIFICATION 

IfitrajueUon 

The  mean  atnosoherlc  CO;  cencentratlan  Is  steAdUy  increasing  Sue 
to  the  release  of  COg  Oy  combustian  of  fossil  fuel.  Forecasts  suggest 
that  atnospheric  CO;  levels  win  have  doubled  by  the  end  of  next  cent- 
ury from  the  ore-industrial  concentration  (Clark,  l9BEb).  It  is  a 
classic  observation  Chat  aerial  CO2  is  necessary  for  olant  growth  and 
that  CO^  enrlchatent  usually  has  a siiiulaCing  effect  on  growth  rate. 
Increased  COj  will  increase  net  CO^  fixation,  leading  to  increased 
olant  dry  «ight  and  probably  increased  leaf  area  in  o'ants  with  C3 
Dhotbsynthesis  (Kramer,  1931).  On  the  other  hand  high  CO;  will  incre- 
ase stomatal  resistance,  reducing  transoiration  and  therefore  increas- 
ing water  use  efficiency.  These  two  effects  of  high  CO;  on  leaf  area 
per  plant  and  on  transpiration  per  unit  leaf  area  may  have  compensat- 
ing effects  Oh  water  use  per  plant. 

Another  anticipated  effect  of  atmospneric  00;  enrictkiienc  is  a 1 
to  5 *C  increase  in  the  global  temperature  with  each  doubling  of  the 
atmospheric  CO^  concentration,  a consequent  decrease  in  precioitation 
and  an  increase  in  aridity  in  many  important  areas  of  the  world  ( 
Horlson  and  Slfford,  1984).  Such  climatic  alterations  will  also  aff- 

There  have  been  many  Investigations  of  plant  responses  to  CO^  and 
atmospneric  conditions.  McNaughton  and  Jarvis  (1983)  emphasized  that 
in  short  vegetation  such  as  arable  crops,  vapor  pressure  deficit  in 
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and  aDowe  tha  plant  canopies  is  influenced  by,  and  partly  depends  on 
tne  evaporation  loss  from  the  canopy.  Thus  a reduction  in  plant  eva- 
poration due  to  Increased  COj  may  increase  vapor  pressure  deficit  that 
will  drive  more  evaporation.  This  feed  bac«  mechanism,  if  It  exists 
ln  most  crops,  could  undermine  the  effect  of  Increased  CO^  on  lowering 
crop  water  use.  Jones  et  a1.  (1905)  experimentally  showed  that  for 

soybeans  grown  in  COg  concentration  of  000  ppm,  increasing  the  ambient 
temperature  from  28  to  35  “C  resulted  In  25  t Increase  in  daily  trans- 
piration while  daily  photosynthesis  was  reduced  by  4 i.  This  apparent 
increase  in  water  use  reduced  the  water  use  efficiency  by  22t. 
Gifford  (1979)  found  that  relative  yield  of  water  stressed  wheat  grown 
In  twice  the  normal  environmental  OOg  concentrations  was  as  large  as 
well  watered  wheat  grown  in  normal  COg  concentrations.  However  the 
total  water  loss  was  similar. 

The  actual  long  term  water  use  of  COg  enriched  crops  under  field 
conditions  needs  to  be  Investigated.  Plants  grown  in  high  COg  envi- 
ronments will  nave  high  leaf  area  which  tends  to  Increase  transpirat- 
ion. They  will  also  have  a larger,  more  vigorous  root  system  and  be 
able  to  extract  more  water,  resulting  in  less  water  stress,  higher 
yields,  but  possible  higher  total  water  use  (Kimhall  and  Idso, 
1983),  Any  decrease  In  transpiration  will  make  more  energy  available 
for  soil  evaporation.  Mechanistic  models  of  soil-plant-atmosphere  can 
be  used  to  study  tne  individual  and  combined  effect  of  environmental 
conditions.  Therefore,  tne  objectives  of  tnis  paper  were  to; 

1.  Evaluate  the  ability  of  a mechanistic  evapotranspiration 
mode!  to  predict  water  use  by  soybeans  grown  under 
different  COg  environments  in  outdoor,  controlled  growth 
chambers,  and 


2.  Apply  the  moPel  to  estloiate  the  indlviPual  and  conpined 
effects  of  COj  induced  cnanges  in  stonetai  resistance, 
leaf  area  development,  and  temperature  on  water  use  by 
soybeans  in  the  field. 

Description  Of  The  tindel 

The  schematic  of  the  soil-canopy  model  developed  in  Chapter  II  Is 
shown  In  Fig.  4.1.  The  canopy  was  divided  Into  tones  of  sunlit  and 
shaded  leaves  with  soil  forming  the  lower  boundary.  Hie  division  of 
canopy  Into  sunlit  and  shaded  zones  was  based  on  the  light  extinction 
characteristics  of  the  canopy  and  leaf  area  index.  The  resistances 
considered  were  the  boundary  layer  resistances  (DB),  the  surface  resi- 
stances (RS),  and  Che  vertical  aerodynamic  diffusion  resistance  R. 
The  subscripts  si,  01,  and  ss  were  used  to  denote  the  canopy  sunlit 
zone,  the  shaded  zone,  and  the  soil  surface.  The  canopy  microclimate 
differed  from  ambient  conditions  above  the  canopy  and  was  describeO 
with  an  average  temperature  and  vapor  pressure.  Similarly  properties 
of  buH  air  above  the  canopy  were  represented  by  air  temperature  T^. 
and  vapor  pressure  e^.  The  dez^  denoted  the  height  of  momentum  trans- 
fer from  the  air  to  the  canopy.  The  model  was  described  using  a 
steady  state  approach  and  principles  of  energy  balance  by  a series  of 
simultaneous  equations.  Two  types  of  solution  algorithms  were 
presented  based  on  whether  soli  surface  temperature  was  known  or 
unknown.  The  climatic  input  to  the  model  included  above  canopy  tempe- 
ratures, nimiidlty  or  dew  point  temperatures,  net  radiation,  and  wind 
speed.  Besides  tnese,  leaf  area  Index,  stomatal  resistance  and  canopy 
light  extinction  characteristics  were  also  needed.  The  model  predicts 
diomal  fluxes  of  latent  and  sensible  heat  fnom  the  soil,  the  shaded 


l?fl 


Fig.  4.1 


Schenidtic  of  the  evapotransolratfon  nodel  for 
incomplete  end  complete  canopies  presented  as  an 
electrical  analog 
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canopy  an<1  tne  &un1H  canopy,  along  with  their  surface  temperatures 
and  mean  canopy  humidity  and  temperature. 

Studies  were  conducted  between  September  and  November  of  1934  at 
the  Irrigation  Besearcn  and  Education  Pare  of  the  University  of 
Florida  In  Gainesville,  Florida  to  obtain  model  Inputs  and  a data  set 
to  validate  the  model. 

Soybeans  [ Glycine  max(L.)  Herr.  cv.  Bragg]  were  grown  In  COj 
concentrations  of  330.  360  and  990  ppm  ( Pants  Pen  HlHlon.  mole  fra- 
ction basis,  or  moles  of  CO2  per  mole  of  atmospheric  air  ) In  natural- 
ly sunlit  controlled  environment  chambers  (OoneS  et  a1.  1934b)  located 
In  the  field.  In  this  study  three  levels  of  COj  were  considered:  330, 
660,  and  990  ppm,  respectively.  The  daytime  dry  bulb  and  dew  point 
temperature  controls  were  set  at  31  and  13  and  were  maintained 
within  1 - 1.6  *C.  These  experimental  conditions  were  continuously 
monitored  by  a dedicated  computer  operating  1n  real  time.  The  envi- 
ronmental data  were  oolleoted  every  6 minutes  throughout  the  growing 
season  and  consisted  of  dry  bulb  temperature,  dew  point  temperature, 
soil  surface  temperature,  net  radiation,  photosynthetIcaTly  active 
photon  flux  density  (PPFB).  and  COj  concentrations.  Canopy  transpira- 
tion rates  were  measured  every  6 minutes  by  collecting  water  chat 
condensed  on  cooling  colls  which  were  used  to  control  a1r  humidity. 

Bragg  soybeans  were  planted  in  North-South  oriented  rows  with  0.5 
m row  spacing  on  14  September  1934,  Because  of  late  planting,  four 
Incandescent  lights  were  placed  above  each  chamber  to  extend  phptoper- 
lod  and  prevent  premature  flowering.  As  the  canopy  developed  shade 
cloths  were  maintained  at  the  canopy  height  to  simulate  natural  canopy 


exposure  in  Incident  radletion,  T>ie  soybeans  vere  Kept  well  watered 
during  tne  entire  study  period.  Irrigation  was  frequent  and  on  demand 
during  first  two  weeks  after  pianting,  however  thereafter  a fixed 
irrigation  schedule  ( irrigation  on  Monday  and  Thursday  at  1600)  was 
followed  with  occasional  supplemental  Irrigation  on  Sunday  evenings. 
The  soybeans  were  hand  watered  until  Oct  11  (27  DRPL)  and  thereafter  a 
trickle  Irrigation  system  was  used.  This  change  was  necessary  as  Che 
wind  was  causing  lodging  when  the  acrylic  cops  were  removed  during 
irrigation.  Starting  on  Oct  3,  during  each  subsequent  week,  the  soil 
was  covered  with  polythene  sheets  beginning  Wednesday  morning  [800 
EST)  until  just  before  Irrigation  on  Thursday  (1600  EST)  to  eliminate 
soil  evaporation.  This  sequence  was  repeated  over  four  consKutIve 
weeks  (last  sequence  ended  on  Oct  24)  and  thereafter  the  soil  surface 
was  completely  sealed  off  for  ZOf  exchange  measurements.  8y  this  time 
most  treatments  had  reached  full  cover  with  UI  of  5.4g,  4.53,  and 
6.54  for  COg  treatments  330.  660,  and  990  ppm.  respectively.  Data 
collection  for  this  study  continued  thereafter  until  Hov  7. 

Curing  the  season  the  plant  height  and  width  were  measured  on  13, 
22,  26.  31.  34,  39  and  46  OAPl  and  LAI  was  measured  on  13,  24.  34,  46. 
and  66  OAPL,  Stomata!  resistances  of  Individual  sunlit  and  shaded 
leaves  were  measured  on  21,  28,  31  and  33  OAPL  using  a LI*C0R  1600 
steady  state  porometer.  Measurements  were  made  on  both  the  abaxial 
and  adaxial  leaf  surfaces  of  three  randomly  selected  sunlit  and  shaded 
leaves.  Ambient  canopy  air  temperature,  relative  humidity,  sunlit  and 
shaded  leaf  temperatures  were  also 
dieted  estimates. 


recorded  to  verify 
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Results  Of  Expgriiieital  Study 

The  time  cojrse  of  the  cunijletlve  LAI  for  the  three  CO2  treatreen- 
ts  Is  shown  in  FI9.  4.?.  The  LAI  fncreased  rapidly  In  the  990  treat- 
ment whereas  the  rates  of  development  were  almost  same  In  the  330  and 
660  treatment  for  the  first  36  OAPL.  There  was  some  ohewino  Insect 
damage  early  In  the  660  ppm  chamher  which  may  have  delayed  LAI  develo- 
pment. The  naxlmem  LAI  reached  by  the  three  treatments  were  6.61, 
6.38  and  6. 73,  respectively. 

Stomatal  resistances  of  Individual  leaves  were  calculated  as  the 
reciprocal  of  the  sum  of  the  reciprocals  of  the  diffusive  resistance 
for  both  surfaces.  A nan-1lnear  relationship  was  developed  between 
leaf  conductance  and  PPFO  for  Individual  COj  treatments.  There  was  a 
lot  of  scatter  In  the  measured  values  especially  for  the  660  OOj  trea- 
tment. The  resistances  used  in  simulations  are  listed  1n  Table  A.I. 
Hodel  Calibration 

The  aerodynamic  boundary  layer  resistance  at  the  soil  and 

the  leaf  surface  (33,]  and  PBd))  ■ and  the  vertical  diffusion  resista- 
nce (3)  shown  In  Fig.  4,1  were  related  by  the  classical  theory  to  the 
logarithmic  wind  speed  profile  above  the  canopy  and  Its  exponential 
extinction  within  the  canopy.  Oue  to  the  air  drculatloh  patterns  and 
size  of  the  chambers,  these  theories  could  not  be  applied  Co  compute 
these  resistances.  Therefore  the  model  was  calibrated  using  selected 
days  (Table  4.2)  to  characterize  the  variation  of  these  resistances 
with  the  leaf  area  development.  The  leaf  boundary  layer  resistance 
was  assumed  constant  In  this  study.  Under  natural  field  conditions  It 
should  be  under  100  s/m  and  will  probably  range  From  6 to  50  s/m 
(Jackson.  19B2i  hontelth,  1963).  The  values 


reported  by  Oeanmead 
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Fig. 


Ledf  area  developaient  as  a function  of  days 
after  planting  for  Bragg  soybean  gro*m  in 
controlled  environment  chamber  under  three 
levels  of  COj  concentrations 


12S 


TABLE  4.1  Psramsters  In  the  stomatal  conductance  function. 

S BPFD 

C . ( 1 - e ) 

•here  1n  n/s  , PPfD  1n  E/n^-sec  and  s has  units  of  r?t  E. 


TABLE  4.2 


Tne  results  of  calibration  for  the  vertical  aerodjmanic 
resistance  R and  the  soil  resistance  in  the  control 


RSss.  s/m 


Oct  10  26  1 49  10 

Oct  17  33  2.41  15 


Seo  14  0 0 153.8 
Oct  6 22  1.06  200.0 
Oct  12  28  1.7S  225.0 
Oct  20  36  3.03  300.0 
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[1976)  rangeil  Fron  5-7  s/m  for  tenpersce  cereals.  In  the  chamber  the 
air  was  continuously  circulated  at  a constant  rate  so  a constant  value 
of  10  s/n  on  a leaf  basis  was  asscwied.  The  bare  soil  boundary  layer 
resistance  value  of  154  s/t  reported  in  Chapter  II  was  used  as  the 
initial  value  for  RBjj  before  LAI  started  developing. 

Uith  these  initial  selections  the  resistances  R and  95^^ 
calibrated  using  three  days  during  early  development  (Table  4.2)  To 
calibrate  R »e  chose  days  from  the  330  ppm  treatment  when  the  soil 
surface  was  covered  and  varied  R until  measured  and  simulated  values 
agreed  reasonably  well  for  most  of  the  day.  Due  to  scatter  in  the 
stomatal  resistance  measurements  the  of  0.0190  m/sec  consistently 
gave  better  fit  than  the  estimate  (0.0208)  made  using  the  steady  state 
ponomater  values.  Therefore  this  value  was  used  in  the  simulations. 
RBjj  was  calibrated  using  days  when  both  soil  evaporation  and  transpi- 
ration occurred  simultaneously.  Days  selected  for  calibration  of  bgtn 
R and  R855  were  within  one  or  two  days  of  each  ocher  depending  on  the 
data  availability.  The  dallbraCed  values  of  R and  R8jj  varied  until 
LAI  became  3.03.  Thereafter  Che  chambers  floors  were  closed  permanent- 
ly and  no  further  calibration  for  RSj,  was  needed  and  the  value  of  R 
and  RB;,  was  assumed  constant  there  after.  No  adjustment  In  the  leaf 
boundary  resistances  were  made.  The  sensitivity  analysis  in  Chapter 
ill  had  demonstrated  that  SOI  change  in  leaf  boundary  layer  resistanc- 
es would  result  in  less  than  1%  change  in  daily  ET  from  a well  watered 
canopy.  The  results  of  calibration  are  suimarized  in  Table  4.2. 
Simulation  of  diurnal  ET 

Data  collected  during  the  first  two  weeks  after  planting  were  not 
considered  for  simulation  due  to  irregular  irrigation  schedules. 
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Selected  days  fteglnnfng  wltn  11  days  after  planting  (OAPL)  until  58 
WPL  were  simulated.  On  many  occasions  during  the  study,  the  chambers 
were  opened  for  plant  measurements,  irrigation,  spraying,  and  emergen- 
cy maintenance.  These  days  as  well  as  days  when  there  were  control 
and  instrumentation  problems  were  omitted  from  the  analysis.  The 
instantaneous  net  radiation  data  were  available  from  700  to  1900  hours 
each  day  and  therefore  all  the  results  reported  here  are  for  this 
time  period  only.  Hereafter  , the  three  COj  treatments  will  be  ref- 
erred simply  as  330.  660.  and  990  treatments,  respectively. 

The  individual  00^  treatments  as  well  as  the  time  course  of  simu- 
lated and  measured  evapotranspi ration  {ET)  during  early  (19  OIPL.  Oct 
9).  middle  (28  OAPL,  Oct  121  , and  at  near  full  canopy  (35  OAPL,  Oct 
201  are  compared  in  Figs.  1.3-4, 5.  On  Oct  3,  the  soil  was  covered 
with  plastic  sheets  to  measure  transpiration  (T)  only,  >diile  on  Oct  12 
and  Oct  20  both  soil  evaporation  (El  and  T were  occurring  simulta- 
neously. The  results  clearly  show  the  differences  in  total  water  loss 
due  td  the  CO^  treatments  and  in  general  the  water  use  was  reduced 
with  increased  OO^  concentrations. 

It  was  ^served  that  simulated  values  consistently  exceeded  mea- 
surements in  the  mornings  but  were  similar  during  the  rest  of  the  day. 
There  were  three  probable  explanations:  First,  transpiration  was 
measured  by  condensing  water  vapor  from  outgoing  air  that  passed 
through  each  chamber  cooling  coils.  In  the  mornings,  the  coils  were 
known  to  accumulate  condensate  fcr  a certain  amount  of  tine  before 
water  was  released  td  be  measured.  As  a result  the  measurements  lagg- 
ed the  true  water  loss  rate  from  the  canopy.  This  problem  was  more 
noticeable  especially  when  the  canopy  was  small.  Secondly,  1t  was 
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Fig. 


Hally  time  course  of  simulsted  anO  measured 
transolration  races  or  Oct  3 (19  OAPL)  for  the 
330,  660,  arid  990  (left  to  right)  COj 
treatments.  The  LAI  »ere  Q.Rl,  0.85,  and  l.OB 
respectively  and  the  soil  was  covered 
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Diilj'  tine  course  of  simulated  and  measured 
e*apotratisp1ratloe  rates  on  Oct  12  (28  OAPl.)  for 
tde  330.  660.  ariO  990  (left  to  right)  CO? 
treatnents.  Tne  L8I  *ere  1.75.  1.73,  ana  2.40 
respectively 
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noticed  that  on  cool  tiornlrigs  water  was  condensing  on  the  side  walls 
which  would  reduce  water  vapor  passing  through  the  cooling  coils. 
Thirdly,  the  chanbers  tended  to  loose  sensible  energy  to  the  envi- 
roraient  in  the  mornings  and  gain  sensible  energy  in  tne  afternoons. 
This  would  tend  to  decrease  ET  in  the  morning  and  increase  ET  In  the 
sfternoonn.  Pinally,  at  each  tinte  step  in  the  simulation,  the  model 
assimes  steady  state  and  point  values  of  net  radiation  and  PPED  were 
used  to  compute  energy  balance  components.  The  measurement  system  was 
responding  to  average  conditions  over  the  whole  time  step  (IS 
minutes).  This  was  particularly  important  during  mornings  when  radia- 
tion was  Increasing  rapidly.  This  is  an  important  issue  In  using 
steady  state  models  and  the  errors  could  be  reduced  to  some  extent  Py 
using  inputs  averaged  over  the  whole  time  step. 

Simulated  and  measured  values  agreed  well  during  noon  hours  when 
none  of  the  system  related  prpolems  as  well  as  the  assumption  of 
steady  state  was  approximately  correct.  On  days  when  there  was  varia- 
ble cloud  cover,  the  model,  again  due  to  the  assmption  of  steady 
state,  tended  to  underestimate  ET  (see  fig.  4.5}  late  In  the 
afternoon.  However,  on  the  clear  days,  even  the  late  afternoon  errors 
were  small.  On  Oct  20.  which  was  the  third  day  after  Irrigation,  the 
model  predicted  condensation  on  soil  surface  in  the  the  990  treatment. 
Simulation  of  Canopy  Temperatures 

The  model  also  predicted  temperatures  of  the  sunlit,  and  shaded 
canopy  as  well  as  mean  canopy  humidity  and  temperature,  A representa- 
tive daily  sample  output  comparing  leaf  temperatures  with  air  tempera- 
ture is  shown  In  Fig.  4.6  for  Oct  IZ  which  was  a relatively  clear  day 
except  for  some  clouds  at  the  beginning  and  end  of  the  day.  Leaf 
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I’ 


Fig.  4.6 


Comparison  or  daily  course  oF  sinulated  sun-1H 
and  snaded  leaf  ensemble  tenperacures  wUn  tbe 
air  CenperaCpre  on  Oct  13.  The  leaf  temperature 


134 


temperatures  Increased  adpve  air  temperature  witn  an  Increase  in  COj 
levels,  nils  response  was  the  direct  result  of  COj  effect  an  stomatal 
beKavior  causin9  lower  transpiration  and  STmuUaneous  rise  in  the  leaf 
temperatures.  The  measured  leaf  temperatures  with  the  1.1-1600  steady 
state  porometer  were  usually  lower  than  the  simulated  temperatures. 
This  error  was  attributed  to  the  physical  operation  of  the  porometer. 
First,  during  leaf  temperature  measurements  a leaf  is  enclosed 
between  the  metal  rings  which  dampens  the  temperature  and  additional- 
ly the  fan  In  the  cuyette  blows  air  directly  on  to  the  leaf  causing 
euaporative  cooling  that  lowers  leaf  temperatures.  J.  H.  Sennett, 
(Personal  Communication)  conpared  Infrared  thermometer  measured  leaf 
temperatures  with  the  steady  state  porometer  measured  values  and  found 
that  the  values  were  lower  by  as  much  as  5-7  “C.  The  simulated  tempe- 
ratures response  to  COg  were  similar  to  the  single  leaf  temperatures 
reported  by  Valle  et  al.  <I9SS). 
long  Term  Mater  Use 

Simulation  of  dally  cycles  end  their  agreements  witn  measurements 
gave  confidence  in  the  model  structure  and  the  essumptions  made  In  its 
development.  Tne  model  was  next  used  to  simulate  ET  on  the  remaining 
days  to  evaluate  the  ability  of  the  model  to  simulate  daily  water 
use.  The  summary  results  comparing  simulated  and  measured  ET  are 
presented  in  Tables  4.3-4. 5.  In  summary  the  low  COg  treatment  had 
higher  water  loss  compared  to  Its  counterparts.  The  climatic  condi- 
tions were  so  variable  during  the  study  period  that  clear  trends  in 

have  stabilized  when  the  canopies  reached  an  LAI  around  3-3.5.  The 
correlational  diagrams  (Fig.  4.7)  compering  simulation  and  measure- 
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OAPL  LAI  oAI 
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Measured  ET.  simulated  Et.  7 and  E and  errors  Detween 
the  neasurement  and  simulation  for  the  990  COj  treatment 
on  various  days  after  plantlns  (OAPL)  and  days  after 
Irrigation  (DM).  The  eateoslon  -C  following  the  OAl 
column  Indicates  Chat  Che  soli  surface  was  covered, 
whereas  as  the  extension  «S  Indicates  suDDlenenta! 
Irrigation, 


MEASWED  SIHUUTED  ERROR  SIHUUTEO  SIMULATED 
OAPL  Ul  ET.rnm  ET,mm  m T,mm  E,mm  OAI 
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Co-relatlonal  illigrani  os  slnulatefl  and  measured 
evapotransol ration  for  the  330,  S60,  and  990  CO7 
treatments  conpared  with  1:1  line.  Days  used  In 
the  calibration  are  indicated  with  large  soilld 
dots.  The  difference  between  simulated  and 
measured  values  were  insignificant  at  95S  level 
of  confidence 


revealed 


that  there  nere  ho  signiffcant  differe- 
ecea  between  the  sfnelation  and  neaaurenents  at  951  confidence 
teyel.  The  n«»1iniin  errors  during  the  entire  simulation  were  12t.  23t, 
and  18X  while  the  average  absolute  errors  were  6X.  3.K  and  4.71  for 
the  330,  660,  and  990  treatments  respectively. 

Soybean  Water  Use  With  Increased  CO?  Levels 

The  development  of  the  computer  simulation  has  made  possible  rapid 
evaluation  management  alternatives.  In  this  section,  the  model  was 
used  to  forecast  water  use  of  a field  grown  soybean  under  the  present 
CO^  level  (330  ppm)  and  future  anticipated  CO^  level  of  990  ppm.  In 
the  esoerlmental  study  it  was  found  that  the  stomatal  resistance  was 
increased  abd  leaf  area  development  was  rapid  In  the  high  COn  treat- 
ments compared  over  the  low  COj  level.  It  Is  anticipated  that  the 
rise  In  envi rcnnental  CO^  will  lead  to  global  Increase  In  average  a1r 
temperature  by  1 to  5 ”C.  The  water  use  under  the  influence  of  each 
of  these  three  factors  was  examined. 

Assumptions 

Soybean  crop  was  assumed  to  Oe  groan  in  the  constant  weather 
conditions  shown  In  Fig.  4.8.  The  leaf  area  Index  data  fdr  the  330 
and  990  treatments  presented  In  Fig.  4.2  was  assumed.  The  soybeans  in 
this  study  were  harvested  115  DAPL.  therefore  between  66  and  115  DAPL. 
LAI  was  decreased  linearly  with  time.  The  assumed  LAI  during  66  to 
115  day  Interval  was  similar  to  recorded  field  measurements  (E, 
Bishal,  personal  communication).  At  full  canopy  the  height  of  the 
crop  was  0.5  m,  light  extinction  coefficient  0.5,  velocity  attenuation 
coefficient  3,  and  uniform  wind  speed  of  1.5  m/sec.  Refer  to  Chapter 
111  for  detail  model  inputs  and  soil  characteristics. 


Fig.  4.3  Tde  iiuml  uiriadon  in  air  tenperatura,  dew 
point  tenparatijra,  net  radiation,  and  P*F*FD  osad 
on  all  days  in  the  simulation  oF  induced 
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StoiMtal  Resistance 

Tne  CO2  concertratlons  affects  leaf  stomatal  resistance.  If  all 
ccndltlons  such  as  light,  water  anfl  temperature  are  at  optimum  values, 
the  minimum  stomatal  resistance  occurs  when  the  external  COj  level  Is 
about  egual  to  the  COg  compensation  point,  the  point  where  the  photo- 
synthetic  COj  uptake  Just  balances  COp  released  by  respiration  (Ting, 
1982).  Uhereas  low  CQp  tends  to  cause  stomatal  opening,  high  COp  will 
tend  to  close  stomata.  The  stomatal  resistance  data  presented  1n 
Table  4,1  were  used  to  simulate  the  relative  contribution  of  this  COj 
induced  change  in  crop  water  use.  The  results  of  the  simulation  runs 
are  sumnariaed  in  Table  4.S,  In  the  following  discussion  the  abbrevi- 
ation SR  was  used  to  indicate  stomatal  resistance,  and  the  two  COj 
levels  and  corresponding  LAI  end  S«  were  referred  as  low  (330 
treatment)  and  high  (990  treatment)  respectively.  The  simulation 
result  snowed  that,  with  LAI  and  SR  that  occurred  in  330  ppm  COj.  over 
the  115  days  of  simulation,  the  crop  transpired  348.5  mm  of  water  (80* 
of  total  water  consumption)  compared  to  289.3  nm  ( 75  1 of  total  water 
consumption)  with  the  same  LAI  but  with  SR  for  the  990  COj  plants. 
There  was  a net  reduction  In  water  use.  Hore  realistically,  however, 
increases  in  COp  will  cause  higher  LAI  due  to  higher  photosynthesis. 

The  true  comparison  must  be  made  between  low  LAI  and  low  SR  with 
high  LAI  and  hign  SR,  The  total  water  loss  from  soybeans  grown  in 
high  COp  environmenl  was  reduced  from  436.5  nm  to  394.1  mi  ( a 19* 
reduction)  whereas  transpiration  loss  was  reduced  by  13*.  The 
results  of  simulation  compare  well  with  water  saving  reported  by 
Jones  et  al.  (1985).  They  reported  that  a 10*  savings  in  transpirat- 
ion when  soybeans  were  grown  In  660  COp  concentration  compared  with 
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Effect  of  scenatal  resistance  and  leaf  area  ctevelopnenC 
of  weekly  irrigated  soybeans  in  330  and  990  opm  COn 
concentrations  on  seasonal  water  use.  ‘ 


330  oom  LAI 

SR-330  SR-990 

( LOW)  (HIGH) 
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33D  CO^  grown  soybeans  and  Cbaa  the  total  water  loss  was  346  nn  conpa* 
red  to  400  rm  over  the  lOS  day  growing  season.  The  cunulative  trans- 
piration loss  due  to  changes  in  stomata!  resistance  are  compared  In 
Fig.  4.9.  Increased  stomatal  resistance  led  to  1.5  to  1.75  °C  rise 
both  Che  sunlit  and  shaded  leaf  temperatures  (Flg.  4.10  and  Fig. 
4.11). 

Leaf  Area  Development 

The  effect  of  leaf  area  development  by  Itself  on  transpiration  at 
low  and  nigh  SR  Is  compared  1n  Fig.  4.12.  In  this  case  1t  can  be  seen 
that  the  water  use  by  transpiration  was  not  very  sensitive  to  leaf 
area  development,  probably  because  leaf  areas  were  similar  under  both 
CO^  environments. 

Air  Temperature 

To  simulate  the  effect  of  Increased  temperature  on  water  use, 
each  hourly  temperature  was  increased  by  5 *C  and  simulations  were 
performed  for  high  LAI  under  both  the  low  and  high  The  water  use 
components  are  sunnarlted  in  Table  4.7.  Incnease  in  ain  temperature 
caused  Increased  vappr  pressure  deficit  and  seasonal  transpiration 
[Fig.  4.13)  and  caused  about  3 'C  ( Figs.  4.14  -4.15)  highen  leaf 
temperature.  In  the  previous  simulations  with  increased  stomatal 
resistance,  the  effect  of  increase  in  temperature  on  stomatal  resista- 
nce was  ignored.  Therefore  a true  test  to  see  the  effect  of  COj 
induced  environmental  changes  on  water  use,  would  be  to  compare  the 
effect  of  stcmatal  resistance,  leaf  area  and  temperature  together. 
This  comparison  can  be  derived  fron  the  first  column  of  Table  4.6  and 
the  last  columh  of  Table  4,7.  According  to  the  model,  the  total  water 
loss  will  be  Increased  to  520.4  nil  compared  to  436.5  by  increasing 
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BAYS  AFTER  p2ntd«'** 


ftg.  4.10  CoopiH5on  of  sirajiated  a»ei-age  daily  sufi-lU 
laaf  BosemDle  temperature  at  low  {330  CO,)  and 
nign  (990  CO,)  LAl  under  the  influence  of 
measured  stonatal  resistance  of  the  low  and  nigh 
Co,  treatments 
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Table  4.7  Comparison  seasonal  water  use  of  soybeans  with  stomata! 

resistance  corresponding  to  330  and  990  ppm  CO; 
concentrations  under  present  and  future  global  air 
temperature.  LAI  of  990  CO;  treatment  was  used. 


SR  -990  COp 
Tain  Tain 


ET  446.0  583.6 


394.1 
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Comparison  of  sinulaCod  average  daily  spaded 
leaf  ensemble  temperature  under  present  and 
future  { S *C  higher)  air  temperature  simulated 
using  leaf  area  and  stomata  resistance 
corresponding  to  990  CO^  treatments 


atnospN«rlc  CO2  levels  CO  990  conpared  to  present  330  oofn.  Tnis  anal- 
ysis used  the  extreme  value  predicted  for  a Cenperature  change  {5 
and  no  change  In  atmospheric  vaoor  pressure.  However.  It  is  11kei> 
that  the  atmospheric  vapor  pressure  would  be  Increased  with  an 
Increase  in  temperature.  These  results  show  a wide  range  in  predicted 
water  use,  varying  from  a decrease  By  10%  If  atmospheric  temperature 
does  not  change,  to  an  Increase  by  19  % if  atmospheric  temperature 
changes  by  5 and  no  change  takes  place  in  atmospheric  vapor 
pressure. 

Summary  And  Conclusions 

The  principle  goals  of  this  paper  nere  to  evaluate  Che  evapotra- 
nspiration  model  of  Chapter  III  witn  data  from  soybean  crop  grown  in 
controlled  environment  chamOers  under  three  CG2  concentrations  of  330, 
660  and  990  ppm  and  constant  dry  bulb  and  dew  point  temperatures,  and 
to  apply  the  model  to  study  the  possible  effects  of  changes  in 
atmospheric  COg  concentration  on  water  use  of  field  grown  soybeans. 

The  model  simulated  results  agreed  well  witn  the  measurements  at 
all  times  of  the  day  evcept  in  the  mornings.  The  probable  reasons  for 
this  were  attributed  to  the  natural  lag  in  the  measurement  system, 
cold  outside  temperature  relative  to  the  controlled  chamber  temperatu- 
re causing  condensation  on  the  side  chamber  walls,  and  the  validity  of 
the  steady  state  assumptions  during  rapidly  changing  morning  envi- 
ronment, heat  sinks  In  chamber  walls,  and  to  the  use  of  Instantaneous 
rather  than  average  net  radiation  values  over  tne  time  step  used  in 
the  simulation.  The  average  dally  errors  over  the  simulated  period 


respectively. 


Tne  Sunlit  and  snaded  leaf  tenperatures  Increased  with  CO2 
levels.  The  slnjlated  values  agreed  well  with  tne  single  leaf  temper- 
atures responses  to  COg  reported  hy  Valle  et  a1 . {19B5)  Put  Mere  loner 
tpan  tne  steady  state  porometer  values. 

The  model  provided  an  excellent  tool  to  analyze  changes  In  the 
microclimate  and  water  use  by  a developing  soybean  crop.  Altnougb  it 
Is  difficult  to  extend  tne  model  findings  to  field  situations  due  to 
Che  constant  environmental  conditions  of  the  chamber,  it  is  certainly 
a valuable  tool  to  analyze  climate  and  crop  interactions.  Simulation 
studies  were  conducted  to  study  the  effect  of  CO^  induced  changes  in 
stomatal  resistance,  leaf  area  development,  and  gloDal  air 
temperature.  For  a fixed  leaf  area,  Increased  stomatal  resistance 
reduced  water  use  and  increased  leaf  temperatures.  Transpiration 
water  use  under  low  and  high  COg  treatments  differed  by  only  13  t and 
agreed  well  with  experimental  findings  of  Jones  et  a1.  (19851.  The 
individual  effect  of  a 5 °C  In  air  temperature  on  leaf  temperature  and 
water  use  was  dramatic.  The  leaf  temperatures  may  be  higher  by  as 
much  as  3 “C,  and  water  use  was  predicted  to  Increase.  When  the  proj- 
ected effects  of  COg  on  stomatal  resistance,  leaf  area  develomenc,  and 
temperature  were  combined,  model  predictions  of  water  use  were  Increa- 
sed by  as  much  as  19t  due  to  COg  enriclment  plus  5 X temperature 
increase. 

The  study  provided  an  excellent  first  test  and  validated  many  of 
the  model  predictions.  The  overall  assessment  of  the  model  was  satis- 
factory and  further  test  of  the  model  under  the  field  conditions  is 
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CHAPTER  V 

ESTIHATING  CROP  COEFFICIENTS  FOR  SOYBEAN  UNDER  OIFFERENT 
IRRI5ATION  HANASEHENT  PRACTICES 

Introduction 

In  many  parts  of  ttie  worlrt  «nere  »ater  holding  capacity  of  the 
soil  Is  low  and  procIpUatlon  falls  short  of  water  requi rements  or  is 
irregular  during  the  growing  season,  irrigation  is  practiced  to  avoid 
drought  and  loss  of  yields.  A critical  problem  in  Irrigation  is  to 
determine  Just  when,  how  and  how  much  water  to  apply  to  mayimiie  the 
benefit  to  cost  ratio.  There  are  at  least  three  approaches  to  deter- 
mine the  optimal  timing  and  quantity  of  irrigation.  One  way  Is  to 
observe  crop  Indicators  such  as  change  of  color  or  leaf  angle  in  an 
attempts  to  detect  early  signs  of  stress.  Secondly,  the  direct  mea- 

available  is  adequate  to  permit  water  extraction  by  roots  at  potential 
rates.  Finally  the  weather  can  be  monitored  to  determine  evapotrans- 
piration  (£T)  demand  of  a crop  and  dompare  It  with  the  supply  of  water 
from  rainfall  and  Irrigation  events. 

Irrigation  scheduling  using  ET  confuted  from  weather  record  is 
appealing  because  the  approach  is  relatively  simple  compared  to  on- 
site ET  measurements.  In  one  approach  to  estimating  ET,  the  whole 
realm  of  crop  responses  to  the  soil-  plant-climate  and  management 
factors  Is  simply  taken  care  by  " crop  coefficients  " for  various 
crops  at  different  stages  of  growth.  In  this  approach  ET  of  a 
disease-free  crop  grown  in  large  field  under  optimum  soil-water  and 
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fortuity  conditions  and  achisvin9  full  production  potential  under  a 
given  environment  Is  determined  using  reference  crop 
evapotransplratlon  (ETp)  and  a coefficient  (X,.)  to  account  for  stage 
of  growth  of  the  crop.  Accordingly, 

ET  ■ X(.  ETr  (5.1) 

where  Xj  Is  the  crop  coefficient.  Crop  coefficients  for  many  crops 
nave  oeen  developed  In  field  eiperlments  and  are  siininarued  by 
Ooorenbos  and  Pruitt  (1977).  Crop  coefficients  are  sensitive  to  soli 
evaporation  especially  when  canopies  are  small  with  Incomplete  ground 
cover  and  If  Irrigation  or  rainfall  occurs.  The  data  presented  by 
Ooorenbos  and  Prult  (1977)  snowed  as  oucn  as  70  - 80  I variation  In 
crop  coefficient  during  the  early  phase  of  crop  development  depending 
on  whether  the  soil  surface  was  dry  or  wet.  As  the  crop  cover 
develops,  crop  coefficients  become  less  sensitive  to  variation  1n  the 
soil  surface  conditions.  Therefore,  attempts  have  been  made  to  syn- 
thesise the  overall  crop  coefficient  Into  plant  and  soil  coefficients 
(Mright,  I9B2;  Phene  et  a1..  1985).  Accordingly.  Eqn  (5.1)  can  be 
rewritten  to  Include  soil  specific  coefficient  and  plant  transpira- 
tion coefficient  X.  as: 

Xj  . ET/ETr  ■ E/ET,.  * T/ET..  ■ Xj  * K„  (5.2) 

It  Is  argued  that  with  differences  In  weather  accounted  for  by  ET^. 
the  transpiration  coetficient  Xp  would  be  universally  valid  for  a 
specific  well-watered  crop.  Or  the  other  hand,  Kj  would  vary  locally 
depending  upon  Irrigation  frequency,  method  of  Irrigation,  soil  chara- 
cteristics and  degree  of  soil  cover  (Hright.  1982;  Phene  et  al..  1985; 
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It  Is  common  practice  to  use  mean  climatic  data  for  determining 
ITf.,  However  due  to  changes  In  weather,  ET,.  will  vary  from  year  to 
year,  for  each  period  within  a year  and  also  due  to  changes  In  surface 
conditions.  Smalstrla  et  al.  (19041  have  presented  probability  dis- 
tribution of  monthly  ET,.  values  for  nine  Florida  locations  using  2S 
years  of  weather  data.  According  to  Ooorenbos  and  Pruitt  (1977)  the 
monthly  values  could  vary  by  50%  or  more  from  one  year  to  the  neat. 
The  daily  values  can  also  vary  drastically,  with  low  values  on  the 
days  that  are  rainy,  cloudy,  humid  and  calm  with  high  values  on  dry, 
sunny,  and  windy  days.  There  are  no  clear  guidelines  to  adjust  crop 
coefficients  due  to  weather  changes.  Therefore  the  accuracy  of  crop 
coefficients  Is  questionable  unless  both  the  ETr  and  crop  coefficients 
are  obtained  under  the  same  climate. 

ET  Models 

Determination  of  crop  coefficients  with  a direct  method  Is  a time 
consuming  and  expensive  procedure.  To  determine  crop  coefficients, 
generally  a minimum  of  two  sensitive  weighing  lysimeters  are  required 
In  such  a study  so  that  both  the  ET  and  ET^  are  measured  simultane- 
ously, A1  ternatively  crop  coefficients  could  also  be  determined  using 
a computer  simulatipn  model  that  computes  evaporation  and  transpira- 
tion separately  In  response  to  soil,  climate  and  cnop  variables, 
ilodels  have  been  developed  that  compute  evaporation  and  transpiration 
separately  (Ritchie,  1972;  Tanner  and  Jury,  1975;  Chapter  III).  The 
models  of  Ritchie  (1972)  and  Tanner  and  Jury  (1975)  have  related  tra- 
nspiration to  the  leaf  area  Index  (UI)  and  potential  transpiration 
using  empirical  functions.  Similarly,  soil  evaporation  was  related  Co 
Che  potential  evaporation  when  soil  surface  was  wet  followed  by  an 
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eitip1r1cdl  time  function  es  soil  dries.  On  trie  other  hand  Che  model  of 
Chapter  III  treated  the  soil-plant-atmosphere  as  an  energy  exchange 
system.  This  approach  cakes  into  account  the  discrlpucion  of  energy 
sources  in  the  canopy  and  soil  as  well  as  physical  characteristics  of 
soil,  plant  and  atmosphere,  and  their  interactions  to  compute  soil 
evaporation,  transpl ration,  canopy  microclimate  (temperatures  and 
hunldny).  spil  temperatures,  and  soil  heat  flux.  The  model  is 
diurnal  1h  nature  and  requires  half  hourly  values  of  air  temperature, 
dew  point  temperature,  net  radiation,  and  wind  speed  as  inputs.  It 
also  needs  LAI.  crop  light  extinction  cnaracteristics,  and  LAI  and 
height  at  full  canopy 

The  objectives  of  this  paper  were  to  use  the  model  of  Chapter  II 

1.  Determine  the  effect  of  variation  in  climatic  variables 
such  as  temperature,  vapor  pressure  deficit  (VPD), 
radiation  and  wind  speed  on  crop  coefficients  K^.  kp,  and 
<5  under  three  irrigation  frequencies;  irrigated  daily, 
once  every  three  days  and  weekly. 

2.  Determine  the  effect  of  planting  dates,  crop  development 
rates,  and  length  of  the  growing  season  on  water  use  and 
crop  coefficients  for  a soybean  crop  grown  in  sandy  soil 
and  using  actual  weather  data  in  humid  Florida  climate. 

3.  Discuss  possible  applications  and  linkages  of  the  mode! 
with  crop  growth  models. 

Waterials  And  Wethods 

*010  simulatidn  study  was  divided  into  two  parts  to  fulfill  the 
first  two  objectives.  During  the  first  part,  a set  of  input  values 
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for  soil,  crop  aod  meteorological  parameters  were  selected  to 
represent  a soybean  crop  gro»n  on  sandy  soil  1n  Florida.  The  soil  was 
represented  by  eight  layers  ( of  10,10.10,30,30,40,40,40  on  thickness) 
with  each  layer  having  field  capacity  water  content  of  101  and  lower 
extraction  limit  of  it  on  the  volumetric  basis.  Tne  thermal  conducti- 
vity and  heat  capacity  were  0.574  J/m-sec-’C  and  1.126E6  J/m^-“C  (at 
0»  moisture  content)  and  4.522  J/in-sec-*C  and  2.663E6  J/n^-’C  (at  13.1 
I moisture  content)  based  on  data  reported  by  Herva  (1975)  for 
granitic  sand.  The  thermal  properties  at  Intermediate  water  contents 
were  linearly  interpolated.  Soil  temperatures  were  simulated  using  a 
one  dimensional  neat  transfer  model  with  uniform  initial  temperature 
of  26  ’C.  after  each  Irrigation  new  Initial  temperatures  were  computed 
as  the  weighted  average  temperature  of  soli  and  water  by  assiming  that 
added  water  had  average  air  temperature. 

Soil  evaporation  occurred  from  the  top  30  cm  of  soil  while  tra- 
nspiration loss  occurred  throughout  the  soil  profile.  When  the  canopy 
was  small  (LAI  < 0,5)  all  of  Che  trahspiration  loss  occurned  from  the 
top  30  cm,  however  afterward  this  proportion  was  reduced  to  51  for  the 
remaining  growth  period.  The  percent  extraction  values  of  transpirat- 
ion were  selected  arbitrarily  to  Include  the  effect  of  root  growth, 
(n  Chapter  III  It  was  found  that  seasonal  water  use  of  a well  watered 
crop  was  insensitive  to  this  ratio.  The  leaf  area  Index  (LAI)  with 
days  after  planting  (OAPL)  was  confuted  using  the  S0V6R0  model 
(Wilkerson  et  al,,  1984),  For  the  verlety  CobO  planted  on  June  26  and 
August  15,  and  the  variety  Hilllans  planted  on  March  16  In  Gainesville 
Florida  with  a planting  density  of  36  plants/m^,  and  row  and  plant 
spacing  of  0.762  and  0.0365  m respectively  (Fig.  5.1).  The  crop  light 
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eitinctisn  coefr!c<9tit  of  0.5  [spherlcil  leaf  angU  OUtMOutlon)  and 
stoiDatal  conductance  as  a function  of  PPFD.  front  Chapter  1)1,  1$  given 
by: 


"bere  C,,,,  «as  nailnum  conductance  (0.018  m/sec)  and  s {38.14*10"*  ) 
was  Initial  slope  of  C vs  PPFD  ( E/u^/see). 

The  climatic  data  typical  for  tne  June  planting  date  was  obtained 
from  Jones  et  a1.  (198Aa)  and  is  listed  In  Table  5.1,  The  hourly  air 
temperatures  were  ccmputed  using  the  diurnal  algorithm  in  Che  SOYGPO 
model.  The  dew  point  temperature  and  wind  speed  were  assumed 
constant.  Dally  net  radiation  was  computed  using  the  procedures  of 
Jones  et  al.  (1984)  and  It  was  distributed  dlurnally  between  sunrise 
and  sunset  using  a sine  function.  The  PPFD  between  sunrise  and  sunset 
was  obtained  by  muUlplylng  net  radiation  with  0,10,  a locally  determ- 
ined empirical  constant. 

To  assess  the  Impact  of  weather  on  water  use  and  crop  coeffic- 
ients the  sensitivity  tests  were  carried  out  by  varying  standard  air 
temperature,  relative  humidity,  radiation  and  wind  speed  as  per  the 
summary  In  Table  5,1,  The  temperatures  were  varied  such  that  the 
maximum  and  minimum  daily  temperature  difference  was  always  10  *C.  In 
a natural  system  cnanges  in  temperature  would  bring  changes  in  the 
VPD.  However,  during  this  phase  of  the  simulation  simulation  eiperi- 
ments,  the  reletive  humidity  was  adjusted  to  maintain  a constant  VPD 
of  7.5  mb  (millibar).  The  range  of  temperature  variation  was  typical 
of  moderate  to  warm  agro-eco-system.  The  relative  humidity  levels 
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TABLE  5,1  Wesclter  rlats  used  for  the  sensftWHy  analysts  and  the 
range  of  variations  considered. 


Parameter  Symbol  Units  Standard  Simulated 

value  values 


Haiimun  temperature  Tnax  °C 

Minimum  temperature  Tmin  'C 

Radiation  RAO  cai/cn* 

Relative  humidity  RH  t 

Hind  speed  MS  Wday 
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Mere  characceristlcs  of  Ounid,  suD-ruenid  and  arid  regiorvs.  For  this 
region  the  clear  sty  dally  radiation  averages  as  high  as  744  ca1/cn^  ( 
Jones  et  a1.,  19B4a].  Therefore  to  study  the  effect  of  variability  in 
the  cloud  cover,  situations  with  clear  sty  with  maxiniuni  possible  radi- 
ation and  60t  of  average  seasonal  radiation  sere  considered.  Wind 
speed  ranged  from  86,4  to  518.4  Kn/day, 

Water  holding  capacities  of  the  sandy  soils  are  low,  Uierefore, 
snail  frequent  applications  early  in  Che  growing  season,  and  greater, 
less  frequent  applications  later  are  suggested  to  inprove  irrigation 
scheduling  (Clark,  ig82a].  Therefore  daily,  once  every  three  days, 
and  weekly  irrigation  frequencies  were  considered.  Since  the  pritnary 
goal  of  this  paper  was  to  detemine  crop  coefficients,  soybeans  were 
assumed  to  be  well  watered  at  all  tines.  Irrigation  was  applied  at 
the  beginning  of  each  day,  and  no  irrigation  was  applied  following 
physiological  maturity  of  the  crop. 

To  fulfill  the  second  objective  of  this  study,  two  additional 
experimental ly  feasible  planting  dates,  spring  (Narch  15.  variety 
Wililans)  and  fall  (August  IS,  variety  Cobh)  were  considered  in  addi- 
tion to  the  June  planting  date.  Monthly  historical  weather  data  repo- 
rted by  Jones  et  al.  (1984a)  was  used  for  the  simulation  (Table 
5.2).  Daily  values  were  obtained  by  linear  interpolation  of  each 
variable,  assuming  the  tabular  values  occurred  on  the  middle  of  the 
month.  Soil  and  crop  characteristics,  irrigation  frequencies  as  well 
as  all  other  procedures  were  identical.  The  ETr  was  computed  using 
the  Penman  method  (Jones  et  al..  1984a)  and  crop  coefficients  were 
computed  using  Eqn.  (5.2). 
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Averdge  monthly  cHmatologlcal  Oata  used  to  simulate 
effect  of  planting  dates  on  water  use  and  crop 
coefficients  of  soyOeans  grown  In  Florida.  The  data  is 
for  the  Belle  Glade.  Florida,  weather  station  (Jones  et 
al..  19S4). 
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— [ 0.263  ( e,  - ) (0.5  ♦ 0.006  U ) (5.3) 

nhere  ■ ^7*  (0.56  - 0.08  ) (1.42  - 0.42  ) 

and  T » average  dally  air  temperature  in  •*  (*C  ♦ 273) 

ej  = vapor  pressure  ■ ♦ ®niln'^^  ■ 

e„  " vapor  pressure  at  dew  point  temperature,  mb 
®max*  "'Saltnum  vapor  pressure  of  air  during  a day,  mb 

e|„jn*  minimum  vapor  pressure  of  air  during  a day,  mb 

Rjg  • total  daily  cloudless  sty  radiation 
6 • slope  of  saturated  vapor  pressure  curve  of  air,  mO/"C 
7^  * psychrometric  constant.  0.66  nb/*C 
> a latent  beat  of  vaporization  of  water  58  caT/tn^ 

" albedo  of  surface  for  Rj 
^ ■ Stefan-BoTtzmann  constant 
Rj  = total  Incoming  radiation  cal/cm^.day 
U ■ wind  speed  at  height  of  2 m,  Km/day 
RESULTS 

Simulations  were  performed  on  the  PRIME  550  ccmputer  in  the  agri- 
cultural Engineering  Bepartment  at  University  of  Florida.  The  length 
of  growing  season  varied  with  planting  dates;  117  days  for  the  March 
15  planting,  134  days  for  the  June  26  planting  and  110  days  for  the 
august  planting.  Therefore,  to  compare  the  simulated  results  the 
duration  after  planting  was  normalized  to  give  time  elapsed  as  a 


percenta9^  of  growing  season  unless  mentioned  otherwise. 

Irrioation  Freouency 

Evepotranspiretlon  (ET).  soil  evaooration  (E]  and  transpiration 
(T)  along  with  crop  coefficients  tor  the  standard  constant  weather 
condition  and  Che  June  slanting  dace  are  shown  in  Fig.  E 2.  With  the 
daily  irrigated  treatment,  tne  surface  was  re-wetted  each  morning  and 
allowed  soil  evaponation  to  occur  at  potential  rate.  Soil  evaporation 
was  less  than  ETr  and  caused  increase  in  soil  surface  temperature.  As 
the  soil  surface  temperature  Increased,  the  saturated  vapor  pressure 
was  higher  causing  nwre  soil  evaporation  until  soybeans  einerged  on  the 
ninth  day.  Similar  phenomena  was  observed  during  the  three  and  seven 
days  Irrigation  frequencies.  The  seven  day  frequency  allowed  soil 
surface  to  warm  even  more  than  the  other  two  treatments  and  therefore 
resulted  in  relatively  higher  evaporation  following  an  irrigation 
event.  As  the  canopy  developed,  the  energy  availaole  at  the  soil 
surface  became  limited  and  soil  evaporation  was  reduced.  At  high  LAI 
during  Sn-119  OAPL,  elmost  a constant  portion  of  the  Incoming  radia- 
tion was  captured  by  the  canopy,  allowing  only  a small,  constant 
amount  to  reach  the  soil  surface  as  illustrated  Oy  the  constant  evapo- 
ration during  this  period.  Following  crop  maturity  (119  OAPL)  irriga- 
tion was  stopped  and  soil  evaporation  increased  as  LAI  dropped  and 
more  radiation  reached  the  soil  surface.  Soon  this  trend  reversed  as 
the  soil  became  dry. 

As  tbe  irrigation  interval  was  Increased  to  3 and  7 days,  the 
soil  surface  went  through  cycles  of  wetting  and  drying  as  illustrated 
by  tne  varying  amplitudes  of  E and  ET.  The  soil  evaporation  occurred 
at  the  potential  rate  following  irrigation,  but  decreased  rapidly  as 
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tne  surfice  dried.  As  wltn  We  dally  irrigation,  the  canopy  develop- 
ment reduced  available  energy  at  the  soil  surface,  the  amplitudes  were 
damped  out,  hut  mere  not  completely  eliminated  even  at  maximum  LAI, 
It  mas  also  Interesting  to  see  the  “dips"  in  the  T rates  following  an 
Irrigation  rfilch  Increased  soil  evaporation  and  himidity  surrounding 
the  canopy  or  reduced  the  vapor  pressure  gradient  betmeen  the  leaves 
and  canopy.  The  Interaction  of  £ and  T seemed  to  nave  mutually  compe- 
nsating effects.  This  relationship  could  only  be  observed  and  studied 
using  a physical  model  of  this  nature  mnich  is  one  of  the  major  advan- 
tages of  this  model  over  Ritchie's  [1972)  or  Tanner  and  Jury's  (1975) 
models.  The  feedback  of  E an  T became  less  significant  at  lomer  mag- 
nitudes  of  £.  The  dally  variation  in  K_,  K_,  and  K_  mere  similar  to 
£T,  E and  T as  the  El,,  was  not  affected  by  the  Irrigation  cycles. 

Unlike  ET.  by  definition  ET^  mas  not  influenced  by  the  Irrigation 
management  practice.  Therefore  the  hehavior  of  crop  coefficients  mere 
similar  to  the  corresponding  water  loss  components.  Excluding  the 
first  nine  days  (bare  soil)  and  the  period  following  crop  maturity,  ET 
exceeded  ET,.  ( Kg>l  ).  The  soil  evaporation  was  lower  than  ET,., 
except  during  the  second  and  third  drying  cycles  with  the  three  and 
seven  day  frequencies,  therefore  Kj  was  less  than  unity.  Similarly, 
transpiration  was  lower  than  ET,.  until  LAI  of  1.63  ( Up  < 1)  and  then 
exceeded  ETr  ( Xp  > 1)  up  to  the  beginning  of  leaf  senescence.  ET,. 
decreased  gradually  with  canopy  development  due  to  the  effect  of  dif- 
ference between  soil  and  crop  albedos. 

Tne  seasonal  mater  use  (Table  5.3)  mas  affected  by  the  irrigation 
method,  mainly  due  to  changes  in  soil  evaporation  during  early  develo- 
pment- During  the  first  fifty  days  cumulative  evaporations  from  the 
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TABLE  5.3 

Cumulative  reference  crop  evapotranspi ration  (cETr) 
computed  uting  Che  Penman  Method,  simulated  soybean 
evapotranspiradon  (cET),  soil  evaporation  (c£)  and 
transpiration  (cT),  using  the  standard  values  of  air 
temperature,  radiation.  VPD.  and  wind  speed  (Table  1), 
under  daily,  once  every  three  days  and  weekly  irrigation 

Variable 

1 Day  cycle  3 Day  cycle  7 day  cycle 

iiil 

588.3  589.8  589.8 

627.4  608.5  585.2 
170.8  147.5  113.9 
455.6  461.0  471.3 

dsy  frequencies 


Here  132,  100  and  71  nm  respect- 
ively compared  to  30.0,  30.3  and  43  nn  during  rest  of  tde  season. 
Therefore  the  difference  In  and  iC,  occurred  mainly  during  the  inco- 
mplete canopy  pfiase,  while  for  tne  major  part  of  the  peak  growihg 
season,  these  coefficients  were  practically  the  same  among  all  irriga- 
tion treatments.  On  tne  other  hand  transpiration  was  little  changed 
due  to  differences  in  the  irrigation  treatments  and  Kp  remained  pract- 
ically unchanged. 

In  Che  following  discussion  the  effects  of  climate  will  he 
analysed  to  determine  what  changes  to  aspect  in  water  use  pattern  and 
crop  coefficients.  To  improve  the  clarity  of  graphical  presentations 
all  daily  water  loss  components  and  crop  coefficients  were  averaged 
over  the  interval  between  irrigations. 

Temperature 

The  maximum  and  minimum  daily  temperatures  were  changed  by  a 5 "C 
over  the  standard  values.  Examining  Eqn.  (5.3),  the  availaole  net 
radiation  increased  with  increase  in  temperature.  To  maintain  the  VPD 
constant  at  7.5  mh.  the  dew  point  temperature  and  consequently  vapor 
pressure  at  dew  point  temperature  increased  and  reduced  long  wave 
thermal  radiation.  Therefore  both  ET  and  ETp  increased  with 
temperature  ( Table  5.4).  The  change  in  temperature  also  influenced 
the  sensible  heat  exchange  between  the  canopy  and  the  surrounding 
air.  Uith  a daily  irrigation  cycle,  lower  air  temperature  reduced 
leaf  temperature  and  transpiration  by  14*  while  an  increase  by  tne 
same  amount  caused  a 17*  higher  transpl ration  (Table  5.4).  The 
unequal  change  in  transpiration  was  attributed  to  the  non-] ineartty 
between  temperature  and  saturated  vapor 


pressure  relationship. 
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Slnllar  behavior  with  3 and  7 day  frequencies  was  observed.  At  high 
LAI.  Increasing  teriperaCure  caused  higher  leaf  Cenperatures  and  higher 
vapor  pressure  gradient  between  leaves  and  canopy  air  causing  more 
transpiration.  This  was  illustrated  by  the  shift  between  T and  in 
Fig  5.3a,  On  the  other  hand,  higher  transpiration  reduced  soil  evapo- 

The  changes  in  Kp  and  K,  associated  with  temperature  change 
are  shown  in  Fig  5.3b.  As  net  radiation  Increased,  bare  sol!  evapora* 
tion  increased  over  a short  tine  but  soon  became  limited  by  available 
water.  At  low  temperature  and  before  plant  emergence  it  exceeded  ET., 
but  with  further  rise  in  temperature  and  net  radiation  it  tailed 
behind  ET^.  Therefore  Kj  decreased  with  temperature.  On  the 
otherhahd,  transpiration  of  the  well-watered  soybean  Increased  propor- 
tionately and  1C  were  similar  among  all  temperature  simulations. 
Later,  when  the  canopy  reached  its  peal:  radiation  absorption  characte- 
ristics, vapor  pressure  deficits  induced  due  to  low  soli  evaporation 
and  temperature  brought  additional  transpiration.  Hence  Up  Increased 
with  temperature.  Early  in  the  growing  cycle,  the  combination  of 
relatively  similar  Kp  among  all  treatments,  but  lowering  Kj,  caused  a 
net  reduction  in  Kp,  Uhile  at  high  LAI,  the  Increase  in  Kp  was  compe- 
nsated by  decrease  in  Kj  and  Kept  Kp  constant.  Therefore  with  changes 
in  tenperature  Kp  was  aimpst  constant,  although  Kp  and  Kj  changed 
during  the  peak  growing  season. 

Vapor  Pressure  Deficit  (VPDl 

The  Fenmah  method  is  not  very  sensitive  to  changes  in  VPO  in  the 
humid  climate.  Low  VPO  is  accompanied  by  low  relative  humidity,  and 
according  to  Eqn.  (5.3)  lower  net  radiation.  CUe  to  these  counteract- 


F1g.  5.3  Effect  of  tempersture  on  water  use  and  crop 
coefficients 

a)  Variation  In  daily  water  use  conponents  under 
various  conOinatlons  of  maxinian  and  nlninun  air 
temperature  in  ‘C.  Top  (27,17),  niiadle  (32,22). 
and  bottom  (37,27) 

b)  Comparison  of  crop  coefficients  due  to 
cnanges  In  natimum  and  mlninim  air  temperatures. 
Solid  line  (27.17),  dotted  line  (32,22).  and  dot 
and  dasn  combination  (37,27). 
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Ifig  effects,  tfie  seasondl  Increeserl  by  merely  6%  for  e three  fold 
Increese  in  the  VPD  (Table  S.S)  Changes  In  the  VPO  showed  some  intere- 
sting changes  in  the  water  use  and  crop  coefficients  (fig.  5.4).  The 
soil  evaporation  increased  more  so  than  transpiration  with  Increase  In 
VPD.  When  the  canopy  was  small,  energy  was  not  the  limiting  factor 
for  soil  evaporation  and  although  VPO  increased,  it’s  influence  was 
minor.  However,  as  the  canopy  developed  and  energy  became  limiting, 
VPO  was  more  important  and  caused  increased  soil  evaporation,  A three 
fold  increase  in  VPO  increased  seasonal  soil  evaporation  by  48*  and  T 
by  9*  compared  with  the  standard  weather  case.  Host  of  this  Increase 
occurred  during  full  canopy.  For  transpiration,  on  the  otherhand  at 
low  LM  ( high  canopy  stomatal  resistance).  VPO  increased  transpira- 
tion. wnile  at  high  LAI  It  had  very  little  Influence.  Also  at  high 
LAI.  the  effect  of  the  Increase  in  VPD  was  neutralized  by  higher  eva- 
poration. Similar  behavior  was  seen  among  all  irrigation  treatments. 

The  crop  coefficients  are  shown  in  Fig.  5.40.  K,  remained  relat- 
ively constant  at  high  LAI  under  alt  three  Irrigation  frequencies 
whereas  Aj  increased  with  VPO.  At  low  LAI  both  Kg  and  K,  Increased  as 
transpiration  and  sol!  evaporation  increased  with  VPO.  The  compensa- 
tory effect  of  soil  evaporation  on  transpiration  was  neutralized  by 
higher  VPD.  Therefore  was  high  for  high  VPO  treatment  by  as  much 
as  15  to  20  % throughout  the  growing  period.  These  findings  suggests 
chat  determined  under  humid  climatic  conditions  will  underpredict 
water  needs  of  a dry  region  and  vice  versa.  However.  If  the  surface 
wetting  Is  minimized  then  adjustment  in  or  Kg  may  not  he  required. 
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Cumulative  reference  crop  evapotranspiration  (cETr), 
sifflulateO  eoybean  evapotranspiration  (cET),  soil 
evaporation  (cE)  and  transpiration  fcT),  for  three  VPD 
levels  under  dally,  once  every  three  days  and  rteHy 
I rrigation  cycles. 


I day  cycle  3 day  cycle  7 day  cycle 


7.5  15  22.5 


550  5ll  6^) 
627  6>18  752 
171  209  254 
456  479  498 


7.5  15  22.5 

590"  611  6^3 
611  667  725 
148  182  221 
463  405  504 


Effects  of  vapor  pressure  deftcU  on  water 
and  crop  coefficients 


a)  Variation  in  daily  water  use  components  under 
various  vapor  pressure  deficits  in  mb.  Top 
(7.5),  middle  (15),  and  bottom  (2Z.5) 


b)  Comparison  of  crop  coefficients  due  to 
changes  in  vapor  pressure  deficits.  Solid  line 
(7.5  mb),  dotted  line  (15  mb),  and  dot  and  dash 
combination  (22.5  mb) 
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EFFECT  OF  TPD  ON  WATER  USE  AND  CROP  COEPnOENTS 


3 Bir  CYCLE  7 0»Y  CYCLE 


US 

Radiation 

Increaset  or  decreases  In  radiation  caused  proportional  rises  and 
falls  in  dailji  and  seasonal  ET  and  ET^  (see  Fig.  5. SB  and  Table 
S.6}.  Soil  evaporation  during  the  incofnpiete  canopy  stage  was  also 
higher  initially  hut  decreased  sharply  under  high  radiation  condi- 
tions. The  effect  of  increased  radiation  on  soil  evaporation  was  less 
significant  at  longer  Irrigation  intervals.  Sinilarly,  transpiration 
increased  more  rapidly  with  radiation  and  was  little  influenced  By  the 
irrigation  frequency.  M low  LAI  and  high  radiation,  soil  evaporation 
was  limited  by  available  water,  while  at  high  LAI  it  was  limited  by 
both  the  energy  and  the  humid  environment  in  the  canopy.  As  a result 
Ks  decreased  with  increasing  radiation.  On  the  other  hand,  lower 
evaporation  at  fu11  canopy  was  compensated  by  Increased 
transpiration.  Therefore  H was  similar  among  all  treatments  early  in 
the  season  while  it  was  higher  during  the  peat  growing  season.  Hue  to 
the  counteracting  behavior  of  Kj  and  K„  at  low  LAI,  K was  lower  for 
high  radiation  but  at  high  LAI,  radiation  did  not  mate  any 
difference.  Therefore  from  a practical  point  of  view,  changes  in 
radiation  could  Influence  considerably  during  the  incomplete  canopy 
stage  but  vrauld  nave  little  effect  during  tne  peak  growing  season. 

Kind  Speed 

Wind  speed  Influenced  Che  canopy  tP  hulk  air  boundary  layer  resi- 
stance as  well  as  the  soil  boundary  layer  resistance.  The  boundary 
layer  resistance  Is  inversely  proportional  to  Che  wind  speed,  there- 
fore both  the  soil  evaporation  and  transpiration  increased  at  higher 
wind  speed  (Table  6.7  and  Fig.  5.6).  The  effect  of  higher  wind  speed 
on  soil  evaporation  was  noticeable  during  full  canopy  stage  when 


Effect  of  ridtatfori  on  water  use  and  crop 
coefficients 


a)  Variation  1n  dally  water  use  components  under 
various  radiation  conditions.  In  ca1/cn2/day. 
Top  (Z22),  middle  (444),  and  bottom  (722). 


b)  Comparison  of  crop  coefficients  due  to 
changes  in  radiation.  Solid  line  (222),  dotted 
line  (444),  and  dot  and  dasfl  cofflOination  (722). 
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Curnul  acive  reference  crop  evapocransplretlon  (cETr). 
sinulated  soybean  e«apotranspiratfon  (cET),  soil 
evaporation  |cE)  and  transpiration  [cT),  for  tbree 
radiation  levels  under  daily,  once  every  three  days  and 
eeekly  irrigation  cycles. 


1 day  cycle 

3 day  cycle 

7 day  cycle 

1 ' ^1 

RAD > 

2ZE  444  m 

222'  444  '?22 

222  444  722 

ce!  m 
cT.  nm 

4DS  627  904 
121  171  224 
2S4  4S6  680 

402  611  863 
117  148  164 
2SS  463  699 

374  m 86Q 
389  535  833 
IDO  114  120 
289  471  713 

1S2 


TABLE  5.7  Cumulative  reference  crop  evapotransfil  ration  (cETr), 
sinulaced  soybean  evapotranspiration  IcET).  soil 
evaporation  (cE)  anil  transpiration  (cT).  for  three  Kind 
speeds  under  dally,  once  every  three  days  and  weekly 
1 rrigation  cycles. 


1 day  cycle 
4 259.2  518.4 


cETr,  im 
cET,  ™ 
cE,  nun 
cT.  ran 


456 


590  661  767 
585  638  696 
114  146  174 
471  492  522 


fig.  5.6  Effect  of  wind  sped  on  eater  use  and  crt^ 
coefficients 

a)  Variation  in  daily  water  use  components  under 
various  wind  speeds,  in  Kn/day.  Top  (B6.4), 
middle  (259.2),  and  bottom  (S18.4). 

0)  Ccvnparison  of  crop  coefficients  due  to 
changes  In  wind  speeds.  Solid  line  (86.4), 
dotted  line  (259.2),  and  dot  and  dash 
combination  (518.4). 
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energy  wae  Hmitfng  soil  evaporation.  U1th  a daily  Irrigation  cycle, 
seasonal  soil  evaporation  was  more  than  doubted  at  wind  speeds  of 
518. a Kn/day  white  transpiration  increased  by  merely  9%.  The 
Influence  of  wind  speed  on  soil  evaporation  reduced  with  length  of 
Irrigation  cycle  as  soil  began  to  limit  E.  The  feedbacn  and  compensa- 
tory effects  between  soil  evaporation  and  transpiration  was  worth 
notihg  An  Increase  in  soil  evaporation  with  wind  speed  Increased  A,, 
hot  decreased  Kp.  However  as  the  Increase  in  Kj  was  greater  than  the 
decrease  in  Kp.  increased  during  the  daily  Irrigation  treatment. 
As  the  time  between  irrigation  events  increased,  soil  became  the  E 
Uniting  factor  at  low  LAI,  and  Kp  was  lower. 

In  sunvnary,  radiation,  VQO,  and  wind  speed  affected  Kp.  Radia- 
tion had  large  impact  early  in  Che  season  more  so  than  VPD.  Late  In 
the  season,  Kp  were  insensitive  to  radiation  but  increased  with  VPO  by 
as  nejch  as  15  to  20  1.  Hind  speed  affected  Kp  throughoot  the  growing 
season  for  1 and  3 day  cycles,  where  as  during  the  7 day  cycle  was 
affected  only  early  in  the  season. 

Planting  Dates 

Management  factors  affecting  the  value  of  crop  coefficient  are 
mainly  the  crop  characteristics,  crop  planting  dates,  rate  of  crop 
development,  and  length  of  growing  season.  In  the  earlier  section 
the  effects  of  climatic  factors  was  considered.  In  this  section 
results  of  using  actual  data  with  spring  planting,  summer  planting  and 
fall  planting  are  analyzed  to  evaluate  the  effect  of  soybean  develop- 
ment and  length  of  growing  season  on  water  demand  and  crop 


coefficients. 
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Uilliams  soybeans  pUrted  during  larch  (naturUy  Group  III  ) had 
slow  vogecative  growth  due  to  relatively  lower  ambient  temperatures. 
By  contrast  temperatures  were  high  for  the  June  and  August  plantings 
and  the  canopies  attained  ntazlmum  LAI  within  60%  of  growing  season 
compared  to  661  for  the  March  planting  (Fig.  6.1).  It  is  Interesting 
to  note  that  although  both  the  June  and  August  plantings  dates  reached 
niaaimun  LAI  at  the  same  relative  times,  their  magnitudes  were  signifi- 
cantly different  doe  to  the  photoperiod  response  of  the  cultivar.  Due 
to  snort  days,  soybeans  planted  on  August  15  completed  their  Hfe 
cycle  in  UO  days  compared  to  13A  days  for  the  June  planting  and  117 
days  for  the  March  planting.  The  dally  ET,  E,  T and  ET,.  for  all  three 
planting  dates  under  the  three  irrigation  frequencies  is  shown  In  Fig, 
5.7a.  The  daily  rate  of  transpiration  increase  was  highest  with  the 
June  planting  date  and  slowest  witn  the  March  planting  date,  while  for 
the  August  planting  it  was  intermediate.  Put  it  reached  a lower  peak 
value  due  to  low  LAI.  The  trend  Ih  the  soil  evaporation  rates  were 
opposite  that  of  transpiration. 

Under  all  planting  dates  seasonal  soil  evaporation  decreased  with 
the  time  between  irrigations  and  transpiration  Increased  (Table 
5.8).  Seasonal  water  use  revealed  that,  with  the  August  planting  and 
daily  irnigatlon,  soil  evaporation  eiceeded  transpiration  due  to  lower 
LAI  development  compared  with  the  March  and  June  planting.  The  incom- 
ing radiation  Increased  with  time  during  the  March  plahclng,  while 
radiation  decreased  for  tne  other  two  planting  dates  (see  Table 
5.2].  As  a result  ET^  steadily  Increased  for  the  March  planting  and 
decreased  during  the  June  and  August  planting  dates. 

The  crop  coefficients  were  significantly  different  among  planting 


Fig.  6.7  Effect  of  pUntIng  Pates  on  water  use  and  crop 
coefficients 

a)  Variation  in  daily  water  use  rates  under 
different  planting  dates;  top  (Mar  IS),  middle 
(Jun  26),  and  bottom  (Aug  16). 

b)  Variation  in  crop  coefficients  under  to 
different  planting  dates. 
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EFFECT  OF  PUHTIHG  DUES  ON  tilEB  USE  IKD  CHOP  COEFFICIENTS 
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TABLE 


Cumulative  reference  crop  evapotransplratlon  (cETr). 
simulated  soySean  evapotransplratlon  (cET),  soil 
evaporation  (cE)  and  transpiration  (cT) . for  three 
planting  dates  under  dally,  once  every  three  days  and 
weekly  Irrigation  cycles. 


1 day  cycle 

HAft  JUN  AUG 

IS  26  IS 

cETr,  i-tn  507  511  355 

cET,  rm  514  651  338 

cE.  ran  210  179  181 

cT.  inn  304  372  157 


3 day  cycle  7 day  cycle 


MAR  JUN 

507  511 
442  S04 
125  117 
317  387 


AUG 

15 
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Oates  (Fig  S.7b).  ts  t»ie  tleie  between  irrigations  increased,  E ano 
consequently  £T  was  reduced.  TiiereFore  according  to  Eqn  (5.2)  both 
tde  and  varied  with  planting  date.  Since  E was  always  less  than 
ET,..  etcept  for  the  first  few  days  under  the  dally  Irrigation  cycle, 
Kj  was  always  less  than  1.  The  relative  magnitude  of  before  the 
time  of  maximum  LAI  among  different  planting  dates  was  directly 
related  to  the  soli  evaporation  rates  (or  inversely  related  to  the 
rate  of  LAI  development).  Aj,  for  the  August  planting  date,  due  to 
low  LAI,  exceeded  the  other  two  planting.  Kp  For  the  June  planting 
date  was  consistently  higher  then  the  other  planting  dates  due  to 
higher  LAI  and  gradually  decreasing  ET,..  Peak  T rates  for  the  June 
and  August  planting  occurred  at  the  same  time  during  Che  growing 
season,  while  it  was  delayed  for  the  March  planting.  This  behavior 
was  duplicated  by  Ap  values  which  also  peaked  at  different  times. 

ET  peaked  at  different  times  during  each  planting  depending  upon 
the  magnitudes  of  T and  £.  Therefore  Kp  also  varied  accordingly.  For 
the  daily  Irrigation  cycle  was  almost  constant  between  planting 
and  physiological  maturity  irrespective  of  the  planting  dates. 
However  as  the  time  between  irrigations  Increased,  Kp  increased 
gradually  before  attaining  a relatively  staole  value.  Therefore  the 
magnitude  of  kp  was  primarily  governed  by  E during  Incomplete  canopy 
state,  wM1e  It  depended  more  on  T as  the  canopy  attained  full  cover. 

In  summary,  planting  dates  could  change  the  magnitude  and 
relative  behavior  of  crop  coefficients  in  the  dynamic  climate  of 
Florida.  The  rata  of  development  as  well  as  variations  In  the  ETp 


governed  the  magnitude 


ConclniUnq  Remarks 


Most  of  the  crop  water  requirements  have  been  aetermined  by 
specific  field  experiments  (DoorenOos  and  Pruitt,  1977],  These 
studies  provided  information  on  water  required  for  ^T  when  crop  growth 
was  not  limited  by  soil  water.  Extrapolation  of  these  published  coef- 
ficients to  different  environments  has  been  practiced  repeatedly 
wftnout  knowing  how  vaiid  they  are  for  these  different  environments. 
In  this  paper,  the  modei  deveioped  in  Chapter  III  was  used  to  study 
the  effect  of  varying  climate  and  management  factors  on  water  require- 
ments and  crop  coefficients  for  a soybean  crop  growing  on  a sandy 
soli.  A sensitivity  analysis  of  crop  coefficients  to  changes  in  temp- 
erature. radiation,  vapor  pressure  deficit,  and  wind  speed  was 
performed.  Since  the  method  computes  soii  evaporation  and  transpira- 
tion separately,  the  overall  crop  coefficient  K^,  was  divided  into 
soil  evaporation  coefficient  and  transpi ration  coefficient  hp. 
Tne  analysis  showed  that  Kp,  the  transpiration  coefficient  was  insens- 
itive to  climatic  changes  when  compared  to  Kp  and  K«,  and  K, 
changed  with  variation  in  radiation,  VPf),  wind  speed  and  irrigation 
management  practice  but  were  insensitive  to  temperature  changes. 

Simulation  with  actual  weather  data  during  three  growing  seasons, 
spring,  sumer  and  fall  showed  some  interesting  results.  Due  to 
changes  in  the  ambient  air  temperature  and  photoperlod  response  of 
soybean  cultivars,  different  canopy  development  rates  occurred.  As  a 
result,  Kp,  Rq,  and  all  varied  with  planting  date. 

Irrigation  frequency  is  a management  practice  that  can  nave  sign- 
ificant impact  on  water  use  and  crop  coefficients  for  the  sandy 
soil.  Seasonal  water  use  for  tne  daily  irrigation  was  predicted  to  be 
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6 i rtigrter  thsn  that  for  the  7-flay  Interval  [ Taols  5.1) 
Irrigations.  Croo  coefficients  near  planting  redoceO  front  about  1.0 
for  the  daily  Irrigation  cycle  to  less  than  0.4  for  the  7-day  cycle. 
This  difference  dissipated  as  LAI  Increased  and  for  all  cycles  were 
similar  after  full  canopy.  Therefore  In  planning  for  water  needs  of  a 
crop,  Irrigation  frequency  Is  very  Important.  Crop  managers  could 
perhaps  conserve  water  for  crops  by  developing  strategies  chat 
maintain  adequate  water  for  transpiration  but  reduce  soil 
evaporation.  Possible  methods  Include  longer  Irrigation  Intervals, 
partial  soil  surface  wetting,  or  Irrigation  below  the  soil  surface. 

There  Is  growing  evidence  that  computer  models  relating  yield  to 
transpiration  are  more  sound  than  those  relating  yield  to  evapotrans- 
piraclon  (Hanks  and  Hill.  1980),  Therefore  the  model  used  in  this 
study  when  Integrated  with  a crop  growth  model  could  enhance  Us  appl- 
ication potential.  The  weather  Inputs  required  for  Che  model  are 
similar  to  many  crop  growth  models,  while  the  drop  parameters  are 
usually  estimated  during  the  simulation.  However  the  time  step  needed 
by  this  model  may  be  of  some  disadvantage.  Tbe  monel  needed  1.8  CPU 
seconds  to  simulate  one  day  on  a P9IHE  550  computer.  Ifhen  using  this 
model  with  other  models,  Che  CPU  requirements  will  be  less  since  many 
of  the  inputs/outpucs  would  become  unnecessary.  The  entire  model  Is 
less  Chat  1000  lines  of  code  and  can  be  run  on  a micro-computer. 


CHAPTER  VI 

SUHMART  AKI5  CONCLUSIONS 

Tne  need  to  predict  evaporation  and  transpiration  components  of  a 
developing  crop  was  identified  as  an  important  issue  in  determining 
appropriate  water  management  strategies.  Evapotranspi ration  for  a 
specific  crop  is  governed  by  prevailing  climate,  crop  and  stage  of 
growth,  and  availability  of  soil  water  content.  In  this  study  a me- 
chanistic model  of  evapotranspiration  for  well-watered  crop,  where  the 
soil  may  go  through  cycles  of  wetting  and  drying,  was  developed  based 
on  the  physical  characteristics  of  soil,  plant  and  atmosphere. 

A physical  model  of  soil  evaporation  was  developed  by  assuming 
that  soil  evaporation  was  occurring  from  wet  soil  below  a dry  soil  of 
increasing  thickness.  A ten-  week  soil  evaporation  study  was  condu- 
cted in  outdoor  controHed-environment  chambers  to  determine  an 
empirical  function  of  soil  resistance  with  drying.  Cumulative  soil 
evaporation  over  a 7*day  period  was  less  under  high  radiation, 
cloudless  conditions  than  for  cloudy  conditions.  This  result 
demonstrated  the  importance  of  sensible  heat  transfer  in  forming  a 
barrier  to  evaporation  on  the  sandy  soil. 

The  soil  evaporation  model  was  integrated  with  a two  aone  canopy 
model  to  form  a ccmiprehensive  evapotranspiration  model  valid  for  any 
stage  of  crop  growth.  The  canopy  was  divided  into  a zone  of  sun-lit 
leaves  and  a zone  of  shaded  leaves.  Algorithms  to  allow  this  distinc- 
tion were  presented.  The  proposed  model  included  a number  of  rarely 
measured  canopy  and  spil  attributes,  so  solution  schemes  were  develop- 
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ed  using  noraially  measured  cllmace,  crop  and  soil  variaQle  so  that  Ctie 
model  could  be  used  for  practical  applications  and  evaTuations, 

Bragg  soybeans  were  groan  in  330.  660  and  990  COg  concentration 
levels  in  three  outdoor,  envi ronmental  controlled  (temperature  and 
hunidity  ) chambers.  Hater  use  (evaporation  and  transpiration) , plant 
groatn  cnaracterlstics,  and  weather  inputs  were  monitored.  The  simu- 
lated results  agreed  well  with  the  measurements.  The  average  daily 
errors  were  less  than  61  . The  cumulative  leaf  area  developed  faster 
and  water  use  decreased  with  Increasing  COj  levels.  The  sunlit  and 
shaded  leaf  temperature  increased  with  COj  levels. 

Simulation  studies  were  conducted  to  study  the  effect  of  COj 
induced  increases  in  stomatal  resistance,  leaf  area  development  and 
air  temperature  on  crop  water  use.  Seasonal  water  use  was  reduced  at 
high  COj  concentration  If  the  Increase  in  air  temperature  projected 
for  higher  COj  environment  was  ignored.  Sowever,  water  use  under 
high  CO2  conditions  was  predicted  to  Increase  under  a combined  effect 
of  higher  stomatal  resistance,  faster  leaf  area  development,  and  a 5 
®C  higher  air  temperature. 

In  another  simulation  study,  variation  in  soybean  crop 
coefficient  due  to  climatic  changes  were  studied.  Pmong  the  climate 
factors  considered,  water  use  and  crop  coefficients  were  most 
sensitive  to  Increased  radiation  and  lowering  relative  humidity.  The 
overall  crop  coefficient,  was  divided  into  a soil  evaporation 
coefficient  and  a transpiration  coefficient.  For  a fixed  crop 
development,  the  transpi ration  coefficient  was  Insensitive  to  changes 
in  climate  compared  to  soil  evaporation  and  overall  crop 
coefficient.  However,  all  crop  coefficients  variefl  when  soybeans  were 
grown  at  different  tines  of  the  year. 


cliinate  dnd  management 


This  study  showed  that  varied  with 
factors.  Use  of  published  values  for  varying  cHnate  and 

management  may  over  or  under  predict  crop  water  needs.  The 

application  model  demonstrated  how  values  can  he  adjusted  to  more 
closely  fit  local  weather  and  management  conditions  when  better 
accuracy  is  needed. 

The  model  developed  in  this  study  will  provide  an  useful  tool  to 
study  water  use  by  a crop  under  a variety  of  climatic,  soil,  plant  and 
management  scenarios.  The  model  can  provide  a useful  research  tool  in 
Identifying  important  and  critical  research  needs  by  eliminating  less 
important  situations.  Dynamic  crop  growth  models  are  being  developed 
to  predict  growth  and  yield  of  crops.  These  models  are  attractive 
because  crop  growth  stresses  due  to  pest,  nutrient  and  or  water  can  be 
included  to  provide  a more  comprehensive  tool  for  crop  production 
management.  An  evapotranspiratlon  model  such  as  the  one  developed 
here  will  greatly  enhance  application  potentials  of  crop  simulation 
models. 

A model  of  this  nature  opens  doors  for  tOe  prospects  of  automated 
Irrigetion  scheduling.  A crop  growth  model  integrated  with  this  model 
will  simulate  plant  growth  and  water  requirmnents  with  inputs  obtained 
from  sensors  In  remote  fields.  WitO  innovations  in  the  integrated 
circuits,  microprocessors  and  computer  technology,  a farmer  can 
install  an  automated  weather  and  crop  reporting  system  for  less  than 
toe  cost  of  a garden  tractor.  Such  self  sustained  systems  are  already 
providing  weather,  crop  and  soil  information  necessary  for  the  model 
implementations. 


■5?PEN0IX  « 

OESCRIPTION  OF  SUBROUTINES 


FUEOI 

S4USSJ 


INir 

IRIBAT 

NATADB 

MATPRO 

OUTPUT 


SniLTEMP 

SUNANB 

SOLUTION 


This  is  the  nein  driver  program. 

Computes  soil  houndary  layer  and  canopy 
aerodynamic  resistance  to  heat  and  vapor  for  hare 
soil  and  fully  developed  canopy.  Nain  program 
uses  these  values  and  compotes  reslstancves  for 
intermediate  growth  stages. 

Reads  time  Invariant  inputs  needed  by  the  model. 
This  rdutine  is  called  repeatedly  by  suhroutine 
SOLUTION  to  solve  a system  of  three  simultaneous 
1 inear  eguations. 

Initialises  variables  each  day 

Updates  soil  water  content,  and  soil  temperatures, 
fpllowing  an  irrigation. 

Natrih  addition  routine 
Matrix  product  routine 

Rdutine  to  integrate  diurnal  flutes  of  sensible, 
latent,  and  soil  heat  fluxes  during  the  season  and 
write  diuma!  and  daily  outputs. 

Routine  to  comoute  soil  temperature  by  layers. 
Computes  solar  elevation  as  a function  of 
geographical  co-ordinates,  month,  and  time  of  day. 
Routine  to  compute  sensible  and  latent  heat 
fluxes,  and  soil  surface  temperatures. 


19h 


19? 


'Ok-up  function 


Function  to  compute  setureterj 
elr  at  a given  tenpereture 


pressure 


APPENDS  8 

DEFINITION  OF  VARIABLES  IISEO  IN  T8E  COMPUTER  COOE 

VARIARLE  NAME  [ VARIABLE  TYPE,  LINKAOE  (IF  LOCAL  VARIABLE)  OR  /COMiflN 
STATEMENT  NAHE/,SU8R0UT1NE(S)  IN  IIHICK  VARIABLE  IS  USED, 
SUBROUTINE  cons  ] • DEFINITION  Of  VARIARLE 


COOE  VARIABLE  TYPE 


I - SHORT  INTEGER(INTEOER*2) 

J - LONG  INTEGER  (INTEGERM) 

L - LOGICAL 

R -REAL  (REALM) 

TNe  following  code  eupears  after  the  suNrouclne  name  In  the  glossary; 
CODE  SUBROUTINE  CODE 

A - VARIABLE  IS  CONTAINED  IN  THE  ARGUMENT  LIST  OF  THIS  FUNCTION  OR 
SUBROUTINE 

D - VARIABLE  IS  DEFINED  IN  THIS  FUNCTION  OR  SUBROUTINE 

I - VARIARLE  IS  INITIALIZED  IN  THIS  FUNCTION  OR  SUBROUTINE  ( RATA 
STATEMENT) 

M - VARIABLE  IS  MODIFIED  IN  ONE  OR  MORE  ASSIGNMENT  STATEMENTS  OF 
THIS  FIINOTION  OR  SUBROUTINE 

S - VARIABLE  IS  IN  A DATA-MOOE  SPECIFICATION  STATEMENT  OF  THIS 
FUNCTION  OR  SUBROUTINE 

ALPHA  R LINKAGE  MAIN  I - Constant  in  the  soil  resistance  equation 
to  compute  effect  of  net  radiation  on  soil  resistance. 


BETA  R LINKAGE  MAIN  A • Solar  elevation  angle  in  radians 


CAP  R LINKAGE/  FIXEDl  M:IfiIGET  M ; SOILTEMP  M - Heat  capacity  of 
soil  at  certain  soil  water  content  (J/m3-  C) 
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: - No.  of 


or  soft 


S M:  SOLUTION  S : PsychPO(*ietr1c 


(P*/C1 


liL  » /NISC/  NAIN  S I;  OUTPUT  S N : Soil  Heat  flux  from  last 
Simulation  step 

MCAP  R /STMP/  FIXEfll  SNA:  IRIGET  S A : SOILTEHP  S A ; Soil 

Heat  capacity  {J/m3-C) 

HEC  R /HISC/  MAIN  S:  OUTPUT  S M : Hourly  transpiration,  im/H 

HES  R /MISC/  MAIN  S:  OUTPUT  S M : Hourly  evaporation,  mm/n 

lERROR  J ARRUMENT  SAUSSJ  S N:  SOLUTION  A ; Error  code 
ITER  J LINKAGE  SOLUTION  M:  Iteration  number 
JULIAN  J /UETH/  NAIN  S N:  HTHR  SNA:  Julian  day 

KET  R /fIXI/  FIXEOI  S M:  SOLUTION  5 : KEY-1  soil  surface 

temperature  is  an  input;  KEY-2  soil  surface  temperature  need 
to  Oe  computed  in  tne  model 

KK  R LINKAGE  SOLUTION  M ; Error  code 

LAI  R /MISC/  MAIN  S M:  OUTPUT  5 : Leaf  area  index 

LAIOL  R LINKAGE  MAIN  N:  LAI  of  shaded  leaves 

LAIMAX  R /Fill/  FIXEOI  S M:  MAIN  S : LAI  when  canopy  acnieves  full 


LAISL  R LINKAGE  IIAIN  N:  LAI  of  sunlit  leaves 

LDAILY  R /FIXI/  FIXEOI  S M:  OUTPUT  S ; flag  for  dally  output 
C LOAILY-1  i.e.  YES.  LDAILY-1  I.e.  NO  j 

LMOUfi  R /FIXI/  FIXEOI  S M:  OUTPUT  S : Flag  for  Hourly  output 
r LNOUR-1  I.e.  YES.  LH0UR--1  i.e.  NO  ] 

LHV  R /MISC/  MAIN  S N:  OUTPIJT  S : MTHR  S H : Latent  Heat  of 
vaporUatlon  of  water  (J/gm] 

LL  R /IRR/  FIXEOI  S M:  MAIN  S : Lower  limit  of  soil  water 
extraction,  volumetric  basis. 

NDATAP  J /METH/  MAIN  S:  MTHR  S H : Time  steps  per  day 

NELAVR  J /FIXI/  FIXEOI  5 M;  MAIN  S : NumOer  of  soil  layers 

(counted  from  soil  surface)  from  which  water  can  be  lost  as 
soil  evaporation. 

NFLOH  R LINKAGE  SOILTEMP  S II:  Net  Heat  flow  from  a soil  layer 
(J/s-C) 


WINU  INIT  S M;  OUIOUT  S H : Daily  sun  of  sensible  heat 

,".;s  s 
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WTHR  S H : Aip  temoerature,  C 

TAVG  R LINICAGE  MAIN  M A ; Average  of  soil  and  afr  tefngeratvre  C4 
compute  slope  of  saturated  vapor  pressure  vs  teaiperature 


TC  fi  /SGL/  OUTPUT  S:  SOLUTION  S M : >dean  canopy  temperature,  C 

TD  R /WETH/  MAIN  S A:  SOLUTION  S A ; UTHR  S H : Oe»  point 

temperature,  0 

TEMP  S /STNP/  FIXEOI  S H:  IRIGET  S M ; SOLUTION  S : SOILTEMP  S 
1 :Tenperature  of  soil  by  layers,  C 

TNCONO  R ySTNP/  flXEDI  S M;  SOLUTION  S A ; SOILTEHP  S A ; Themal 
conductivity  of  soil,  J/m-s-C. 

THETA  R /STMP/  FUEOI  S H A:  IRIGET  S A : SOLUTION  S A : 

SOILTEMP  S A ; Water  content  of  soil  at  wMcti  NCAP  values 


TINE  R /WETH/  MAIN  S A;  OUTPUT  S : UTHR  S '1  : Tine  of  day,  n 

TOTAL  R /WETH/  MAIN  S:  WTHR  S H : Total  daily  inconing  net 

radiation,  mn 

TS  R /SOL/  MAIN  S Mr  OUTPUT  S r SOLUTION  S H : Computed  surface 
temperatures,  C [ TS(l)-  sunlit  leaves.  TS(2)-  snaded  leaves, 
rs(3)-  sol  1 Surface  ] 

TS3MAX  R /SOL/  INIT  S ; OUTPUT  M : riaeimun  soil  surface  temperature 

TS3MN  R /SOL/  INlT  $ : OUTPUT  ri  : Minimum  soil  surface  temperature 

TSESN  R /MISC/  MAIN  S I:  OUTPUT  S M : Transpiration  for  season,  ran 

TSS  R /WETH/  MAIN  S:  HTHR  5 M r Soil  surface  temperature,  C 

TSTEP  R /STMP/  OUTPUT  S M:  SOILTEMP  S r Tine  step  of  simulation,  h 

TSilRF  R /SOL/  SOLUTION  S M:  SOILTEMP  S : Intenmefllate  value  of  soil 
Surface  temperature 

U R /WETH/  WTHR  S M : Mind  speed,  m/s 

UAVG  R /FUI/  FIXEOI  S M:  MAIN  S A ; BLRRES  S : Average  wind 

speed  betmeen  LAI*0  and  LAI«LAIMAx,  m/s 

V R /SOL/  OUTPUT  S;  SOLUTION  S A ; Latent  neat  flux.  U/n2-s 
[S[l)-  sunlit  leaves,  S(2)-  snaded  leaves,  S(3)-  soil  ] 

VHCAR  fi  /STMP/  fUEDI  S H;  IRIGET  S N : SOILTEMP  5 M : Volumetric 
heat  capacity  of  soil,  (J/m3) 


UA7 


zn  li  LINKAGE  SLIRES  N - Crop  rougtiness  Height,  o 

Z1  R LINKAGE  BLRRES  N - Soft  roughness  neignc,  m 


UPfENDI*  C 

SOURCE  CODE  OF  SUBROUTINES 


C***  tlAIN  PROGRAM 
SSTORAGE:E 
IMOFUOATCAUS 
SINCLUOE;  'CONHON.DAT' 

DATA  ALPHA. RITA, ETAn.094PE-3. 3. 01R7E-3, 1.36757/ 


INPUT  SECTION 


OPEN  (UNIT  = 6 .FILE 
OPEN  (UNIT  ■ 7 .FILE 
OPEN  (UNIT  » 8 .FILE 
OPEN  (UNIT  ■ 9 .FILE 
OPEN  (UNIT  - in, FILE 


INPUT. OAT' .STATUS  ■ 'OLO‘) 
WEATHER. DAT'. STATUS  » 'OLD') 
PLANT.DAT' .STATUS  ■ 'OLD') 
OUTPUTl. OAT'. STATUS  • 'NEH'l 
0UTPUT2.DAT'. STATUS  • 'NEW') 


DAIR  - 1.1AOOE3 
J/jir-C 


3SILQN  • O.B22 

RNTOT  • 0. 

TSESN  • 0. 

ESESN  - 0. 

ETSE5N  > 0. 

RNNET  • 0. 

CALL  FIKEOI 

CALL  eLRRES(CHIGHT,Uj 


3.RSl.RS2,fiCl,RC2) 


.END  • 1000)  JULIAN, Lj 


IF  (OAT  .GT.  1)  no  TO  200 
TS(3.1)  . TA(1) 

G ■ RN0(1)*EXP(-E)!TINC*LAI) 

0 RNTOT  ■ RNTOTtRNNET 

IF  (ORYDAT  .LE.  FREQNC  .OR.  OAPL  .SE.  NR7)  G 
m IF  (LAI  .ST.  LAINAX)  GO  TO  400 
58N0RY  • RS1*(RS2-RS1)*LAI/LAINAX 
R ■ RC1*(RC2-RCI)*LAI/LAII1AX 
0 RNNET  ■ TOTAL*EX'>(-EtTINC*LAI) 

RNA  . (RNTOT  ♦ RNNET)/  ORYOAY 
IF  (RNA  .LT.  RNAHIN)  RNA  ■ RNANIN 
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CULL  INIT 


TLAST  . riHEU) 
M 9f-  ' - - 


■ 2500.9  - 2.345*TA(NHDUR) 

SAW  = P*CP/(LHV*aSILON] 

TAVG  ■ (TS(3.n*rA(NH0UR))*0.6 
DELTA  - VPDF(TAVS+1.)-VPDF(TAVG) 

CALL  5UNANG(XLAT,XL0KG.NH0NTH.T1«(NH0UR),BETA,AZANGL) 
59ETA  • SIN(BETA) 

IF  (SBETA  .LE.  01  SBETA  • 1. 

XX  " EXTIKC/SBETA 

LAISL  ■ (l./XX)*(l.  - EXP(-XX*LAI)) 

LAIDL  ■ LAI  . LAISL 
RB(U  • RMAX 

R3(3)  • SBNORX 

IF  (LAISL  .ST.  O.QODl)  RB(1)  ■ RA/LAISL 
IF  [LAIDL  .GT.  0.00011  RB(21  = RA/LAIOL 
RS(I)  . RnAX 
RS{2)  ■ RMAX 

CSL  ■ CMAX'(1.-EXP(SLOPEOPFO(IIHOUR)/CHAX)) 

COL  ■ CHAX*(1.-EXP(SLOPE*PPFD(NHOUR)*FRACTN/C'1AX)) 

IF  (LAISL  .OT.  0.001  .AND.  CSL  .GT.  0)  RS(1)  ■ 
l./rCSL'LAlSL) 

IF  (LAIDL  .GT.  0.001  .AND.  CDL  .GT.  0)  RS(2)  ■ 
l./(CDL*LAIOL) 

RS(3)  ■ 0. 

CEN  ■ (C£/CEH)*10n. 

IF  (CEN  .OT.  100)  CEN  . 100. 

IF  (CEN  .ST.  CEC)  RS(3)  > 

IF  (RS(31  .LT.  0)  RS{3)  • 0. 

IF  (FREQNC  .GT.  99)  RS(3)  ■ RMAX 
RN(l.l)  * RNO(NHOUfi)*XX*LAISL 
RN(2,1)  • RNO(NHOUR)'FRACTII  *XK«LAIOL 
RN(3,1)  • RNO(NHOUR)*£XP(-XK*LAI)-«N(2.1) 

IF  (RN(3,l)  .LT.  0)  RN(3,1)  • 0. 

TS(3,1)  ■ TSS(»HOUR)-TA(NHOUR) 

00  son  1 ■ 1,3 

DO  500  J • 1.3 
RH|I,J)  > R 
RV(I,J)  . R 

IF  (I  .NE.  J)  SO  TO  500 
RH(I,J)  . R .R6(I) 

RV(I.J)  . R tRB(l)*RS(I) 

CONTINUE 

0(1.11  - VP0F(TA(NH01IR1)-XP0F(TD(NH0UR)) 

CONTINUE 


SOLVE  THE  THREE  ZONE  SOIL-PUHT  HQOEL  RT  THIS  TIME 


CRLL  SOILTEHP 

UPDATE  SOIL  MATER 

eVZONE  • HES  ♦ FTEZ'HEC 
CE  • CE  ♦ EVZONE 
EVZONE  = EVZONE‘0.001 
DO  500  IJ  » l.NELAYR 
IF  (EVZONE  .L£.  0)  SO  TO  BOO 
IF  (SM(IJ)  .L£.  LL[IJ))  SO  TO  800 
XTRACT  . (SW(IJ)-LL(IJ))*OUYR(IJ) 
IF  (EVZONE  .LE.  XTRACT)  SO  TO  700 
SM(IJ)  . LL(IO) 

EVZONE  ■ EVZONE-XTRACT 

GO  TO  flOO 

10  SM(IJ)  - SH(IJ)  - EVZONE/OLAYRdJ) 

ifl  CONTINUE 

10  CONTINUE 
CALL  0UTPUT(2) 

DAPL  ■ 0APL*1 
OAY  ■ DAY  t 1 
ORYOAY  • DRYOAYtl 


CLOSE  (S) 
CLOSE  (9) 
CLOSE  (10) 
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SSTOMBE:? 

SNOFLOATCALLS 

SURROUTINE  RLRRES  (CHlGHT,UAtfG,RSl,RS2.RCl .RCJ) 
EN  - 3. 

X ■ 

0 ■ 0.7*CHIGHT**0.979 

ZO  • 0.13*CHIGHT”0.997 

El  ■ E*P(ZN*n.-(D+ZO)/CHIGHT)) 

C E2-EXRtEN*(CHISHT-Zll) 

E2  • EXP(EN) 

E3  = CKIGHT/(EK*1CHIGHT-D) ] 

EA  - AL0G((X-0)/Z0) 

E5  ■ 0.4*n.flniAVG 

SS2  ■ £«*E3*(E2-E1)/E5 

EG  ■ ALOGMX-Dl/ICHISHT-O)) 

RCZ  ■ E4*(E6»E3*(E1-1))/E5 

RSI  • AL0G(X/Zl)*AL0G((ll+20)/Zl)/E5 

RCl  . (AL0G(X/Z1)**2)/E5 

RCl  ■ RCl-RS: 

RETURN 
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55T0RA0e;3 

SNOFLOATCALLS 

SUBROUTINE  FIXEOT 
SHCLUOE:  'COMUOK.OAT' 

READ  (6.')  DAPL,O«mY,CE,LOAILY,LH0UR,FRE0NC,CHIGHT,U!'1AX 
‘ UAVfi.RA.ERRBK0,RMAX.RNA>II«,XUT,«1.0>IS,<Er 

READ  (5,»)  EXTINC.FRACIN.CTAX, SLOPE, NR7 
READ  (6.*)  NELAVR.NTLAYR.CEM.CEC.FTEZ 
REAO  (6,*)  (OLAYR(I).I  . l.NTLAYR) 

READ  (6,*)  (TEMRdl  ,I  . l.NTLAYR) 

READ  (6.*)  (SU(I)  .1  . l.NTLAYR) 

REAO  (6.*)  (SAT[1)  ,1  = l.NTLAYR) 

REAO  (6.*)  [OilL(D  ,1  • l.NTUYR) 

REAO  (S,*)  (LL  (I)  .1  = l.NTLAYR) 

REAO  (6.*)  (THCONO(I),I  - 1.3) 

REAO  (6.«)  ( HCAPdl.I  • 1,3) 

REAO  (6.*)  ( THETA(I).I  • 1.3) 

DO  100  I • l.NTLAYR 
XSN  .Slid) 

CAP  . TA8EL(HCAR.TKETA,XSW.3) 

100  YHCAR(I)  . TENP(I)'CAp.OLAYRd) 

RETURN 
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SSTOASE:2 

MOFLOATCALLS 

SUBRDUTIKE  SAUSSJ(B.  U.  NAflU  ,HAAS,IEI)iiaR) 

C MATRIX  INVERSION  BY  THE  GAOSS  JORDAN  METHOD 


lERROR  LOGICAL 


CURRENT  MATRIX  SIZE. 

DIMENSION  MATRIX  SIZE. 

FLAG  .TRUE.  IF  MATRIX  SINGULAR 


DIMENSION  INOEX(15,3) 

INTEGER  NROU,MAXD,I,J,K,L,IROll,ICOL,Ll 

REAL  B(MAXO.MAXD).  U(MAXD.MAXD) .RIG.SUM.T.PIVOT.OETERH 

on  ino  I ■ 1,  NROu 
INnEX(I,3)  = 0 
inn  CONTINUE 
OETERM  • l.n 

DO  9DD  I • 1.  NROU 
C SEARCH  FOR  LARGEST  ELEMENT. 


IF  (INDEX(X,3)  .GT.  1)  GO 
IF  (IN0£X(X,31  .Eg.  11  GO 
IF  (A0S(8(J,K)1  .LE.  BIG) 
IRON  • J 
ICOL  • K 

BIG  • ABS(B(>I,X)1 
200  CONTINUE 

300  CONTINUE 

INDEX(IC0L.3)  • 1NDEX(IC0L,3)  * 1 
INDEI(l.l)  ’ IROU 
INDEX! 1.2)  ■ ICOL 

INTERCHANGE  RONS  TO  RUT  RIVOT  ON  DIAGONAL. 


: RIVOT  RSI  BY  PIVOT  ELEMENT 

PIVOT  ••  B(ICOL.ICOL) 
DETERH  ■ OETERM  • PIVOT 
8(IOOL,ICOL)  = 1.0 
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C REDUCE  NOK  P: 


IF  (LI  .EO.  ICDL 
T . B(U.ICOL) 
8(U.IC0L)  - 0.0 


B(L1,L)  ■ 8(L1,L)  - 8(IC0L.L)  • T 
700  CONTINUE 

800  CONTINUE 

900  CONTINUE 
END  OF  I LOOP. 

INTERCNANOINS  COLUNNS . 

00  UOO  I ° 1.  HROU 
L • N80W  -1*1 

IF  (INOEJa.l)  .EO.  INDEX(L.2)1  SOTO  1100 
IROU  • lNnEX(L,l) 

ICOL  • IN0EX(L,2) 


CALL  SWAP(B(X.IR( 
CONTINUE 
CONTINUE 

DO  1200  K ■ 1.  NRON 
IF  (INDEX(K.3)  .NE.  V 


, B{K,ICOL)) 


RETURN 
0 lERROR  > 
RETURN 


9 ° HOLD 

RETURN 

END 


SST0SA(!E:2 

SNOFLOATCALLS 

SURROIITINE  IKIT 
JINCLUBE;  'CWWON.WT' 


SHSIH 

ETSINL 

ETSCL 

ETSSL 

SNSCL 

SHSSL 

SHSIML 


TS3MX 

RETURN 
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SST0RARE;2 

SNOFLOATCALLS 

SURROIJTIKE  [RISAT 
SINCLUOE:  ’CWMN.DAT' 

00  11)0  I ■ l.NTLAYR 

CAP  • TA6£L(HCAP,TH£TA.SU(H.3) 

VHCAP(I)  . ( CAP*TE«P{I)  » a.l87E6*(0UUl)-SH(I))*TA(l)  )* 
♦ OLATS(l) 

S0(I)  ■ OUL(l) 

CAP  - TA0EL(HCAP.THETA.SU(O,3) 

TEHP(I)  ■ VHCAP(!)/(CAP-OUYR(I)) 

100  CONTINUE 
CE  • 0. 

RNTOT  = n. 

ORYOAT  • 1. 

RETURN 

END 
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$ST0RAGE:2 

JNOFLOATCAUS 

SUBROUTINE  HATSDO(A,  B,  C.  U.  lY) 

IMPLICIT  »EAL*4(A-H,0-2) 

REAL  A(U.IY),R(IX,IY),C(IX.IY) 

00  100  IROH  ■ l.IX 

DO  lOQ  JCOL  = l.IY 

C(IROU.JCOL)  ■ A(IROW,JCOL)  * 8(IR0W.JC0L) 
ino  CONTINUE 

RETURN 
END 
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SSTORME:? 

SNOFLOATCALLS 

SUBROUTIKE  MTPRO(A,  R,  C,  IX.  lY.IZ) 
REAL  A(IX.IY),B(1Y,IZ).C(IX.IZ) 

60  300  NCOL  • l,tZ 

DO  200  IRDU  • l.IX 

BO'lflO  JCOL  • l.It 
SUM  = SUM  + A(IROW.JCOL) 
100  CONTINUE 

CIIROW.NCOL)  ■ SUM 
200  CONTINUE 

300  CONTINUE 
RETURN 


R(JCOL.NCOL) 
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SST0RAGE:2 

INOFLOATCALLS 

SUBROUTINE  OUTPUT(III) 
tlNOLUOE:  'COtlttON.DAT' 

C 

C 

IF  (III  .£5.  2)  GO  TO  300 
TSTEP  • (TIHE(NH0UR)-TLAST)*3S00. 

NEC  ■ (£TSCUV(1,1)*V(2,11)*0.5*TSTEP/(1000.*LK») 

HES  ■ (£TSSL*V(3.1))’0.5*TSTEP/(inOO.*LHV) 

ETSC  . ETSC  ♦(V(1.1)+V(2.1)-'ETSCL)  • 0.5nSTEP/(10n0.*LHVl 
ETSS  • ETSS  ♦(V(3,l)*£TSSL)*0.5*TSTEP/(10nn.*LHV) 

SHSC  • SHSC  HSU  1)*S(2.1)*SHSCL)  • 0.5»TSTEP/(  1000.*LHV) 
SHSS  = SHSS  *(S(3.n*SHSSL)*0.5*TSTEP/(in00.*lHVl 
OG  ■ DG  ♦(G*SL)*0.5*TSTEP/(1000.‘LHV) 

ETSIN  = ETSIHt|V(l.l)tV(2.1)tV(3,l)t£TSmL)*n.5* 

♦ TST£P/(lflOn.*LH«) 

SKSIM  ■ SHSI>l*(S(l,ll*S(5,U*S(3,l)*SH51Hl.l*0.5* 

♦ TST£P/(1000.*LHV1 
a ■ e 

ETSia  ■ V(l,n*»(2,0*V(3,l) 

SHSIML  • S(l.l)*S(2,U*S(3.1) 

ETSCl.  = VU.U+HE.l) 

ETSSL  • V(3.1) 

SHSCL  ■ S(1,1)*S(2,1) 

SH55L  ■ S(3,l) 

TLAST  ■ TIHE(NHOUR) 

C**'  CONVERT  DIURNAL  VALUES  TO  MH/HR 
C 

SUNV  • SUNV*3SQ0./(LHV*1000.) 

SlPIH  • SUNH*3600./(UIV*1000.) 

r.  . G*36DO./(LHV*lflOO.) 

SUNSK  • (RK(l,ll^N(2.n*RN(3,l))*3600./(LHV*10fl0.) 

DO  100  IJ  • 1,3 

vriJ,l)  . V(IJ.1)*3600./(LHV*1000.) 

S(IJ,U  ■ 5(IJ.l)*3600./(LHV*100fl.) 

RN(IJ,1)  > RNriJ.l)'3600./(LHV*1000.) 

100  CONTINUE 

IF  (TS(3,1)  ,LT.  TS3MN)  TS3MN  • TS(3,I) 

IF  (TS(3,1)  ,GT.  TSW*)  TS3HX  • TS(3.1) 

IF  (TC  ,LT.  TACHN)  TACMN  » TC 

IF  (TC  .GT.  TACNXl  TACMX  - TC 

IF(TIN£(NH0UR1.GT.9.75.AN0.TW£(NH0UR1.LT.10.25) 

♦ TCTA-TC-TA(NNOUR) 

C 

C***  DIURNAL  OUTPUT 

RAONET  . RNIT  * RN2T  ♦ RN3T 
VPRES  • VPOF(TD(NHOUR)) 

IF  (LHflUR  ,L£.  0)  RETURN 

WRITE  (9,200)  OAy,OAPL,TIME(NHOUR).RN(l,l),V[l,l),S(l,l), 

» RN(2,1),V(2,1),S(2,11,RN(3.1),V(3.1).S(3,1),G,TC.TS(1,1). 
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‘ TS(2.1).TS(3,1).VPR£S.EC 
200  F0RHHT(2F5.n.U,F5.1,lX,10(fS.2  .IX) .6(F5.2,U1 ) 
RETURN 


”•  OAII.X  OUTPUT 

300  IF  (LDAILY  .LE.  0)  RETURN 
TSESN  ■ TSESN  t ETSC 
ESESN  • ESESN  » ETSS 
ETSESN  ■ ETSESNt  ETSC  ♦ ETSS 

WRITE  (10,400)  OAPL, LAI, ETSIM, ETSC. ETSS, ETSESN, TSESN, ESESN, 

» TCTA,TS3MX.TS3NN.DAPL 

400  F0RIIAT(F4.0,U  ,F3. 1 .U.3(F6.3,1X) ,3(F6.1,1X) .3(F5. 1 ,U) ,F4.0) 


SST0RAGE:2 

5NOFLOATCALL5 

SUBROUTINE  SOILTEHR 
SINCLUOE:  'CO'KION.OAT' 

REAL  EL0H(20).NFL0W(20) 

COND  ■ IABEL(THC0«O.THETA,SH(l).3) 

OELTOZ  > (TSURF-TENP(r)l/(0.5*0LAYR(l)) 


Q IDO 


DELTDZ 

FLOU(I) 

NFLM(I-1 

»HCAF(I- 


DELTDZ'COMD 
I • Z.NTLAYR 

= TAREL(THCOND.THETA.SU(M),3) 
' TABELITRCOND, THETA. SW(I)  ,3) 

■ OLAYR(M)tnLAYR(I))*0.5 

■ (TEMP(I-1)-TEUP[I))/THI« 

■ DELTDZ  • (C0ND1*C0ND2)*D.5 

■ (flom(i-i)-flo«(i))*tst£p 

■ VHCAP([-l)»NFLO»(I-l) 


F (VHCAPd-l)  . 

10  CONTINUE 

NFLOM(NTLAYR)  • FLOUINTLAYRl’TSTEP 
VNCAPfNTLAYRl  . YHCAPINTLAYKl  t NFLOH{NTLAYR ) 
DO  200  I » l.NTLAYR 
CAP  . TABEL(HCAP.Th£TA,SM(I1,3) 
TEMPd)  ■ VKCAP{I)/(OLAYRd)*CAP) 
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5ST0tAS£;2 

SNOFLOATCAUS 

SUBROUTINE  SOLUTION 
SINCLUDE:  'COWON.BAT' 

DIMENSION  A(3,3).N(3,3),RHC1(3.3),R«C2(3.3).ARV(3.31.AVFD(3,1), 
♦ARN(3.1).TSCl(3,l),TSC2r3.n.»VVPO(3.L).RVTSC2[3,l) 


C2 


SO  TO  1600 
100  CONTINUE 

CALL  GAUSSJ(RV,U.N,N,IEIIftORl 
IF  (lERROH  .LT.  0)  GO  TO  200 

GO  TO  1600 

200  00  300  I - l.N 

00  300  0 • l.N 
RHCUI.J)  ■ RH(I,J)  * Cl 
RVC2(I,J)  • RV(I,J)  ' 02 
300  A(I.J)  • RHCld.J)  ♦ R»C2(I, 

CALL  GAUSSJ(A,H,N.N. [ERROR) 

IF  (lERROR  .LT.  0)  GO  TO  400 

SO  TO  1600 

400  CALL  MATPRO(R\l,O.RVVPO.N.N.l) 

00  500  [ ■ l.N 
RVVFDd.U  ■ RVVR0d,ll»C3 
500  CONTINUE 

CALL  KATPRO(A,R»VPD.AVI>O.N,N.l) 


SO  TO  (( 

10  SUM  -0 

00  700  J > 
SUH  - SUM  ■ 
■n  CONTINUE 

call’matpro(a.rn 


N.J]  • RK{J.l) 

- ( SUM  1 t AVPO(N.l)  ) / A(N.N 
.N.N.!) 

■ ARNdJ.D-AKPOdJ.l)  + TA(NHOUfi) 


TS(lJ.l)  " ARNdJ.ll-AVPDdJ,l)tTA(NHOUR) 
■ TAREUTNCONO.TMETA.SVKD.S) 

• TEMP[1)  ♦ 0LAYR(l)*O.5-*S/C0ND 
» ABS(TS|3.l)-TSURF) 


nz 


IF  (CSROR  .LE.  ERRBNO)  SO  TO  1200 
OENOHl  ■ DLAYR(1)'0.6/C0ND 
RV3  . R«(3)  * R t RS(3) 

RR3  . RB(3)  ♦ R 

DEW  ■ (-GW‘iR/Cl)*R83*RV3yiRAWW*RV3*DELTA»RB3) 
OENOH  ■ OENOHl  - nEN0M2 
DELS  ■ (TS(3.1)-TSURF)/DEN0H 
GO  TO  1000 
1200  TSURF  • TS(3.1) 

RN(3,1)  • RN3 
1300  CONTINUE 

00  IROO  I = l.N 

TSCl(l.l)  ■ (ARNlI.l)  - AVPOd.D)  • Cl 
TSCEd.l)  ■ (ARN(I.l)  - AVPOd.D)  *02 
CONTINUE 

HATPR0(RH.TSC1,S,N,N,D 
HATPR0(RV,TSC2.RVTSC2.N.N.U 
MATA00(RVTSC2.RVVPD,V,N.l) 

• 0 

■ 0 

00  ISOO  I • l.N 
SUHH  • SUHH  Sd.l) 

SUHV  • sum  ♦ V(I,D 
CONTINUE 

(R’SIP1V/C3)  ♦ VPOF(TO(NHOURD 
(R*5UHH/CD  ♦ TA(NKOUR) 

IF  (KEY  .EO.  2)  50  TO  1700 
6 ■ RH3  - RN(3.D 
RN(3.D  ° RN3 

GO  TO  1700 

Ifion  WRITE  (*,•)  'ERROR  CODE  '.  KK 
STOP 

1700  CONTINUE 
RETURN 
END 
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SST09AGE;2 

ANOFLOATCALLS 

FUNCTION  T4BEM»AL.ARU,0UMHy,K) 

OINENSIOIt  VAL(K1,ARG(K) 

DO  100  J ■ 2. A 

IF  (OUHNY  .GT.  ARG(J11  GO  TO  100 
GO  TO  200 
100  CONTINUE 
J ■ K 

200  TABEL  • (OUMMY-ARGU-l)  )*(VAL(J)-VAL(J-1)  )/{ARG(J)- 
‘ ARG(J.l))tVAL(J-U 


SST0«ASE;2 

SNOFLOATCALU 

FUNCTION  VPDF(T) 

WPUCIT  S£ALM(A-H,0-Z) 

VPDF  ■ 2.7M9  • 1.E8  * EXP(-4278.6/(T  ♦ 242.791) 

RETURN 

END 


SSTQR«GE:2 

tNOFLOATCALLS 

SUBROUTINE  '.OHR 
SINCLUOE:  'COHiON.OAr 


100  REAO  (7. 200. END  • 400)  NYR, JULIAN, TINEI I) ,PRFD( I) ,RNO( I) 
♦ TA(n,T0(I).TSS(I),U(I) 

200  FORMAT!  13,1X,  13. 1».7(F8.2.H)) 

IF  (TIME(I)  .LT.  0)  00  TO  400 
IF  (I  .RT.  l)  00  TO  300 
HAST  ■ TIME(I) 

NHOKTH  ■ MONTH(OULIAN.NYR) 


2500.9  - 2.345*TA(I) 

TIKEin-TLAST 

TOTAL  . TOTAL  *(»N0( I]->RNLAST1*0.5*ITIH£1I )-TUST)*3 
RKUST  ■ RNO(I) 

TLAST  ■ TIME! I) 


400  CONTINUE 
NOAIAP  ■ I 
REIIINO  7 
RETURN 


FUNCTION  MONTH(JULIAN.NYR) 

OIMENSION  H0N(12) 

OATA  HOH/31, 29, 31, 30, 31, 30. 31 ,31,30,31,30.31/ 
M0N(2)  " 28 

IF  !(NYR/4)*4  .£0.  NYR)  H0N(2)  ■ 29 
NSUH  • 0 

00  100  J ■ 1.12 

IF  (JULIAN-NSIIM  .LE.  HON(J))  00  TO  200 
10  NSUN  = NSUM-HIDN(J) 
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COIMJN  BLOCS 


C0HHON/UETH/Tl»IE(5O),PPF0(50),2N0(5O),TA(50).TD(5O). 

2 TSS(50),  U(50) , TOTAL. JULtAN.NYR.NHOMTH.HDATAP 
C0MI«N/S0L/m0UR,R.RHl3.3).RV(3,3).RN(3.1).G,TS(3.1).l)(3.I). 

3 “J^2'^^SANP^M£LTA,RB{3).R3(3).S(3.1).V(3.1).EC.TC. 
C»W0K/STMP/TSTER.THC0IID(3)  ,TH£TA(3)  .HCA[>(3)  ,TEHP(2n)  .VHCAP(2D) 
C0m0N/IRR/SAT(20).0UL[20t.LL(2a),0LAVR[20),SU(20) 
CDHP«N/lNl/B8.£TSIH.ErSC,£TSS,SHSC.SHSS,SHSIM,TACH».TS3«N, 

2 TAC1X,TS3I1X 

C(V»lON/FIXI/DAPL,ORrOAY,FREONC .CHIGHT.LAIMAX.UAYG.RA.ERRBND. 

2 RtlAX.RNAmN.XLAT.XLDTC.KEY.EXTINC.FRACTR.CHAX. SLOPE, 

3 KRX.NELAYR.NTUYR.CEM.CEC.FTEZ.LOAILY.LHOUR 
COM'IOR/'IISC/TSESN.ESESN.CE.RNTOT.TLAST.LHV.OAY.LAI.RNA.ETSINL 

2 SHSWL.ETSCL.ETSSL.SHSCL.SHSSL.SL.KES.HEC 


FILE:  INPUT. OUT  - 

001  01  on  •*!  - 

0.50  0. 

3 8 18. 

n.l  0 


0.0180  -38.14E-6  1 


5 15  0.30 
I 0.3  0.3  U 


5 22.5  22.5  22.5  22.5  22.5  22.6  22.5 
0.10  O.IO  0.10  0.10  0.10  0.10  0.10  0.10 
0.25  0.25  0.25  0.25  0.25  0.25  0.25  0.25 
0.10  0.10  0.10  0.10  0.10  0.10  0.10  0.10 
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